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ABSTRACT 


■^This  paper  describes  an  analytic  model  of  the  strategic  bomber 
penetration  mission.  The  model  covers  all  parts  of  the  mission  from 
take  off  to  kill  or  recovery.  Specifically,  the  model  includes  repre¬ 
sentations  of  base  escape,  aerial  refueling,  forward  defense,  barrier 
SAMs,  area  fighter-intercepter  defense,  random  area  SAMs,  terminal 
SAMs,  weapons  delivered,  and  recovery.  Cruise  missiles  are  included 
i n  the  model ^ 

K 

For  each  part  of  the  mission  listed  above,  a  distribution  is 

\ 

derived  for  the  number  of  bombers  that  survive  given  the  number  that 
begin  that  stage  of  the  mission.  A  computer  program  has  been  written 
to  convolute  these  distributions  to  find  distributions  for  the  total 
number  of  bombers  that  survive,  the  total  number  of  weapons  released, 
and  the  total  number  of  cruise  missiles  that  destroy  targets. 

One  major  advantage  of  this  model  is  its  scope.  It  includes 
descriptions  of  all  the  events  that  occur  during  the  mission  which 
will  have  an  impact  on  the  outcomes  of  later  parts  of  the  mission. 

A  second  major  advantage  of  this  model  is  the  inclusion  of  calcu¬ 
lations  for  the  variances  of  the  results  of  the  model.  This  is  an 
improvement  on  other  analytic  models  which  give  only  expected  value 
results  and  hence  contain  no  indication  of  the  possible  variability 
of  those  results. 

All  the  mathematical  derivation  for  this  model  are  contained  in 
this  paper.  The  paper  also  includes  some  initial  results  from  the  com 
puter  program  which  demonstrate  some  of  the  many  outputs  that  can  be 
obtained  from  the  model. 
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CHAPTER  I 


INTRODUCTION 

Models  of  the  strategic  bomber  mission  have  always  been  of  great 
interest  to  the  Air  Force.  In  recent  years  this  interest  has  been 
heightened  by  contraversies ,  often  of  a  political  nature,  that  have 
developed.  Much  debate  has  centered  on  the  need  to  replace  the  Air 
Force's  aging  fleet  of  B - 52  bombers  with  a  more  modern  weapon  such  as 
the  B-l  or  a  derivative  of  it.  The  problems  of  a  bomber  fleet's  sur¬ 
vivability  and  its  ability  to  penetrate  a  defended  country  raise  the 
question  of  whether  the  manned  bomber  is  still  a  viable  part  of  the 
strategic  deterrent  when  compared  to  such  alternatives  as  cruise 
missiles  or  land  and  sea  based  ballistic  missiles.  If  one  concludes 
that  a  manned  bomber  can  still  be  an  effective  strategic  weapon,  the 
question  remains  whether  a  new  bomber  could  justify  its  cost  by  per¬ 
forming  the  mission  significantly  better  than  the  B-52. 

To  help  answer  these  and  many  other  questions  about  the  strategic 
bomber  mission,  numerous  models  have  been  developed  which  describe  all 
or  part  of  the  mission.  These  models  fall  into  two  broad  categories, 
computer  simulation  models  and  analytic  models.  In  general,  computer 
simulations  can  offer  greater  levels  of  detail  and  potentially  greater 
flexibility;  however,  large  simulations  can  be  both  costly  and  time 
consuming  to  run.  Analytic  models  are  usually  cheaper  and  faster,  but 
they  generally  give  only  expected  value  results  which  may  not  contain 
the  same  amount  of  information  as  a  large  number  of  runs  of  a  simula¬ 
tion  model.  In  particular,  the  analytic  models  seldom  give  any  indi- 
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cation  of  the  possible  uncertainty  (variability)  of  their  results. 

This  is  the  issue  of  primary  concern  in  this  paper. 

The  variance,  or  alternatively,  the  standard  deivation,  of  a 
result  from  a  model  is  an  important  measure.  It  gives  an  analyst 
some  idea  of  the  range  of  values  being  represented  in  his  results. 

This  in  turn  allows  him  to  better  evaluate  the  degree  of  precision 
contained  in  the  results.  This  can  help  an  analyst  avoid  overempha¬ 
sizing  the  importance  or  misinterpreting  the  meaning  of  the  results  of 
his  model.  It  can  in  some  circumstances  be  an  important  factor  in  a 
decision  involving  a  choice  between  alternatives. 

As  an  example  of  how  the  variance  of  a  result  could  affect  a 
decision,  consider  the  following  situation.  A  small  airstrike  is  being 
planned  against  a  defended  target.  Two  alternatives  are  available. 

Plan  A  is  to  send  two  aircraft  against  the  target,  each  carrying  one 
weapon.  These  aircraft  are  estimated  to  each  have  a  60-40  chance  of 
penetrating  the  defense  and  reaching  the  target.  Plan  8  involves  a 
single  newer  aircraft  having  twice  the  payload  of  the  older  aircraft 
and  an  estimated  60-40  chance  of  surviving  the  defense  and  reaching  the 
target. 

In  Plan  A,  if  the  two  aircraft  are  independent  of  each  other,  the 
distribution  for  the  number  of  weapons  that  reach  the  target  is 

P(2  weapons  reach  target)  =  0.36 

P(1  weapon  reaches  target)  =0.48 

P ( 0  weapons  reach  target)  =  0.16 

Under  Plan  B,  the  single  aircraft  will  either  reach  the  target 
and  release  all  its  weapons  or  fail  to  reach  the  target  at  all.  The 
distribution  for  the  number  of  weapons  reaching  the  target  is 
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P( 2  weapons  reach  target)  =  0.6 

P ( 0  weapons  reach  target)  =  0.4 

The  expected  number  of  weapons  that  would  reach  the  target  under 
either  plan  is  1.2.  If  this  is  the  only  information  aviable  to  the 
decision  maker,  he  will  likely  conclude  that  both  plans  are  equally 
effective. 

Now  look  at  the  variances  of  the  two  results.  Plan  A  has  a 
variance  of  0.48  while  Plan  B  has  a  variance  of  0.96.  If  the  decision 
maker's  primary  concern  is  to  have  at  one  weapon  reach  the  tar¬ 

get,  then  knowing  the  relative  variances  involved  in  the  two  plans  he 
will  likely  choose  Plan  A  instead  of  Plan  B.  Even  though  the  expected 
total  number  of  weapons  reaching  the  target  is  the  same  in  both  plans, 
the  smaller  variance  of  Plan  A  indicates  that  it  may  actually  have  a 
better  chance  of  delivering  at  least  one  weapon  to  the  target.  That 
this  is  indeed  the  fact  can  be  verified  from  the  distributions. 

This  is  admittedly  a  very  simple  example  but  it  serves  nonetheless 
to  illustrate  the  importance  of  knowing  the  variance  in  model  results. 
In  a  simulation  model,  an  estimate  of  the  variance  of  the  results  can 
be  obtained  by  running  the  model  several  times  and  calculating  the 
sample  variance  of  the  results.  An  analytic  model,  however,  generally 
gives  only  expected  value  results  and  no  variance  information  obtained. 
This  is  the  case  in  virtually  all  of  the  analytic  models  currently 
existing  that  deal  with  the  strategic  bomber  penetration  problem. 

The  original  suggestion  for  a  thesis  in  the  area  of  strategic 
bomber  penetration  came  from  Captain  R.  Wilkinson  at  the  Strategic 
Systems  Program  Office  of  the  Aerospace  Systems  Division  (Strat  SPO). 
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Of  current  interest  to  the  Strat  SPO  is  an  expected  value  analytic 
model  of  the  bomber  penetration  problem  developed  by  William  J.  Schultis 
at  the  Institute  for  Defense  Analysis  (Ref  13).  Captain  Wilkinson 
suggested  that  the  Schultis  model,  which  currently  does  not  cover  the 
problems  of  base  escape,  aerial  refueling,  and  recovery,  be  expanded  to 
include  the  entire  strategic  bomber  mission.  Because  the  Schultis 
model  is  an  expected  value  model  with  no  variance  calculations,  it  was 
decided  to  broaden  the  scope  of  the  project  to  include  the  variance 
problem.  In  doing  this,  some  parts  of  the  Schultis  approach  were 
employed;  however,  some  concepts  were  added  from  other  models  and  more 
were  developed  specifically  for  this  model.  The  Schultis  model,  along 
with  other  models  in  the  area  of  strategic  bomber  penetration,  is 
discussed  in  more  detail  in  the  next  chapter. 

OBJECTIVES  AND  SCOPE 

The  purpose  of  this  paper  is  to  develop  an  analytic,  national 
level  model  of  the  strategic  bomber  mission  from  take-off  to  kill  or 
recovery.  The  object  is  to  construct  the  model  to  allow  for  both 
expected  value  and  variance  calculations. 

Structure.  The  strategic  bomber  mission  can  be  divided  up  into 
nine  parts,  not  all  mutually  exclusive,  which  define  the  major  aspects 
of  the  mission.  These  are: 

a)  Base  Escape  -  bombers  escape  from  their  own  bases  which  may  be 
under  SLBM  attack; 

b)  Aerial  Refueling  -  bombers  are  refueled  in  flight  during  the 
mission; 

c)  Forward  Air  Defense  -  bombers  are  attacked  by  defensive  inter¬ 
ceptors  before  reaching  the  defended  continent; 
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d)  Barrier  SAMs  -  bombers  are  attacked  by  a  perimeter  SAM  defense 
as  they  reach  the  coast  of  the  defended  continent; 

e)  Random  Area  SAMs  -  bombers  penetrate  to  targets  inside  an  area 
defended  by  randomly  located  SAM  sites; 

f)  Fighter-Interceptors  -  bombers  are  attacked  by  defensive 
fighter-interceptors  while  penetrating  to  their  targets; 

g)  Weapons  Del i vered/Termi nal  SAMs  -  bombers  deliver  weapons  at 
targets  which  may  be  defended  by  terminal  SAMs; 

h)  Recovery  -  bombers  leave  the  defended  area  and  attempt  to 
find  a  friendly  base  at  which  to  recover; 

i)  Cruise  Missiles  (CMs)  -  during  the  mission,  bombers  may 
release  CMs  which  will  penetrate  to  targets  on  their  own. 

Philosophy.  Because  of  the  size  and  complexity  of  this  model  of 
the  bomber  penetration  problem,  the  philosophy  of  this  model  is  to 
sacrifice  resolution  where  necessary  in  favor  of  breadth.  In  the 
final  chapter  of  this  paper  some  areas  are  suggested  where  it  should 
be  possible  to  build  further  detail  into  the  model.  However,  it  should 
be  recognized  that  there  are  practical  limitations  to  the  level  of 
detail  that  can  be  obtained  in  an  analytic  model  of  this  scope.  To 
achieve  finer  detail,  it  will  likely  be  necessary  to  return  to  large 
simulation  models  such  as  the  Boeing  Company's  Advanced  Peneteration 
Model  (see  Chapter  II) . 

APPROACH 

The  approach  used  to  model  the  strategic  bomber  mission  is  to 
first  separate  the  mission  into  parts  that  can  be  modelled  individually 
and  then  to  integrate  the  various  parts  into  a  complete  model.  In 
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terms  of  the  parts  of  the  mission  previously  defined,  it  will  be 
assumed  that  the  first  four  parts,  base  escape,  aerial  refueling,  for¬ 
ward  air  defense,  and  barrier  SAMs,  occur  sequentially  and  do  not  over¬ 
lap.  Hence,  they  can  each  be  modelled  separately.  Random  Area  SAM 
and  f igher-interceptor  defenses,  however,  occur  simultaneously.  In 
this  model,  they  are  assumed  to  operate  independently  of  each  other  and 
so  are  initially  considered  separately.  The  results  of  each  are  then 
combined  in  the  model  under  area  defenses.  The  terminal  SAM  and  weapons 
delivered  part  of  the  mission  would  in  an  actual  mission  overlap  the 
area  defense  stage;  however,  in  this  model  it  will  be  treated  separately. 
Recovery  is  the  last  part  of  the  mission  and  will  be  treated  separately 
in  the  model.  The  release  of  curise  missiles  will  be  considered  at 
three  different  points  in  the  mission;  first  between  the  refueling  and 
FAD  stages;  second,  between  the  FAD  and  BSAM  stages;  and  finally, 
during  the  area  defense  penetration  stage. 

A  flow  chart  of  the  model  is  shown  in  Figure  1.1.  Each  part  of 
the  model  will  be  developed  independently  of  the  others  and  the  out¬ 
puts  from  early  parts  will  be  the  inputs  to  later  parts.  An  outline 
of  the  various  parts  of  the  model  is  given  below. 

Base  Escape.  The  initial  number  of  bombers  in  the  problem  is  an 
input  and  the  model  finds  a  distribution  for  the  number  of  bombers 
that  survive  the  SLBM  attack  on  their  bases. 

Aerial  Refueling.  The  initial  number  of  tankers  is  an  input  and 
the  model  finds  a  distribution  of  the  number  of  bombers  that  success¬ 
fully  refuel  for  a  given  number  that  survived  the  base  escape. 

Forward  Air  Defense.  The  number  of  defensive  interceptors  is  an 
input  and  the  model  finds  distributions  for  the  number  of  bombers 
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that  survive  the  FAD  given  the  number  that  escaped  their  bases  and 
the  number  that  were  refueled. 

Barrier  SAMs.  The  number  of  SAM  sites  in  the  perimeter  defense 
and  the  width  of  the  perimeter  are  input  and  the  model  gives  distri¬ 
butions  for  the  number  of  bombers  that  survive  the  5SAM  defense  given 
the  number  that  survived  the  FAD. 

Area  Defense.  The  number  of  SAM  sites  in  the  defended  area  and 
the  size  of  the  defended  area  are  input.  The  model  finds  a  single 
bomber's  probability  of  penetrating  through  the  RASAM  defenses  as  a 
function  of  the  depth  of  penetration. 

An  intercept  intensity  parameter  describing  the  average  number  of 
intercepts  that  could  occur  on  a  detected  bomber  per  unit  time  is  in¬ 
put  and  the  model  finds  a  bomber's  probability  of  penetrating  through 
area  FI  defenses  as  a  function  of  the  depth  of  penetration. 

The  two  penetration  probabilities  are  then  combined  to  give  a 
single  probability  that  a  bomber  survives  the  area  defenses  as  a  func¬ 
tion  of  the  depth  of  penetration. 

Terminal  SAMs  and  Weapons  Delivered.  The  bombers'  target  set  and 
the  level  of  terminal  defenses  at  the  targets  are  input  and  the  model 
gives  a  distribution  for  the  total  number  of  weapons  delivered  by  a 
given  number  of  bombers  that  have  survived  the  BSAM  defense. 

Recovery.  The  model  gives  a  distribution  for  the  number  of  bombers 
that  survive  to  exit  the  defended  area  for  a  given  number  of  bombers 
that  have  survived  the  BSAM  defense. 

Cruise  Missiles.  The  number  of  CMs  on  each  bomber  and  a  release 
policy  are  input  and  the  model  gives  a  distribution  for  the  total  num- 
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ber  of  CMs  that  reach  their  targets.  Cruise  missiles  may  be  released 
at  any  of  three  points  in  the  model. 

In  this  model,  the  defended  area  in  which  the  targets  lie  is 
assumed  to  be  a  rectangular  corridor.  The  BSAM  defenses  lie  at  the 
northern  boundary  of  the  corridor  and  bombers  enter  from  the  north  and 
penetrate  to  their  targets  along  straight  lines  to  the  south.  Once 
the  bombers  have  released  their  weapons  they  do  not  attempt  to  return 
to  the  northern  boundary  but  continue  to  fly  south  until  they  exit  the 
corridor  at  the  southern  end. 

FORMAT 

Chapter  II  of  this  paper  gives  a  brief  review  of  some  of  the 
currently  existing  models  that  deal  with  the  strategic  bomber  mission. 

Chapter  III  presents  a  more  detailed  description  of  the  stages  of 
the  bomber  mission  to  be  modelled.  The  major  assumptions  made  in  the 
model  are  discussed  in  this  chapter. 

Chapter  IV  develops  the  model  representation  of  the  first  four 
stages  of  the  strategic  bomber  mission.  These  are  ba.  »  escape,  aerial 
refueling,  forward  air  defense,  and  the  barrier  SAM  defense. 

Chapter  V  develops  model  representations  for  the  two  forms  of  area 
defense,  random  area  SAMs  and  fi ghter- interceptors . 

Chapter  VI  develops  model  representations  for  the  final  parts  of 
the  strategic  bomber  mission,  terminal  SAMs  and  weapons  delivered, 
and  recovery.  This  chapter  also  includes  the  modelling  of  cruise 
missiles. 

Chapter  VII  shows  how  the  various  parts  of  the  model  that  have 
been  developed  in  the  previous  chapters  are  integrated  to  form  a 
complete  model.  The  chapter  includes  a  brief  discussion  of  the 
computerization  of  the  model. 


Chapter  VIII  briefly  discusses  verification  of  the  model  and 
shows  some  of  the  model  results. 

Chapter  IX  presents  a  summary  of  this  paper  along  with  some  con¬ 
clusions  and  recommendations. 


FIGURE  1.1 

Flowchart  of  the  model. 
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CHAPTER  II 


REVIEW  OF  CURRENT  MODELS 

The  purpose  of  this  chapter  is  to  briefly  review  some  of  the  cur¬ 
rently  existing  models,  both  analytic  and  simulation,  that  deal  with 
all  or  parts  of  the  strategic  bomber  penetration  mission.  Some  of 
the  models  mentioned  below  are  force  on  force  models  that  treat  missile 
exchange  in  addition  to  bomber  penetration;  they  are  included  because 
they  help  to  put  in  perspective  the  penetrating  bomber's  role  in  the 
strategic  triad. 

Simulations 

Advanced  Penetration  Model.  Perhaps  the  most  comprehensive  pene¬ 
tration  model  developed  to  date  is  the  Advanced  Penetration  Model  (APM), 
developed  by  Boeing  for  Headquarters  USAF,  the  Assistant  Chief  of  Staff, 
Studies  and  Analysis  (Ref  11).  The  model  simulates  the  strategic 
mission  of  the  entire  bomber  and  tanker  force  from  takeoff  to  recovery 
(Ref  7 :  IV-4) . 

The  model  consists  of  two  main  parts,  a  Mission  Planner  and  an 
Air  Battle  Simulator.  The  overall  mission  scenario  is  user  defined; 
the  model  will  then  generate  individual  flight  plans  for  each  bomber 
(and  tanker)  in  the  force.  Various  rules  or  constraints  may  be 
imposed  on  the  Mission  Planner.  The  plan  for  each  sortie  will  include 
routing,  refueling,  target  allocation,  and  recovery  (Ref  7 : I V-8) . 

The  Air  Battle  Simulation  is  a  time  sequenced  process  of  the 
events  that  have  been  planned  In  the  Mission  Planner  and  include  both 
deterministic  and  probabal istic  events.  Seven  different  processors 
are  used  to  model  the  interaction  of  bombers  and  defenses  over  the 
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enemy's  defended  airspace.  These  are 

a)  Command  and  Control  -  assigns  interceptors  to  penetrators; 

b)  Radar  -  detects  penetrators;  includes  effects  of  ECM  use; 

c)  AWACS; 

d)  Interceptor  -  vectors  assigned  interceptor  to  radar  range 
of  the  penetrator,  finds  probability  of  detection  and  conversion  and 
proability  of  kill  and  determines  if  the  penetrator  is  killed; 

e)  SAM  -  finds  probability  of  kill  on  a  penetrator  for  a  SAM 
and  determines  if  the  penetrator  is  killed; 

f)  Nuclear  Effects  -  simulates  CEP  of  a  weapon  and  calculates 
the  damage  done  by  the  weapon;  and 

g)  Penetrator  (Ref  7 : I V-20) . 

All  the  results  of  the  air  battle  simulations  are  stored  as  out¬ 
put,  allowing  a  vast  amount  of  specific  information  to  be  obtained 
from  the  model.  The  user  may  specify  the  particular  parts  of  this 
information  he  wishes  to  analyze. 

This  wealth  of  output  data  and  the  relatively  fine  detail  of  the 
APM  are  its  major  advantages.  The  model  can  capture  the  complex 
realities  of  defense  saturation,  command  and  control  limitations  and 
weapon  assignment  policies  (Ref  7 : 1 V-4 ) .  At  the  same  time,  the  vast 
amount  of  data  that  is  inputted  to  the  model  offers  great  flexibility 
to  the  user. 

However,  in  the  ten-year  history  of  the  APM,  this  same  level  of 
detail  and  data  handling  has  also  been  the  model's  single  biggest 
drawback.  Single  replications  of  the  model  for  any  reasonably  large 
force  could  require  over  twenty  hours  of  computer  time  (Ref  7-.IV-26). 
Extensive  studies  including  repeated  runs  for  sensitivity  tests  are 


11 


both  time  consuming  and  expensive.  New  generation  computers,  such 

as  the  IBM-3032  currently  being  used  for  the  APM,  have  dramatically 

reduced  the  model's  run  time;  however,  the  inputting  and  maintenance 

of  the  computer  routines  and  data  base  remain  costly  and  time  consuming 

and  prohibit  widespread  practical  use  of  the  APM  (Ref  7 : 1 V-26 ) . 

SPEED.  SPEED  (Simulation  of  Penetrator  Encountering  Extensive 

Defenses)  was  developed  by  Calspan  Corporation  and  is  a  large  Monte 

Carlo  simulation  of  bomber  penetration  through  air  defense  systems. 

2 

SPEED  simulates  fighters,  SAMs,  C  ,  ECM,  and  ground  controlled  inter¬ 
cepts  (GCI)  (Ref  8:13).  SPEED  is  compatible  with  and  complements  the 
APM;  in  particular,  the  SPEED  model  can  use  the  APM  Mission  Planner 
(Ref  1:3-12,  8:13) . 

FISCHER.  FISCHER  is  an  event  oriented  Monte  Carlo  simulation  of 
an  engagement  between  an  attacking  bomber  force  and  a  strategic  defense 
system.  The  model  was  developed  at  Headquarters  N'ORAD  (Ref  l:B-6). 

The  model  simulates  the  movements  of  bombers  and  air-to-surface 
weapons  and  their  interaction  with  a  defense  system  including  both 
ground  and  airborne  radars  and  interceptors  which  can  be  committed 
from  either  bases  on  CAP  orbits.  The  model  permits  fuel  monitoring, 
with  reattack  and  recommit  logic  available  (Ref  l:B-4).  The  original 
model  has  been  modified  for  use  by  AF/SAS  (Ref  l:B-4). 

ANALYTIC  MODELS 

ARSENAL  EXCHANGE  MODEL.  The  Arsenal  Exchange  Model  was  developed 
by  the  Martin-Marietta  Corporation.  It  is  a  large  expected  value 
model  used  to  study  the  structure  of  total  strategic  forces  (Ref  8:13). 
This  model  includes  a  perimeter  area  defense  with  a  probability  of 
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bomber  penetration  P^g  given  by 


ppB  •  u-fE>  *  p£n-pAI)- 


where 

Pg  -  probability  of  encountering  the  parameter  defense 

P^j  -  probability  that  a  bomber  is  killed  by  single  interceptor 

pass 

I,B  -  total  number  of  interceptors  and  bombers  respectively 
(Ref  11:5). 

Several  different  techniques  are  used  to  optimize  weapon  alloca¬ 
tion  and  the  model  can  treat  full  force  allocations,  a  variety  of 
defenses  and  force  design  problems.  A  number  of  scenarios  ranging 
from  a  single  strike  against  military  and  value  targets  to  three  strike 
games  involving  problems  of  selecting  a  weapons  reserve  or  a  value 
target  reserve  for  the  initiaters  third  strike  may  be  analyzed  using 
the  model  (Ref  l:b-4).  Three  counterforce  exchange  options  are  also 
available  including  limited  one-  and  two-strike  models  and  an  all  out 
counterforce  model  (Ref  1:8-4).  The  model  is  used  at  Headquarters  USAF, 
Studies  and  Analysis  (AF/SA). 

CODE-50.  This  is  a  widely  used  aggregated  model  developed  by  the 
Lambda  Corporation  (Ref  11:4).  It  is  capable  of  handling  a  mixture 
of  offensive  weapons  types.  Bomber  penetration  probabilities  which 
describe  a  bomber's  chance  of  surviving  against  a  specific  type  of 
fighter  can  be  user  input  using  data  from  other  models.  Alternatively, 
for  B  bombers  and  F  fighters,  the  bomber  penetration  probability  can 
be  assumed  to  be  proportional  to  e'a^F^  where  a  and  c  are  model 
inputs  that  describe  each  bomber/fighter  pair  and  incorporate  the 
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effects  of  most  of  the  parameters  effecting  bomber  penetration  (Ref  11:4). 
The  model  is  probably  oversimplified  in  that  it  is  difficult  to  find 
appropriate  values  for  a  and  c  that  will  adequately  represent  penetra¬ 
tion  probabilities  generated  by  more  detailed  simulations  or  extracted 
from  actual  tactical  exercises  (Ref  11:5). 

COLLIDE.  COLLIDE  is  an  aggregated  conversion  model  for  air  combat 
designed  to  assess  the  impact  of  command  and  control  on  fighter-bomber 
engagements.  The  principle  output  of  the  model  is  a  probability  for 
interceptor  detection  and  conversion  under  different  engagement  scenarios 
(Ref  1 : B- 5 ) .  The  COLLIDE  model  incorporates  bomber  electronic  counter 
measures  (ECM)  in  its  calculations.  The  COLLIDE  model  is  used  to  esti¬ 
mate  ECM  effectiveness  for  the  APM.  Output  from  COLLIDE  is  used  as 
input  for  the  APM  (Ref  6.-IV-27). 

COPEM-1  (Corridor  Penetration  Model).  COPEM-1  was  developed  at 
Stanford  Research  Institute  (SRI)  as  part  of  a  study  to  improve  the 
representation  of  airborne  strategic  systems  in  aggregated  effective¬ 
ness  evaluation  models  (Ref  12:v).  The  other  part  of  the  study  was  the 
development  of  an  allocation  model  to  use  the  output  of  COPEM-1. 

COPEM-1  is  a  time  dependent  engagement  model.  Its  purpose  is  to 
generate  average  bomber  penetration  probabilities  as  a  function  of  the 
depth  of  penetration  along  a  single  attack  corridor  into  a  defended 
area.  The  model  includes  only  fighter  area  defenses;  SAMs  are  handled 
separately  in  the  allocation  model  (Ref  5:5).  The  underlying  assump¬ 
tion  in  the  model  is  that  the  number  of  intercepts  that  occur  on  a 
bomber  follow  a  Poisson  distribution  with  a  time  dependent  parameter 
(Ref  5:6).  This  parameter  is  calculated  iteratively  at  discrete  time 
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intervals  during  the  engagement.  It  incorporates  all  the  information 
about  the  defensive  fighter  except  for  the  single  pass  kill  probability 
for  a  fighter  against  a  bomber  which  is  a  separate  input  (Ref  5:6).  The 
basic  concepts  of  the  COPEM-1  model  have  been  employed  in  this  paper 
to  model  bomber  penetration  of  area  fighter  defenses.  The  mathematical 
development  is  presented  in  Chapter  V.  Later  versions  of  this  model 
have  gone  by  the  name  BOMPEN;  the  names  are  frequently  used  inter¬ 
changeably  (Ref  11:20). 

HISTVEC.  This  is  a  fast  running  expected  value  model  of  bomber 
penetration.  Fighters  and  bomber  decoys  are  both  modelled  in  detail 
(Ref  8:13).  In  particular,  the  model  considers  fighter  airbase  deploy¬ 
ment,  different  fighter  types  with  different  AAMs  possible  and  detec¬ 
tion  and  conversion  probabilities  that  are  functions  of  both  fighter 
type  and  penetration  altitude  (Ref  8:15).  Decoy  considerations  such 
as  flight  range,  credibility,  threat  dilution,  and  primary  payload  dis¬ 
placement  are  all  incorporated  in  the  model  (Ref  8:16).  The  model 
includes  no  representations  of  BDMs  or  SAM  defenses  and  does  not  model 
ECM  (Ref  8:16,19) . 

LULEJIAN-MARKOVIAN.  A  model  developed  by  Lulejian  and  Associates 
uses  a  Markov  process  to  model  the  bomber  penetration  problem  (Ref 
8:13).  Both  fighters  and  decoys  are  modelled  in  as  great  detail 
as  in  HISTVEC  with  the  added  feature  that  fighters  can  be  reassigned 
while  airborne  (Ref  8:15).  SAM  systems  are  modelled  with  two 
two  different  SAM  types  allowed.  Firing  rate  limits  and  degradation 
by  chaff  or  ECM  are  also  included  in  the  model's  representation  of 
SAMs.  BDMs  are  modelled  in  detail  including  the  option  of  short-  or 
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long-range  BDMs.  The  target  structure  can  also  be  modelled.  The 
effects  of  a  precursor  IC3M/SLBM  strike  can  be  included  in  the  model 
(Ref  3:15-16) . 

PENETRATION  INTEGRALS.  The  penetration  integrals  were  formulated 
by  Carl  Builder  at  the  RAND  Corporation  (Ref  15:1).  The  basic  approach 
is  to  use  differential  equations  to  model  combat  based  on  random 
encounters  between  interceptors  and  penetrators.  3uilder's  equations 
give  simple  analytic  results  for  the  extreme  cases  where  interceptors 
have  unlimited  weapons  load  or  unlimited  attack  time  available  (Ref 
15:1).  Builder's  penetration  integral  unifies  these  two  extremes  and 
gives  a  general  solution  for  the  cases  where  there  are  restrictions  on 
weapons  loads  and/or  attack  time  (Ref  15:1,  4 : i i i ) . 

PENEX.  PENEX  models  strategic  bomber  penetration  against  manned 
interceptors  (SAM  defenses  are  not  included).  The  basic  assumption 
made  is  that  the  air  battle  can  be  divided  up  into  discrete  "sub-battles" 
which  occur  sequentially  in  time.  Penetrators  and  decoys  that  survive 
one  sub-battle  will  progress  to  the  next.  At  each  new  sub-battle,  fresh 
interceptors  may  be  added  to  those  that  have  capability  remaining  after 
previous  sub-battles.  Both  close  controlled  and  raid  controlled  inter¬ 
ceptors  are  modelled.  The  basic  model  has  been  extended  to  consider 
several  types  of  interceptors  and  to  consider  the  effect  of  command 
and  control  system  confusion  resulting  in  false  tracks  or  failures  to 
track  some  penetrators.  The  model  was  developed  at  AF/SA  (Ref  2). 

THE  SCHULTIS  MODEL.  Next  a  model  which  had  a  substantial  impact  on 
this  thesis  is  discussed.  The  Schultis  model,  titled  "A  National -Level 
Analytic  Model  for  Penetration  of  Various  Combined  Air  Defense  Deploy¬ 
ments  by  Cruise  Missiles  or  Bombers,"  is  a  small  (capable  of  being  run 
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on  a  TI-59)  expected  value  penetration  model  (Ref  13).  Five  types  of 
air  defenses  are  modelled;  forward  air  defense,  barrier  SAMs,  random 
area  SAMs,  fighter-interceptors,  and  terminal  SAMs.  The  model  uses 
value  extracted  as  its  measure  of  effectiveness  (Ref  13:2). 

The  basic  approach  taken  in  the  model  is  to  separate  the  defenses 
into  bands  that  are  penetrated  sequentially  by  the  bombers.  In  this 
approach,  the  FI  defense  is  envisioned  as  a  band  which  lies  inside  the 
barrier  SAMs  but  does  not  actually  include  the  heartland  where  the 
targets  lie.  Random  area  SAMs  are  interspersed  with  the  targets  but 
do  not  interfere  with  the  FIs. 

The  model  deals  with  large  numbers  of  penetrators,  particularly 
CMs ,  and  relies  on  saturation  of  the  defenses,  rather  than  leakage,  as 
the  primary  method  of  penetration.  For  this  reason,  it  is  assumed  that 
the  offensive  penetrators'  best  strategy  is  to  attack  in  files  along 
narrow  corridors  Instead  of  individually  at  random.  As  the  files 
approach  each  line  of  defense,  losses  may  be  expected  at  first;  however, 
the  SAMs  will  exhaust  their  missiles  or  the  fighters  within  range  of  a 
particular  file  will  exhaust  their  AAMs  and  in  effect  a  path  through 
any  particular  layer  of  defense  will  be  cleared  for  the  remaining 
penetrators. 

Terminal  SAMs  are  treated  as  a  fixed  price  per  unit  value  that  must 
be  paid  by  the  penetrators  to  destroy  a  target  of  given  value.  Using 
this  approach,  the  model  can  find  the  total  value  extracted  by  a  given 
number  of  penetrators  that  have  survived  the  bands  of  defenses  for 
various  allocation  schemes  (Ref  13). 

As  mentioned  in  the  last  chapter,  it  was  the  Schultis  model  that 
originally  led  to  this  project.  In  this  paper,  the  same  five  basic 
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air  defenses  are  modelled  and  to  some  extent,  the  concept  of  layers  of 
defenses  is  employed.  Unlike  the  Schultis  model,  however,  this  thesis 
does  not  treat  fighter-interceptors  as  a  band  defense  like  the  barrier 
SAMs.  In  the  Schultis  model,  the  two  defenses  are  treated  exactly  the 
same  way  to  the  point  that  airbases  are  effectively  modelled  as  SAM 
sites  with  wide  coverages  and  individual  FIs  are  modelled  as  SAMs  which 
can  kill  more  than  one  target. 

The  Schultis  approach  to  modelling  one-on-one  penetrator/SAM 
engagements  is  used  in  this  thesis;  the  mathematics  are  presented  in 
Chapter  IV  in  the  secion  on  barrier  SAMs.  The  concept  of  files,  how¬ 
ever,  is  replaced  by  a  random  attack  and  saturation  of  the  defenses  is 
not  considered. 

OTHER  MODELS 

The  above  are  some  of  the  more  important  and  well-known  models  that 
have  been  developed  in  the  area  of  strategic  bomber  penetration.  The 
following  is  a  partial  listing  of  other  penetration  models  that  may  be 
of  interest. 

a)  Beta  Cadens  -  a  very  large  simulation  model  that  includes  all 
strategic  forces  and  detailed  damage  assessment  information  (Ref  8:13). 

b)  NYLAND  (RAND)  -  a  small  expected  value  model  that  includes 
bombers,  decoys,  BDMs,  interceptors  and  SAMs  (Ref  8:13-16). 

c)  OMEGA  (SRI)  -  an  allocation  model  that  optimizes  bomber  force 
allocations  against  targets  that  may  be  locally  defended  by  SAMs  under 
a  fighter  umbrella  (Ref  11:2). 

d)  PEGASOS  (Academy  for  Interscience  Methodology)  -  an  expected 
value  model  which  includes  FI  and  SAM  defenses,  bomber  decoys,  and 
ICBMs.  Does  not  cover  BDMs  (Ref  8:13-16). 
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e)  QUICK  (Lambda  Corp.)  -  a  large  simulation  that  generates  war- 
plans  and  will  optimize  weapon  allocations  within  restraints  (Ref  11:7). 

f)  STRAT  CRCW  -  a  computer  simulation  to  model  the  effects  of 
various  types  of  ECM  measures  (Ref  1 : B - 1 5 ) . 

g)  STRAT  SPLASH  -  a  model  used  to  estimate  single  shot  kill 
probabilities  for  AAMs  (Ref  1:B~22). 

h)  STRAT  SURVIVOR  -  a  detailed  simulation  model  of  the  base 
escape  problem  (Ref  1 : B - 2 3 ) . 

i)  UNCLE  (RAND)  -  a  model  specifically  designed  to  account  for  the 
large  number  of  uncertainties  the  offense  and  defense  may  have  about 
each  other  (Ref  11 :6) . 

This  is  by  no  means  a  complete  list  of  the  models  that  currently 
exist  to  describe  the  strategic  bomber  penetration  mission.  It  does, 
however,  cover  most  of  the  more  important  models.  Some  of  these 
models,  such  as  COLLIDE  and  STRAT  SPLASH,  may  be  useful  when  choosing 
values  for  some  of  the  inputs  of  the  model  described  in  this  paper. 
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CHAPTER  III 


PROBLEM  DESCRIPTION  AND  ASSUMPTIONS 

This  chapter  has  two  main  purposes.  The  first  purpose  is  to 
describe  the  basic  parts  of  the  mission  of  the  strategic  bomber  force. 
The  second  purpose  is  to  present  some  of  the  assumptions  concerning 
the  mission  that  will  be  made  in  this  model. 

An  overview  of  the  strategic  bomber  mission  is  shown  in  Figure 
3.1.  The  bombers  must  first  escape  from  their  bases  which  are  assumed 
to  be  under  SL8M  attack  at  the  outset  of  a  war  with  the  USSR.  At  some 
point  enroute  to  the  USSR,  the  bombers  will  receive  aerial  refueling 
from  tankers  which  have  escaped  from  their  own  bases  when  the  SLBM 
attack  was  launched.  The  bombers  will  then  encounter  a  forward  air 
defense  system  before  reaching  the  northern  coast  of  the  USSR  and  a 
perimeter  SAM  defense  immediately  after  reaching  the  coast. 

If  the  bomber  survives  to  this  point,  it  will  then  penetrate  to 
its  targets  through  area  defenses  which  may  consist  of  SAMs  and  fighter- 
interceptors.  At  each  target,  the  bomber  may  encounter  terminal  SAM 
defenses.  Finally,  after  releasing  all  its  weapons,  the  bomber  will 
attempt  to  leave  the  USSR  and  recover  at  a  friendly  base. 

The  mission  which  has  been  outlined  above  can  be  divided  into  eight 
separate  stages.  These  are  a)  base  escape,  b)  aerial  refueling,  c)  for¬ 
ward  air  defense,  d)  barrier  SAMs,  e)  random  area  SAMs,  f)  fighter- 
interceptors,  g)  weapons  delivered  and  terminal  SAMs,  and  h)  recovery. 

In  the  remainder  of  this  chapter,  each  of  those  stages  will  be 
described  in  greater  detail  and  the  major  assumptions  associated  with 
each  stage  will  be  discussed.  The  last  section  of  the  chapter  will 
discuss  cruise  missiles  and  the  assumptions  made  by  this  model  con¬ 
cerning  them. 
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FIGURE  3.1 

The  strategic  bomber  penetration  mission. 


BASE  ESCAPE 


It  is  a  commonly  accepted  assumption  that  a  preemptive  Soviet 
nuclear  strike  against  the  U.S.  would  include  a  SLBM  attack  on  U.S. 
strategic  bomber  bases.  The  bombers'  ability  to  escape  the  base  and 
avoid  the  blast  effects  of  the  Soviet  weapons  will  depend  largely  on  the 
amount  of  warning  they  receive  of  the  attack,  their  reaction  time,  and 
the  number,  size,  and  placement  of  the  warheads  targeted  on  the  base. 

At  the  time  of  the  attack,  some  part  of  the  bomber  force  stationed 
at  any  base  may  be  on  airborne  alert.  The  number  of  bombers  on  air¬ 
borne  alert  will  depend  greatly  on  the  degree  of  perceived  crisis 
immediately  prior  to  the  attack.  For  example,  increased  tension  between 
the  US  and  USSR  during  another  war  in  the  Middle  East  could  be  expected 
to  lead  to  an  increase  in  the  number  of  bombers  on  airborne  alert.  The 
situation  immediately  prior  to  attack  is  shown  in  Figure  3.2. 
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FIGURE  3.2 

Situation  at  base  prior  to  SLBM  attack. 
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When  the  first  warning  of  the  attack  is  received,  the  bombers  on 
the  ground  will  begin  to  take  off  and  dash  away  from  the  immediate 
vicinity  of  the  base.  The  number  of  bombers  that  can  get  airborne  will 
depend  on  the  amount  of  warning  they  receive,  their  reaction  time,  and 
the  rate  at  which  they  can  take  off  from  the  base.  Bombers  that  take 
off  between  the  time  the  first  warning  of  the  attack  is  received  and  the 
time  the  warheads  arrive  and  detonate  are  not  guaranteed  of  successfully 
escaping.  Should  a  bomber  fly  too  close  to  a  detonating  warhead,  it 
may  still  be  destroyed.  The  pattern  of  warhead  detonations  around  the 
base  will  depend  on  how  each  weapon  is  targeted  and  its  CEP  around  its 
aim  point.  The  situation  is  depicted  in  Figure  3.3. 
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FIGURE  3.3 

Situation  at  base  during  SLBM  attack. 
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This  model  makes  the  following  assumDtions  about  the  base  escape 
problem.  First,  it  is  assumed  that  all  bombers  on  airborne  alert  at 
the  time  of  the  attack  are  far  enough  away  to  be  immune  and  are  there¬ 
fore  assured  of  surviving  the  attack.  All  bombers  not  on  airborne 
alert  are  assumed  to  be  on  the  ground  at  the  base  when  warning  of  the 
attack  is  first  received. 

To  avoid  the  problems  of  warhead  targeting  and  the  location  uncer¬ 
tainty  of  each  warhead  due  to  its  CEP,  this  model  assumes  that  the 
detonating  warheads  will  be  located  randomly  in  the  airspace  around 
the  base.  However,  it  is  further  assumed  that  at  least  one  warhead 
detonates  close  enough  to  the  base  to  destroy  all  bombers  remaining  on 
the  ground  at  that  time. 

REFUELING 

Bombers  will  require  airborne  refueling  to  effectively  accomplish 
their  mission  as  planned.  Refueling  requirements  depend  on  the  bombers' 
payload,  fuel  capacity,  and  characteristics  of  its  flight  plan. 

Another  factor  is  the  amount  of  time  a  bomber  may  have  spent  on  air¬ 
borne  alert  prior  to  the  enemy  attack;  he  may  not  be  able  to  complete 
his  mission  on  the  fuel  load  he  has  remaining. 

Altitude  is  one  important  aspect  of  the  bomber's  flight  plan  that 
affects  his  fuel  requirements.  At  low  altitude,  a  bomber  stands  a 
better  change  of  escaping  detection  by  most  radars  and  may  have  a 
better  chance  of  surviving  an  attack  by  either  a  fighter  interceptor 
or  a  SAM  site.  This  is  particularly  true  in  the  case  of  many  older 
fighters  which  do  not  possess  the  so-called  "look-down,  shoot-down" 
capability.  Hence,  at  some  time  during  the  mission  a  bomber  would  be 
expected  to  go  to  low  altitude. 
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Unfortunately,  at  low  altitudes,  fuel  consumption  is  much  higher. 

In  order  to  penetrate  great  distances  at  low  altitude  the  bomber  must 
begin  the  penetration  part  of  its  mission  with  a  large  fue'!  load.  This 
necessitates  airborne  refueling  at  some  time  before  encountering  the 
enemy  defenses. 

With  the  current  bomber  force,  it  may  often  be  the  case  that 
multiple  refuelings  are  required.  This  will  become  especially  true 
when  cruise  missiles  are  added  to  the  bombers'  payloads.  Although 
requirements  will  vary  among  individual  bombers  with  different  mission 
plans,  a  typical  bomber  may  expect  to  refuel  twice  in  flight  before 
reaching  the  first  line  of  enemy  defense.  Each  missile  refueling  will 
potentially  degrade  a  bomber's  effectiveness. 

Successful  refueling  by  a  bomber  depends  on  two  main  issues.  The 
first  issue  is  the  availability  of  a  tanker  from  which  to  refuel. 

Tankers  that  are  based  at  the  same  bases  as  the  bombers  will  face  the 
same  problem  of  base  escape.  This  leaves  some  uncertainty  in  the  exact 
number  of  tankers  available  to  refuel  bombers.  The  number  of  bombers 
each  tanker  can  refuel  will  depend  on  the  mission  plan  and  the  fuel 
requirements  of  each  bomber  being  refueled.  Typically,  a  tanker  could 
expect  to  refuel  two  bombers  one  time  each. 

The  second  issue  is  that  of  a  bomber's  ability  to  locate  and 
couple  with  an  available  tanker  to  be  refueled.  Detailed  plans  and 
alternate  plans  exist  for  each  bomber  specifying  when  and  where  to  go 
to  get  refueled.  However,  navigational  error  enroute  to  the  refueling 
point  or  equipment  failure  on  the  part  of  either  bomber  or  tanker 
could  prevent  a  successful  refueling.  The  refueling  problem  Is  pictured 
in  Figure  3.4.  Because  of  the  complexity  of  the  problem,  several 
assumptions  are  made  for  the  purposes  of  this  model. 


25 


FIGURE  3.4 
Aerial  Refueling. 


The  first  assumption  deals  with  the  number  of  tankers  required  to 
serve  each  bomber.  Since  on  the  average,  bombers  require  two  refuelings 
apiece  and  tankers  are  capable  of  two  refuelings,  the  assumption  is 
made  in  this  model  that  one  tanker  is  required  to  refuel  each  bomber. 
This  does  not  imply  that  the  bomber  refuels  twice  from  the  same  tanker, 
which  is  unlikely,  only  that  the  equivalent  of  one  tanker  load  of  fuel 
is  necessary  to  refuel  the  bomber  twice. 

The  next  assumption  made  concerns  a  bomber's  probability  of  suc¬ 
cessfully  refueling,  given  that  there  is  an  available  tanker.  As  long 
as  a  tanker  is  airborne  for  the  bomber  to  refuel  from,  there  will  be 
a  certain  probability  that  it  can  successfully  refuel  once  the  tanker 
is  found.  The  first  probability  could  be  estimated  from  the  results 


of  numerous  simulated  missions  flown  each  year  and  may  be  quite  close 
to  1.  The  second  probability  could  be  obtained  from  equipment 
reliability  data,  perhaps  adjusted  to  allow  for  damage  that  may  occur 
to  tankers  during  base  escape.  Together  these  two  probabilities  define 
a  single  probability  of  a  bomber  successfully  completing  one  refueling. 
The  second  refueling  will  involve  a  different  tanker  and  take  place 
much  later  in  the  mission  and  the  same  process  will  be  repeated.  If  the 
same  probabilities  of  locating  and  successfully  refueling  apply,  then 
the  total  probability  of  refueling  twice  will  be  the  square  of  the 
probability  of  refueling  the  first  time. 

Since  the  bombers  and  tankers  are  following  detailed  plans,  they 
are  essentially  matched  bomber  to  tanker  and  are  not  searching  randomly 
for  each  other.  This  assures  that  when  a  bomber  locates  its  assigned 
tanker  there  will  be  fuel  available  for  it;  the  tanker  will  not  have 
already  off-loaded  its  fuel  onto  other  bombers. 

The  above  reasoning  leads  to  the  assumption  that  a  bomber's 
probability  of  successfully  refueling,  given  an  available  tanker,  can 
be  summed  up  in  one  number  which  takes  into  account  finding  the  tanker, 
successfully  refueling  from  the  ta n<er,  and  repeating  the  process  later 
in  the  mission.  This  probability  will  be  constant  for  all  bombers  and 
each  bomber's  refueling  event  will  be  independent  of  the  other  bombers'. 

The  final  assumption  deals  with  what  happens  to  bombers  that  fail 
to  get  refueled.  In  general,  a  bomber  that  does  not  get  refueled  will 
continue  its  mission.  It  will  still  be  capable  of  reaching  targets 
within  the  Soviet  territory.  It  may,  however,  be  forced  to  fly  much 
or  all  of  its  mission  at  high  altitude  where  it  is  more  easily  detected 
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and  attacked  by  Soviet  defenses  but  where  its  fuel  consumption  rate 
will  be  much  lower. 

For  the  purposes  of  this  model,  it  will  be  assumed  that  all  bombers 
that  fail  to  refuel  will  continue  their  missions  but  will  penetrate  at 
high  altitude  while  bombers  that  successfully  refuel  will  choose  to 
penetrate  at  low  altitude.  After  the  refueling  stage  of  the  mission, 
the  model  will  treat  the  refueled  and  unrefueled  bombers  as  separate 
and  essentially  independent  groups  for  the  remainder  of  the  mission. 

Also,  this  model  does  not  consider  the  case  where  a  bomber  receives 
only  one  refueling;  bombers  are  assumed  to  either  successfully  refuel 
twice  or  fail  to  refuel . 

FORWARD  AIR  DEFENSE 

The  first  line  of  defense  the  bombers  will  meet  is  the  forward 
air  defense  (FAD).  FAD  systems  attempt  to  attack  the  bombers  before 
they  reach  the  defended  continent  and  have  a  chance  to  release  their 
weapons . 

A  FAD  system  will  include  two  major  components.  One  is  some  sys¬ 
tem  of  warning  and  possibly  detection  of  approaching  bombers.  The 
other  component  will  be  a  force  of  fighter-interceptors  sent  out  to 
engage  the  bombers. 

The  system  will  operate  in  one  of  two  general  modes.  In  the  close 
control  mode,  a  sensor  or  sensor  network  is  used  to  detect  and  track 
the  incoming  bombers.  The  fighters  are  then  uniformly  assigned  to  the 
bombers  by  a  control  system  and  vectored  to  appropriate  intercept  points. 
Both  detection  and  control  systems  may  be  part  of  an  airborne  warning 
and  control  system  (AWACS).  An  AWACS  system  could  contain  radars 
capable  of  looking  both  out  and  down  and  would  be  able  to  detect  both 
high  and  low  flying  bombers. 
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The  second  mode  in  which  a  FAD  system  could  operate  is  the  raid 
control  mode.  In  this  case,  warning  is  received  of  the  approaching 
bombers  and  the  defense  will  commit  interceptors  to  an  airborne  search 
over  a  large  area.  Bombers  are  engaged  as  they  are  located  and  the 
interceptors  continue  the  search  while  fuel  and  ammunition  remain. 

This  model  assumes  a  FAD  system  composed  of  early  warning  and  con¬ 
trol  aircraft  (SUAWACS)  and  a  fighter  force  which  is  close  controlled 
as  in  Figure  3.5.  It  is  assumed  that  there  exist  sufficient  SUAWACS 
to  form  a  perimeter  of  overlapping  coverages  so  that  bombers  cannot 
leak  through  undetected.  It  is  assumed  that  the  SUAWACS  are  capable 
of  making  a  highly  accurate  count  of  the  incoming  bombers. 
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FIGURE  3.5 

Forward  Air  Defense. 
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Another  assumption  that  is  made  is  that  each  fighter  can  attack 
at  most  one  bomber.  In  support  of  this  assumption,  consider  the 
following  scenario. 

A  fighter  takes  off  from  its  base  and  flies  to  its  assigned  CAP 
position  where  it  may  orbit  for  a  while  until  a  bomber  is  detected. 

The  fighter  is  then  assigned  to  the  bomber  and  vectored  by  the  SUAWACS 
to  its  interception  point.  Having  arrived  there,  the  fighter  perhaps 
has  to  search  briefly  before  locating  the  bomber.  Once  the  bomber  is 
located,  the  FI  attempts  to  convert  on  and  kill  it.  After  the  attack, 
even  if  the  fighter  still  has  some  weapons  remaining,  fuel  must  be  a 
consideration.  If  the  approaching  bombers  are  separated  by  distances 
large  enough  to  prevent  the  fighter  from  immediately  locating  a  second 
target,  the  fighter  will  probably  have  to  return  to  base.  By  the  time 
he  has  returned,  landed  and  refueled,  the  FAD  battle  will  have  ended. 

(The  fighter  may,  however,  still  return  to  fight  as  part  of  the  area 
FI  defenses  of  the  heartland.)  Note  that  the  above  scenario  implicitly 
assumes  that  the  defensive  fighters  do  not  have  an  aerial  refueling 
capabil ity. 

BARRIER  SAMs 

After  leaving  the  FAD  zone,  the  bombers  will  approach  the  defended 
continent.  Shortly  after  reaching  the  coast,  they  may  encounter  a 
barrier  or  band  of  SAM  sites  (BSAM).  The  coverage  of  this  barrier  may 
not  extend  too  deep  into  the  defended  area,  but  the  barrier  would  pre¬ 
sumably  be  as  wide  as  possible.  This  will  present  the  penetrating 
bombers  with  a  wall  of  defenses  which  they  cannot  go  around  and  must 
therefore  break  through  before  penetrating  into  the  Soviet  heartland. 
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Ideally,  the  defense  would  like  to  place  the  SAM  sites  so  that 
their  coverages  overlapped  at  all  altitudes,  making  the  band  uniformly 
defended.  However,  the  defense  may  be  limited  in  the  total  number  of 
SAMs  he  has  available  to  deploy  and  the  location  of  the  sites  he  has 
will  be  effected  by  terrain  considerations.  In  practice  therefore, 
the  band  may  not  be  uniformly  defended;  in  some  places  site  coverages 
overlap  while  in  other  places  there  may  exist  gaps  in  the  defense. 

The  situation  is  shown  in  Figure  3.6. 


FIGURE  3.6 

Barrier  SAM  defense. 

If  the  offense  knew  with  reasonable  certainty  before  the  attack 
where  the  SAM  sites  were  located,  the  bombers  could  navigate  courses 
through  the  gaps  in  the  barrier.  Alternatively,  they  could  employ  ECM 
to  degrade  the  effective  radius  of  a  site  or  suppress  the  site  entirely 
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with  a  SRAM.  Any  of  these  measures  would  allow  the  offense  to  reduce 
or  avoid  altogether  losses  to  the  BSAMs. 

Unfortunately,  it  is  not  possible  for  the  offense  to  know  with  any 
certainty  the  location  of  SAMs  prior  to  the  attack.  Most  of  the  Soviet 
SAM  systems  are  mobile.  The  missiles  can  be  both  transported  and 
launched  from  tracked  or  even,  as  in  the  case  of  the  SA-3,  amphibious 
vehicles,  while  target  acquisition  and  fire  control  radars  can  be 
deployed  on  other  vehicles.  This  allows  an  entire  firing  battery  to  move 
at  will  around  the  country.  Furthermore,  this  capability  is  not  limited 
to  SAMs  of  small  effective  radius.  The  Mobile  SA-4,  for  example,  has 
a  range  of  36  nautical  miles  and  can  reach  targets  at  altitudes  of 
30,000  feet.  Other  SAM  features  include  guidance  systems  designed  to 
frustrate  ECM  measures.  The  aforementioned  SA-8  employs  an  electro- 
optical  tracker  in  conjunction  with  two  separate  guidance  radars  which 
allow  it  to  launch  two  missiles  at  a  target  on  two  different  guidance 
frequencies  (Ref  9) . 

The  model  makes  the  following  assumptions  about  the  BSAM  defense. 
First,  it  is  assumed  that  the  offense  avoids  or  suppresses  any  SAM 
sites  of  which  he  is  certain  of  the  location  prior  to  the  attack.  This 
implies  that  the  only  SAMs  that  will  effect  the  problem  are  mobile  SAMs 
which,  because  they  can  change  position,  cannot  be  briefed  before  the 
mission. 

The  model  also  assumes  a  homogenous  SAM  force.  Although  the 
different  types  of  SAMs  have  varying  effective  ranges  and  altitudes, 
for  this  model  they  are  all  assumed  to  be  described  by  the  same  average 
capabilities. 
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The  final  assumption  concerns  the  placement  of  the  mobile  SAMs  in 
the  barrier.  It  is  assumed  that  the  defense  cannot  achieve  the  ideal 
deployment  of  overlapping  SAM  coverages  for  the  entire  width  of  the 
barrier,  either  because  he  is  limited  in  resources  or  because  terrain 
prevents  him.  However,  if  the  SAMs  are  mobile,  the  offense  will  not 
know  where  the  gaps  are  at  any  given  time.  Still,  if  the  offense  knows 
some  points  along  the  barrier  are  less  likely  to  be  defended  than 
others,  then  he  will  presumably  choose  to  attack  through  those  points. 
Therefore,  in  order  for  the  barrier  to  be  most  effective,  the  defense 
will  wish  to  defend  all  parts  of  it  as  uniformly  as  possible.  Thus, 
this  model  assumes  that  at  any  given  time,  the  locations  of  the  SAMs 
are  distributed  uniformly  along  the  entire  width  of  the  barrier. 

RANDOM  AREA  SAMs 

After  the  bombers  break  through  the  barrier  SAMs,  they  begin  to 
penetrate  deeper  into  the  defended  territory  enroute  to  their  targets. 
Along  the  way  they  may  encounter  two  different  forms  of  area  defenses, 
random  area  SAMs  (RASAMs)  and  fighter-interceptors  (FI). 

Even  though  these  two  forms  of  defenses  cannot  be  separated  geo¬ 
graphically  like  the  FAD  and  BSAM,  this  model  makes  the  assumption 
that  they  can  be  treated  independently  of  each  other.  This  means  that 
an  individual  bomber's  probability  of  being  killed  by  the  RASAM  system 
is  independent  of  his  probability  of  being  shot  down  by  the  FIs.  The 
remainder  of  this  section  is  addressed  to  RASAMs,  while  FIs  are  dealt 
with  in  the  next  section. 

A  RASAM  defense  differs  in  two  major  ways  from  the  BSAMS.  First, 
the  depth  of  the  coverage  is  no  longer  a  neglible  factor.  It  is 
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assumed  that  the  RASAMs  may  be  encountered  anywhere  within  the  defended 
territory.  Instead  of  a  barrier  which  a  bomber  must  penetrate  once  and 
can  then  ignore  for  the  remainder  of  the  mission,  RASAMs  present  a  con¬ 
tinuous  threat  from  the  time  penetration  begins  until  the  bomber  leaves 
the  defended  area. 

The  second  difference  is  related  to  the  first.  Since  RASAMs  are 
distributed  throughout  the  defended  area,  they  are  necessarily  inter¬ 
spersed  with  the  targets  in  the  area.  A  bomber  need  not  penetrate 
through  the  entire  depth  of  the  RASAM  net  to  reach  his  targets.  Some 
targets  may  be  reached  only  a  short  way  into  the  net;  more  targets 
will  be  exposed  as  deeper  penetration  is  made. 

Although  the  above  difference  between  RASAMs  and  BSAMs  are  impor¬ 
tant  in  the  model's  treatment  of  the  former,  there  are  also  important 
similarities  between  the  two.  Since  bombers  could  be  expected  to  fly 
around  or  suppress  any  known  SAM  sites,  the  assumption  is  made  again 
that,  because  of  their  location  uncertainty,  the  only  SAMs  of  interest 
to  the  model  are  mobile  SAMs.  Also,  it  is  assumed  that  the  RASAMs  are 
dispersed  randomly  and  uniformly  through  the  defended  area.  In  prac¬ 
tice,  the  placement  of  the  RASAMs  might  be  affected  by  the  defense's 
knowledge  of  the  locations  of  the  targets  within  the  defended  area. 
However,  if  the  offense  believed  that  the  actual  distribution  of  RASAMs 
was  significantly  different  from  uniform  they  might  find  it  advantageous 
to  plot  courses  around  suspected  high  concentrations  of  RASAMs. 

The  RASAM  defense  is  shown  in  Figure  3.7. 
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FIGURE  3.7 
Random  Area  SAMs. 

FIGHTER-INTERCEPTORS 

The  second  form  of  area  defense  that  bombers  may  encounter  is 
fighter-interceptors  (FI).  As  a  bomber  penetrates  into  the  defended 
area,  it  will  come  within  range  of  fighters  stationed  at  various  air¬ 
bases  throughout  the  country.  As  the  bomber  slips  within  range  of  a 
base,  any  fighters  stationed  at  that  base  which  are  not  already  engaged 
with  other  bombers  may  be  assigned  to  attack  it.  As  in  the  case  of 
the  FAD,  these  fighters  may  be  either  raid  controlled  or  close  con¬ 
trolled. 

In  the  course  of  the  penetration,  there  may  be  many  such  FI 
attacks  on  the  bomber.  The  number  will  depend  on  the  bomber's  total 
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depth  of  penetration,  the  number  of  airbases  within  FI  range  of  the 
bomber's  course,  and  the  FI  commitment  policy  of  the  defense.  It 
may  also  depend  on  the  total  number  of  offensive  penetrators  since  this 
will  be  one  factor  determining  the  number  of  FIs  available  to  attack  a 
particular  bomber. 

In  the  case  of  close  controlled  intercepts,  a  typical  bomber- 
interceptor  engagement  will  follow  roughly  the  same  scenario  presented 
in  the  FAD  section.  The  fighter  would  be  vectored  by  the  control  system 
to  an  appropriate  intercept  point.  There  it  would  search  for  the  bomDer 
and  after  detecting  it  would  attempt  to  convert  on  it.  The  bomber  can 
to  some  extent  degrade  the  FI's  capability  to  convert  on  it  by  jamming 
the  FI's  radar,  denying  the  FI  range  information.  If  the  FI  successfully 
converts  on  the  bomber,  it  will  then  release  one  or  more  air-to-air 
missiles.  After  the  engagement,  the  fighters'  remaining  fuel  and 
armament  loads  and  the  location  of  other  detected  bombers  will  deter¬ 
mine  whether  it  can  attack  a  second  bomber  or  must  return  to  base  to  be 
recycled. 

To  further  define  the  specific  problem  being  addressed,  several 
assumptions  are  made  by  this  model. 

The  first  assumption  concerns  the  FI  force.  There  are  over  a 
dozen  different  fighters  and  fighter-interceptors  in  the  Soviet  inven¬ 
tory.  They  range  in  capabilities  from  the  subsonic  MiG-17,  dating  back 
to  the  early  1 950 * s  and  having  a  range  of  under  300  nm,  to  the  all- 
weather  SU-15,  with  a  nearly  400  nm  radius  and  capable  of  speeds  in 
excess  of  Mach  2,  and  the  MiG-25  with  a  600  nm  radius  and  capable  of 
speeds  of  Mach  2.8  (Ref  10;  14:178-222).  As  was  done  for  the  SAMs, 
this  model  assumes  a  homogeneous  FI  force  having  capabilities  that  are 
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average  or  typical  of  the  various  fighter-interceptors  that  bombers 
could  actually  expect  to  encounter.  These  capabilities  need  not 
necessarily  be  the  same  as  those  of  the  FAD  fighter  force. 

The  next  assumptions  deal  with  the  defensive  command  and  control 
system.  It  is  assumed  in  this  model  that  all  the  fighter-interceptors 
are  close  controlled,  as  was  the  case  in  the  FAD  defense,  f in like  the 
FAD  system,  it  is  assumed  that  in  the  area  defense  the  interceptors 
are  controlled  from  a  ground  based  detection  and  control  system  rather 
than  SUAWACS.  This  GCI  net  will  detect  approaching  bombers  and  then 
vector  available  FIs  to  appropriate  intercept  points.  The  scenario 
is  shown  in  Figure  3.8.  The  model  assumes  that  the  chance  of  an 
encounter  between  a  previously  undetected  bomber  and  an  airborne  FI  on 
CAP  or  enroute  to  an  assigned  intercept  point  in  neglible.  The  model 
also  assumes  that  the  GCI  net  and  the  RASAMs  operat°  independently. 


FIGURE  3.8 

Area  fighter-interceptor  defense. 
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To  effectively  close  control  the  interceptors,  the  defense  must 
be  capable  of  tracking  bombers  over  the  entire  defended  area.  It  is 
assumed  therefore  that  the  defensive  GCI  net  is  spread  over  the  defended 
area.  As  the  bomber  begins  penetration,  he  may  be  able  to  avoid  detec¬ 
tion  by  the  net  for  some  time  by  flying  at  low  altitudes.  However, 
once  the  bomber  is  detected  by  the  radar  net,  it  is  assumed  that  the 
defense  will  be  able  to  track  the  bomber's  course  and  vector  intercep¬ 
tors  to  it  for  the  remainder  of  the  mission. 

Several  other  assumptions  will  be  introduced  as  the  FI  part  of  the 
model  is  developed.  They  will  be  described  in  the  section  on  FI  area 
defenses  in  Chapter  V. 

TERMINAL  SAM,  WEAPONS  DELIVERED,  AND  VALUE  EXTRACTED 

The  assumptions  and  analysis  involved  in  any  model's  representa¬ 
tion  of  terminal  SAM  defenses,  weapons  delivered,  and  value  extracted 
or  destroyed  by  the  offense  are.  very  interrelated.  These  three  areas 
are  treated  together  in  this  section. 

It  seems  reasonable  to  assume  that  the  offense  and  defense  will 
be  in  rough  agreement  on  what  makes  up  the  target  set.  Differences 
may  arise,  however,  in  the  relative  values  each  places  on  the  targets 
within  the  target  set.  The  defense  will  try  to  defend  the  targets  in 
a  way  it  finds  optimal  according  to  the  relative  values  it  puts  on  the 
targets.  The  offense  will  try  to  optimize  its  attack  using  its  own 
value  system,  which  may  differ  from  that  of  the  defense. 

It  has  been  shown  that  the  optimum  strategy  for  the  defense  is  to 
defend  each  target  relative  to  the  other  targets  in  proportion  to  that 
target's  value  relative  to  the  other  targets  (Ref  13:9).  In  other 
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words,  if  on  some  scale  of  value,  the  defense  considers  target  S  to 
be  twice  as  valuable  as  target  A,  then  they  should  defend  target  3 
with  twice  as  many  SAMs  as  are  used  at  target  A.  This  system  assumes 
that  it  is  possible  to  make  such  relative  value  judgments  and  also 
assumes  that  two  SAMs  can  defend  a  target  twice  as  effectively  as  one. 

If  these  assumptions  are  accepted,  then  for  a  given  number  of  defensive 
SAMs  and  a  given  target  set  of  known  relative  values  the  optimal  de¬ 
fense  can  be  found. 

If  the  offense  places  the  same  relative  values  on  the  targets  as 
the  defense  does,  then  the  above  system  will  present  the  offense  with 
a  constant  price  that  must  be  paid  to  extract  each  unit  of  value. 

Using  the  example  above,  the  offense  must  expend  twice  as  many  resources 
to  be  sure  of  destroying  target  B  as  it  would  to  destroy  target  A. 

This  would  leave  the  offense  indifferent  to  which  target  it  attacked 
since  the  total  value  it  could  destroy  would  be  proportional  to  its 
total  resources  no  matter  which  system  it  used  to  select  targets  to 
attack. 

The  situation  is  different  if  the  offense  places  different  rela¬ 
tive  values  on  the  targes  than  the  defense.  In  this  case,  the  offense 
can  optimize  its  attack  by  attacking  those  targets  which  it  perceives 
as  being  defended  to  a  level  less  than  their  value  compared  to  other 
targets.  Still  using  the  above  example,  if  the  offense  perceives 
targets  A  and  B  as  equally  valuable  then  it  will  choose  to  attack  tar¬ 
get  A  since  it  is  defended  only  half  as  much  as  target  B. 

So  far,  the  entire  analysis  presented  above  has  been  based  solely 
on  considerations  involving  relative  values  of  various  targets. 
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Another  important  factor  that  has  thus  far  been  ignored  is  a  target's 
distance  inside  the  defended  area.  Consider  again  two  targets  A  and 
B  of  equal  value,  but  now  let  target  B  lie  at  a  greater  depth  inside 
the  defended  area  than  target  A.  The  offense  will  be  more  likely  to 
attack  target  A  since  a  bomber  will  stand  a  greater  change  of  pene¬ 
trating  to  A  through  the  area  defenses  than  it  would  of  penetrating 
the  longer  distance  to  target  B.  Knowing  this,  the  defense  may  in¬ 
crease  the  level  of  terminal  defenses  at  target  A  to  try  to  make  the 
offense  indifferent  to  the  two  targets. 

To  make  even  a  reasonably  detailed  probabilistic  model  of  the 
situation  described  above  would  require  several  things.  Some  knowledge 
of  the  target  set  and  the  relative  values  of  the  targets  would,  of 
course,  be  necessary.  Then  it  is  possible  that  one  would  require  some 
form  of  weapon  allocation  scheme  to  determine  the  number  of  weapons  to 
assign  to  each  target.  Finally,  a  distribution  would  have  to  be  found 
giving  the  probability  that  a  target  is  destroyed  as  a  function  of  the 
level  of  defense  and  the  number  of  offensive  penetrators  assigned  to 
attack  it. 

All  that  is  beyond  the  scope  of  this  model.  To  simplify  the  pro¬ 
blem  into  something  that  can  be  more  easily  modelled,  the  following  two 
assumptions  are  made.  First,  it  is  assumed  that  all  targets  are  of 
equal  value  and  are  perceived  to  be  so  by  both  offense  and  defense. 
Second,  it  is  assumed  that  neither  side  considers  the  location  of  the 
target  when  planning  its  attack  or  defense.  These  assumptions  clearly 
differ  significantly  from  reality  but  will  allow  some  representation 
of  the  terminal  defense  and  value  extracted  aspects  of  the  bomber  mis¬ 
sion  without  making  the  model  unreasonably  complicated. 
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The  above  assumptions  lead  directly  to  three  others.  First,  if 
all  the  targets  are  of  equal  value,  then  it  can  reasonably  be  assumed 
that  they  will  all  be  uniformly  defended.  In  other  words,  each  target 
will  be  defended  by  the  same  number  of  terminal  SAMs.  Second,  it  can 
also  be  assumed  that  the  offense  will  allocate  the  same  number  of 
weapons  to  each  target.  Finally,  if  each  target  is  of  equal  value, 
then  it  follows  that  the  total  value  destroyed  by  the  offense  will  be 
directly  proportional  to  the  number  of  targets  destroyed. 

To  determine  the  number  of  targets  destroyed,  two  more  assumptions 
will  be  made.  It  will  be  assumed  that  only  one  weapon  is  released  at 
each  target  and  that  weapon  will  totally  destroy  the  target  if  it 
detonates  where  it  is  aimed.  In  other  words,  the  weapon  either  suc¬ 
ceeds,  in  which  case  the  target  is  destroyed,  or  fails,  in  which  case 
the  target  is  unharmed.  There  is  no  partial  damage  to  the  target. 

RECOVERY 

After  a  bomber  has  released  its  last  weapon,  the  final  part  of 
its  mission  will  be  to  get  out  of  the  defended  country  and  find  a  safe 
place  to  land.  If  the  bomber  succeeds  in  doing  that,  then  a  variety 
of  factors  at  the  time  will  determine  if  it  can  ever  be  refueled, 
rearmed,  and  sent  on  a  second  mission. 

A  complete  analysis  of  all  parts  of  the  recovery  problem  is 
beyond  the  scope  of  this  model.  This  model  will  only  deal  with  the 
bomber's  ability  to  get  out  of  the  defended  area.  All  subsequent 
aspects  of  the  problem  such  as  locating  a  safe  base  at  which  to  land 
will  not  be  considered. 
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This  model  assumes  that  after  a  bomber  has  released  its  last 
weapon  it  will  continue  to  penetrate  south  along  the  corridor  until 
it  exits  at  the  bottom,  or  southern  end.  At  this  point  tne  bomber 
will  be  considered  to  be  outside  the  defended  area  and  its  mission, 
in  terms  of  this  model,  ended. 

It  seems  reasonable  to  assume  that  the  defense  will  not  know  when 
the  bomber  has  exhausted  its  weapons  and  is  no  longer  a  threat.  There¬ 
fore,  attacks  on  the  bomber,  both  by  fighter-interceptors  and  SAMs,  will 
continue  until  the  bomber  leaves  the  defended  area.  In  other  words, 
the  bomber  continues  to  face  similar  defenses  while  leaving  the  defended 
area  as  he  did  when  penetrating  to  reach  his  targets.  The  only  difference 
is  that  there  is  assumed  to  be  no  BSAM  defense  at  the  far  south  south¬ 
ern  end  of  the  corridor.  The  defense  is  not  concerned  about  an  attack 
from  the  south  and  is  not  interested  in  destroying  bombers  that  are 
leaving  the  defended  area  and  are  clearly  no  longer  a  threat. 

CRUISE  MISSILES 

Up  until  now,  the  subject  of  cruise  missiles  and  their  effect  on 
the  mission  of  the  air  breathing  element  has  not  been  addressed.  This 
will  be  the  subject  of  this  section. 

To  begin  with,  several  assumptions  will  be  made.  First,  it  will 
be  assumed  that  each  cruise  missile  is  released  at  its  maximum  range 
from  its  target.  Compared  to  other  weapons,  cruise  missiles  are 
fairly  heavy.  The  farther  they  are  carried  by  the  bomber,  the  more 
fuel  is  consumed  by  the  bomber  which  in  turn  limits  the  amount  of  other 
ordnance  the  bomber  can  carry.  Furthermore,  each  individual  penetrator 
may  have  a  better  chance  or  survival  the  more  penetrators  there  are. 
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If  a  bomber  gets  shot  down,  the  total  loss  to  the  offense  will  be  less 
if  he  is  not  still  carrying  his  cruise  missiles.  These  are  two  reasons 
for  the  bombers  to  release  their  cruise  missiles  as  early  as  possible 
so  the  above  assumption  seems  reasonable. 

Consistent  with  the  discussion  of  value  destroyed  in  the  previous 
section,  the  next  assumption  is  that  only  one  cruise  missile  is  fired 
at  each  target  and  if  the  cruise  missile  detonates  at  its  target,  the 
target  will  be  completely  destroyed.  The  targets  of  the  cruise  missiles 
will  not  be  the  same  as  those  of  the  bombers  themselves. 

The  last  general  assumption  is  that  all  the  defenses  that  threaten 
the  bombers  may  also  threaten  the  cruise  missiles.  The  various  defenses 
may  not  be  equally  effective  against  cruise  missiles  as  against  bombers 
but  they  have  the  potential  to  be  a  threat.  This  means  that  cruise 
missiles  will  face  TSAMs,  fighter-interceptors ,  RASAMs,  and  if  they  are 
released  early  enough,  BSAMs  and  eve'n  FAD  interceptors. 

The  forward  air  defense  part  of  the  model  would  become  significantly 
more  complicated  if  cruise  missiles  were  released  during  the  FAD  engage¬ 
ment  and  so,  for  modelling  convenience,  it  will  be  assumed  that  this 
does  not  occur.  That  is,  the  cruise  missiles  may  be  released  before 
or  after  the  FAD,  but  not  during  the  actual  engagement.  Also,  in  this 
model,  the  depth  dimension  of  the  BSAM  defense  has  been  assumed  to  be 
insignificant,  so  it  is  reasonable  to  assume  that  cruise  missiles  will 
not  be  released  within  the  BSAM  band.  Thus,  within  the  structure  of 
this  model,  there  are  three  different  possible  assumptions  that  can 
be  made  about  cruise  missile  release  points. 

The  first  assumption  would  be  that  all  the  cruise  missiles  are  of 
fairly  short  range  and  are  going  to  targets  reasonably  deep  within  the 
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defended  area.  In  this  case,  the  bombers  would  not  release  the  cruise 
missiles  until  after  they  had  passed  through  the  BSAM  defense  and  begun 
actual  penetration  of  the  defended  area. 

A  second  case  would  involve  cruise  missiles  of  longer  range.  In 
this  case  it  might  be  reasonable  to  assume  that  the  bombers  release 
their  cruise  missiles  after  penetrating  through  the  FAD  zone  but  be¬ 
fore  entering  the  BSAM  band  and  beginning  actual  penetration.  In 
practice,  the  FAD  zone  might  extend  right  up  to  the  coastline  of  the 
continent  and  the  BSAM  band  which  would  leave  no  place  between  the  two 
layers  of  defense  for  the  bombers  to  release  their  cruise  missiles. 
Still,  the  above  assumption  might  be  considered  equivalent  to  having 
the  bombers  release  the  cruise  missiles  so  late  in  the  FAD  engagement 
that  there  is  no  significant  effect  on  the  outcome  of  the  engagement. 

The  third  possible  case  would  involve  cruise  missiles  of  extremely 
long  range  in  which  case  it  could  be  assumed  that  the  cruise  missiles 
are  all  released  before  the  bombers  enter  the  FAD  zone.  At  this  time, 
that  may  not  be  a-  realistic  assumption  except  possibly  in  cases  where 
the  cruise  missiles  are  all  destined  for  targets  in  the  far  northern 
part  of  the  defended  area.  Nonetheless,  it  is  included  here  for  com¬ 
pleteness. 

The  cruise  missiles  currently  being  produced  by  Boeing  have  a 
range  of  over  1500  miles.  Thus,  in  an  actual  mission,  the  most  likely 
case  would  be  some  combination  of  the  above  involving  at  least  the  first 
two  and  possibly  all  three  cases.  In  other  words,  for  cruise  missiles 
of  a  given  range  and  a  target  set  distributed  over  the  entire  defended 
area,  some  of  the  cruise  missiles  may  be  released  before  the  bombers 
encounter  the  BSAM  defense  or  even  before  the  FAD  zone  while  others. 
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those  going  to  southern  targets,  will  have  to  be  carried  through  the 


BSAM  band  and  partway  through  the  defended  area, 
three  cases  are  treated  separately;  however,  with 
tional  assumptions,  it  should  not  be  difficult  to 
tions  of  cases. 


In  thi s  model ,  the 
a  few  minor  add i - 
represent  combina- 
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CHAPTER  IV 


PRIOR  TO  PENETRATION 

Purely  for  convenience,  the  various  parts  of  this  model  have  been 
divided  up  into  three  groups  which  will  be  described  in  this  chapter 
and  the  two  chapters  following.  This  chapter  covers  the  parts  of  the 
model  that  deal  with  the  stages  of  the  bomber's  mission  that  take  place 
before  actual  penetration  of  the  defended  area  begins.  In  this  model, 
actual  penetration  is  somewhat  arbitrarily  defined  to  begin  after  the 
bomber  has  gone  through  the  BSAM  defense.  Therefore,  the  parts  of  the 
model  described  in  this  chapter  are  base  escape  (BE),  airborne  refueling 
(AR),  forward  air  defense  (FAD)  and  the  barrier  SAMs  (BSAMs). 

BASE  ESCAPE 

Consider  a  base  at  which  there  are  stationed  B  bombers.  In  a  pre¬ 
emptive  attack  by  SLBM  foreces,  a  total  of  M  warheads  are  targeted  at 
the  base.  Assume  temporarily  that  at  the  time  of  the  attack,  all  the 
bombers  are  on  the  ground.  The  problem  is  to  find  a  distribution  for 
the  number  of  bombers  that  survive  the  attack. 

When  warning  is  first  received  of  the  attack,  the  bombers  will 
begin  to  take  off  and  fly  away  from  the  base.  It  is  assumed  that  when 
the  warheads  arrive  and  detonate,  all  bombers  that  are  not  airborne 
are  destroyed.  Bombers  that  are  airborne  will  have  some  probability 
of  survival  but  will  not  be  assured  of  surviving.  The  problem  can  now 
be  separated  into  two  parts;  first  finding  the  number  of  bombers  air¬ 
borne  at  the  time  of  warhead  detonation  and  second,  finding  the  prob¬ 
ability  that  an  airborne  bomber  will  survive  the  attack. 
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To  find  A,  the  number  of  bombers  airborne  when  the  missiles 

detonate,  it  is  necessary  to  know  the  amount  of  warning  the  base  has 

received  of  the  attack  and  the  rate  at  which  bombers  can  take  off  from 
the  base.  The  takeoff  rate  is  assumed  to  be  a  constant  depending  on  the 
number  of  runways  at  the  base  and  the  interval  between  takeoffs  of 
consecutive  bombers.  The  warning  time,  T  ,  is  a  random  variable  that 

could  depend  on  the  total  missile  flight  time,  the  time  at  which  the 

missile  is  detected  after  breaking  water,  and  the  command  and  control 
time  required  to  evaluate  the  data,  make  a  decision  and  notify  the  base 
of  attack. 

When  warning  of  the  attack  is  first  received,  there  may  be  a  delay 
before  bombers  can  start  taking  off  due  to  the  reaction  time  of  the  air¬ 
crews.  Call  this  delay  T^.  The  total  amount  of  time  available  for  the 
bombers  to  take  off  in  is  then  T^,  where 

TA  *  Tw  -  TR  ^ 

The  maximum  theoretical  time  in  which  bombers  could  take  off  is 

the  total  of  missile  flight  time,  T^,  while  the  minimum  is  zero.  The 

former  would  require  immediate  detection  of  the  missiles  as  they  break 

the  water  and  crews  already  aboard  the  bombers;  the  latter  would  occur 

if  the  warning  time,  T  ,  was  less  than  the  reaction  time,  T„.  The  most 

likely  time  would  be  somewhere  between  the  two  extremes.  Assume  that 

the  probability  density  function  of  the  time,  pT(  t),  can  be  fully 

'A 

specified.  It  might  have  a  shape  similar  to  the  one  in  Figure  4.1. 
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FIGURE  4.1 

Probability  density  function  for 

With  the  assumption  of  a  constant  take  off  rate,  the  number  of 
bombers  airborne  when  the  warheads  detonate  is  given  by 

A  =  min  (<rT^> ,  B)  4.2 

where  r  is  the  rate  at  which  bombers  take  off.  The  pointed  brackets 
are  read  "greatest  integer  in"  and  account  for  the  fact  that  the  number 
of  bombers  airborne  must  be  an  integer. 

The  probability  that  an  airborne  bomber  is  killed  can  be  found 
by  taking  the  ratio  of  the  volume  of  airspace  around  the  detonating 
missiles  within  which  the  bomber  would  be  destroyed  to  the  total 
volume  of  airspace  in  which  the  bomber  could  be  located. 

It  is  assumed  that  a  detonating  warhead  has  some  lethal  radius, 
rL,  within  which  any  bomber  is  destroyed  and  outside  of  which  a  bomber 
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will  survive  undamaged,  When  the  warhead  detonates,  there  will  be  a 
lethal  volume,  V.  ,  around  it  inside  of  which  any  bomber  is  killed 
with  probability  one.  If  the  warhead  is  assumed  to  be  detonated  above 
the  ground,  then  this  volume  can  be  approximated  by  a  sphere. 

V.  =  (4/3)  -rL2  4.3 

Bombers  outside  this  volume  survive  with  probability  one.  For  the 
case  of  an  attack  involving  M  warheads,  an  upper  bound  for  the  total 
lethal  volume,  V,  ,  of  airspace  is 

VL  =  (4/3)MirrL3  4.4 

This  is  the  case  where  the  lethal  volumes  from  individual  warheads  do 
not  overlap.  It  is  the  most  offense  conservative  case  and  to  avoid 
the  problem  of  trying  to  calculate  the  overlaps,  this  is  the  case 
employed  in  the  model . 

The  total  volume  within  which  bombers  may  be  located  can  be 
described  by  a  cylinder  centered  on  the  base  with  a  radius  R  and 
height  H  given  by 


R  =  Va 

H  =  vcTA 


4.5 


where  v.  and  v„  are  the  horizontal  velocity  and  rate  of  climb  of  a 
he 

bomber.  The  total  volume  of  the  cylinder  is 


VT  =  *JTH 

2  ,  3 

=  *vh  vc  TA 


4.6 
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If  the  assumption  is  made  that  airborne  bombers  are  equally  likely 
to  be  anywhere  inside  Vy  then  the  probability  that  an  airborne  bomber 
is  killed,  given  is  the  ratio  of  the  two  volumes. 


max  (1, 


4.7 


The  function  has  been  written  to  indicate  that  the  probability  of 
being  killed  cannot  be  greater  than  one. 

From  Equation  4.7,  the  probability  that  an  airborne  bomber  sur¬ 
vives  the  attack  given  is 


Pc  I  =  1  -  P 1/  |  T 
S '  '  A  •  A 


=  1  -  max  (1 , 


4Mr, 


,  2  3 

3vh  Va 


4.8 


It  is  assumed  that  for  a  given  TA,  each  airborne  bomber's  sur¬ 
vival  event  is  dependent  of  the  survival  events  of  the  other  bombers. 
The  distribution  of  the  number  of  bombers  that  survive  for  a  given 
T^,  is  therefore  binomial 


pi>e<K!V  * 


<>sitk  '1-pslt>A'k 


4Mr, 


•  (m1n(<Jt>>  B>)  [l-min(l, - ^ - 5)1 


4.9 


3rn  rcTA 


4Mr 


[mind,  l  ))H<n(<rt».  B)  -  k 

*5.#  T  ^ 


3vn  vcTA 


The  distribution  of  the  number  of  bombers  that  survive  can  now 
be  found  by  conditioning  on  T^. 
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4.10 


pie(k> 


s!> !i  -  1MrL 


,  2  .3 

3v,  v,.t 

h  C 


-)i 


[ mi n ( 1 , 


4Mr 


,  2  .3 

3vh  V 


L _ ^min(<rt>,  B) 


In  order  for  k  bombers  to  survive,  at  least  k  must  be  airborne  at  the 
time  the  missiles  detonate.  From  Equation  4.Z,  it  is  found  that  the 
smallest  value  of  which  would  allow  k  bombers  to  be  airborne  is 
k/r.  This  explains  the  lower  limit  on  the  integration. 

It  would  be  possible  to  numerically  integrate  equation  4.10  on 
a  computer  to  obtain  the  exact  distribution  of  the  number  of  bombers 
that  survive  the  attack.  However,  it  may  be  more  convenient  to  use 
the  following  approximation. 

First  break  the  integral  up  as  follows: 


(k+2)/r(.min{<rt>,B)lp  k,,  „  smin(<rt>  ,b)-k  , 

+  ;(k+l)/rl  k  ;PSl  t  U’PS|t;  pTi  t>at 


rr"(<r'B)>psit'1-ps:t>”1n(<rt>,B,'kpTA<t)-“ 


+  'l«r  ("1n(<T'B))PskltU-ps|t! 


min(<rt>  ,B)-k 


pT  (t)dt  4.11 
'A 


Consider  the  value  of  <rt>  in  each  of  the  integrals  in  Equation 


3.11. 


<rt>  =  k+i  Vt:(k+i)/r  <_  t  <(k+i+l)/r 


4.12 
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This  implies  that  max(<rt>).  S)  will  be  constant  in  each  interval  and 
can  be  taken  outside  the  integrals. 


If  one  assumes  that  P<- .  +  is  nearly  constant  within  each  interval, 
then  it  can  be  approximated  by  P<.,*  where  t-  is  the  average  value  of 

•  '•i  1 

t  h 

t  within  the  iLn  interval.  If  this  approximation  is  made,  then 
Ps;^(l-Ps  t.)m1n('rt>’B)-k-  a^so  taken  outside  each  integral, 
giving 


Pie!k>  ■  (k>PS't,k(1-Ps:t>k  .:^I1/rPr  (t)dt 


(k+l\p  ^fi_p  \k+l-k  .( k+2 )  /r  ftldt 
(  k  )Ps;t?(1  Ps;t„]  ' (k+l)/rPT,(t)dt 


4  (B"k1)PSlt  k  (1-PS  t  )B"1‘k  flf l/rPT  (t)dt 
k  5 • tB-k  S,tB-k  B-l/r  Ta 


+  k  (1-Po. t  )B_k  ;Tm  pT  (t)dt 

K  a-k+1  B-k+1  B/r  *A 


Now  consider  the  integrals  remaining  in  equation  4.13.  Each  integral 
can  be  interpreted  as  the  probability  that  T  lies  between  the  upper 
and  lower  bounds  of  integration.  However,  by  Equation  4.2,  this 
probability  can  further  be  interpreted  as  the  probability  that  the 
number  of  bombers  airborne  when  the  warheads  detonate  is  r  times  the 
lower  limit  of  integration.  Mattematically, 

/(kti|/r/r  Pt  "  pf(k+i)/r  <  TA  <  (k+i+1  )/r) 

A 

=  P( k+i  <_  rTA  <  k+i+1) 


P(k+1  <_  <rTA>  <  k+i+1) 


=  P  [max  ( k+i ,  B )  <_  max  (<rT^>,3) 

<  max(k+i+I ,B) ] 

=  P( A  =  max  (8 ,k+i ) )  4.14 

Substituting  equation  4.14  into  4.13  and  simplifying  (substitute 
j+k+i),  one  finds 

pie<k>  (k>PS  ,'-15 

j=k  (  i 

In  deriving  equation  4.15,  the  assumpcion  was  initially  made  that 
all  the  bombers  at  the  base  are  on  the  ground  when  the  attack  occurs. 

In  reality,  it  is  quite  probable  that  some  of  the  bombers  would  be  on 
airborne  alert  at  the  time  of  the  attack.  These  bombers  can  be  assumed 
to  survive  the  attack  on  the  base  with  probability  one.  If  8^  is  the 
number  of  bombers  airborne  and  8^  is  the  number  on  the  ground  at  the 
time  of  attack  where  and  sum  to  B;  then  with  this  assumption, 
the  distribution  of  the  number  of  bombers  that  escape  can  be  written 

Pbe(k)  .0  k  <  ba 

-  pbe(k-eA’ 

■  0  k  >_  8  4.16 

The  above  expressions  apply  to  a  single  base.  If  one  considers 
many  bases  having  different  characteristics  such  as  different  numbers 
of  bombers  assigned  to  them  or  different  distances  to  the  point  of 
SLBM  launch,  then  to  get  the  total  number  of  bombers  that  escape,  one 
must  sum  over  all  the  bases.  For  example,  for  n  bases: 
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B 


P3E( k) 


n-1 '  1)  (2) 

2  Pbe(k1)Pbe(,c2) 
n-1 


(n-1) 

> 

be 


(k 


n-1 ) 


(n) 

P5e(k-k1-k2- 


*  kn-l) 


4.17 


Alternatively,  if  one  assumes  the  bases  are  uniform  in  every 
respect,  including  having  the  same  number  of  warheads  assigned  to  each 
one,  then  the  entire  force  can  be  treated  as  one  large  base  having  nB^ 
bombers  in  the  air,  nB^  bombers  on  the  ground,  and  a  bomber  takeoff 
rate  of  nr.  T^,  the  time  of  flight  of  the  missiles,  will  remain  the 
same. 

It  may  be  that  in  some  cases  the  assumption  of  identical  bases 
would  differ  more  from  reality  than  the  modeller  wished  to  allow,  in 
which  case  Equation  4.17  would  have  to  be  used.  However,  because  of 
the  savings  in  computational  difficulty  to  be  gained,  this  model 
assumes  that  all  bases  are  identical  and  treats  the  bomber  force  as 
one  large  base. 

REFUELING 

The  next  problem  is  that  of  aerial  refueling.  Consider  a  force 
of  BgE  bombers  that  have  successfully  escaped  the  SLBM  attack  on  their 
bases  and  now  desire  inflight  refueling.  These  bombers  are  to  be 
refueled  by  a  force  of  T  tankers.  Some  of  the  tankers,  T$,  either 
because  they  are  airborne  at  the  time  of  attack  or  because  they  are 
stationed  at  a  "safe"  base,  are  assured  of  surviving  the  initial  SLBM 
attack.  The  remaining  tankers,  Tg,  are  on  the  ground  at  bases  that 
are  attacked  and  so  face  a  base  escape  problem  identical  to  that  of 
the  bombers.  Given  these  two  forces  of  aircraft,  the  problem  is  to 
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find  a  distribution  for  the  number  of  bombers  that  are  successfully 
refuel ed . 

Consider  first  the  case  where  BgF  is  less  than  or  equal  to  T<- . 

In  this  case  there  will  likely  to  be  enough  tankers  to  refuel  all 
the  bombers,  and  the  number  of  bombers  that  can  get  refueled  is  limited 
only  by  the  total  number  of  bombers,  Bg^,  requiring  refueling.  Each 
bomber  is  assumed  to  have  some  probability  of  successfully  refueling 
given  an  available  tanker  with  sufficient  fuel.  Since  it  is  assumed 
that  each  bomber's  refueling  event  is  independent  of  all  the  other 
bombers'  refueling  events,  the  distribution  of  bombers  that  success¬ 
fully  refuel  is  binomial. 

PAR^k'BBE  -  V  =  ^  kE)fW1‘lW  E  4,18 

Consider  now  the  case  where  is  greater  than  T^.  In  this  case 
it  is  possible,  although  not  necessarily  inevitable,  that  there  will 
be  more  bombers  needing  refueling  than  tankers  available  to  service 
them. 


If  enough  of  the  tankers  facing  the  base  escape  problem  survive 
to  bring  the  total  number  of  airborne  tankers  up  to  at  least  then 
the  probability  that  k  bombers  get  refueled  will  be  given  by  equation 
4.18  above.  For  this  to  occur,  the  number  of  tankers,  i,  that  survive 
the  base  escape  problem  must  be  greater  than  or  equal  to  the  difference 
between  80E  and  Tg.  The  probability  of  this  is 

TG 

•  PdiBgrV  -  *  PTECi)  4.19 

1”°ar”ts 
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FORWARD  AIR  DEFENSE 

Next  consider  the  problem  of  forward  air  defense.  Temporarily 
set  aside  the  question  of  bombers  penetrating  at  high  altitudes  versus 
penetrating  at  low  altitudes  and  just  look  at  a  homogeneous  force  of 
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force  of  3  bombers  penetrating  through  a  forward  air  defense  zone  de¬ 
fended  by  a  force  of  F  homogeneous  fighters. 

It  is  assumed  that,  with  SUAWACS,  the  defense  can  tell  exactly 
or  make  a  good  estimate  of  the  total  number  of  bombers  in  the  force. 

To  maximize  their  effectiveness,  measured  in  expected  number  of  bombers 
killed,  the  defense  will  attempt  to  assign  a  uniform  number  of  fighters 
to  each  bomber  as  it  enters  the  FAD  zone.  If  F  is  an  integer  multiple 
of  3 , 


F  =  IB  ,  4.22 

then  there  will  be  i  fighters  assigned  to  each  bomber.  Otherwise, 

F  =  iB  +  R  4.23 


where  R  is  some  remaining  number  of  fighters  less  than  B.  In  this 
case,  there  will  be  i+1  fighters  assigned  to  the  first  R  bombers  to 
enter  the  FAD  zone  and  i  fighters  assigned  to  the  other  B-R  bombers. 

Now  consider  one  of  the  bombers  and  the  n  fighters  assigned  to 
it.  Each  fighter  will  have  some  probability  of  detecting  and  con¬ 
verting  on  the  assigned  bomber,  P^,  and  another  probability, 
of  killing  the  bomber  after  detection  and  conversion.  Together, 
these  two  probabil ities  will  give  the  probability  P^  that  an  assigned 
fighter  kills  the  bomber 


Pk] dcPdc 


4.24 


Assume  the  successive  attacks  by  different  fighters  on  the  same 
bomber  are  independent.  This  might  be  because  the  various  fighters 
arrive  in  the  vicinity  and  detect  the  bomber  at  different  times  or 
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simply  because  the  fighter  tactics  against  bombers  do  not  improve  an 
individual  fighter's  probability  of  killing  the  bomber  when  attacking 
in  pairs  or  groups  over  the  probability  when  attacking  by  turn.  If 
this  is  the  case  then  in  order  for  the  bomber  to  survive  the  FAD,  it 
must  survive  the  attacks  of  all  n  fighters  assigned  to  it.  Since  homo¬ 
geneous  fighters  have  been  assumed,  each  will  have  the  same  P^.  There¬ 
fore,  the  bomber's  probability  of  survival,  P  ,  can  be  expressed 

Ps  -  (1  -  Pk)n  4.25 

It  is  assumed  that  each  fighter  can  only  attack  the  bomber  it  was 
assigned  to.  Thus,  there  is  no  link  between  the  killing  of  one  bomber 
and  the  survival  of  another  since  the  fighters  that  have  just  made  the 
kill  cannot  go  to  assist  fighters  attacking  the  other  bomber.  Each 
bomber's  survival  is  therefore  independent  of  what  has  happened  to  the 
other  bombers  after  the  initial  interceptor  assignment.  If  all  the 
bombers  have  the  same  probability  of  surviving,  i.e.,  each  is  attacked 
by  the  same  number  n  of  fighters,  then  the  distribution  of  the  number 
of  bombers  that  survive  is  binomial: 

P(k  bombers  survive  |  B)  =  (® ) P$k( l-?s )B-k  4.26 

It  has  already  been  noted  that  if  F/B  is  non-integer,  there  will 
be  two  distinct  groups  of  bombers,  a  group  of  R  bombers  having  <F/B>+1 
fighters  assigned  to  attack  each  one  and  a  group  of  B-R  bombers  with 
<F/B>  fighters  assigned  to  each.  Within  each  group,  the  probability 
of  survival  for  each  bomber  is  the  same: 

psi  s  Opk)<F/B>+1  4'27a 
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and 


?S2  =  (!-pk)<F/B>  4‘27b 

The  distribution  of  survivors  within  each  group  is  binomial  and  the 
combined  distribution  for  the  total  number  of  bombers  that  survive  from 
both  groups  can  be  found  by  conditioning  on  the  number  that  survive 
from  one  group. 


P(k  survive  FAD  B)  =  TP(k  surv'i  of  R  survive)P(i  of  R  surv) 


=  lP(k-i  of  B-R  surv)P(i  of  R  surv) 
i 

max(k,R)  op  t.  *  R-R-k+i 

z  0PS2 

i=max(Q,k=B+R)  K_1  p 


(i)psi(1"psi)R“1  4-28 

This  completes  the  basic  analysis  of  the  Forward  Air  Defense  pro¬ 
blem.  The  only  question  that  remains  is  how  to  deal  with  two  different 
groups  of  bombers,  those  flying  at  high  altitude  and  those  flying  at 
low  altitude. 

This  model  is  interested  in  the  distribution  of  bombers  from  each 
group  that  survive  the  forward  air  defense,  as  opposed  to  a  distribu¬ 
tion  of  the  total  number  of  bombers  from  both  groups  that  survive. 

Thus,  the  problem  becomes  one  of  deciding  how  to  assign  fighters  to  the 
two  groups.  The  probabilities  a  fighter  detects,  converts  on  and  kills 
a  bomber  may  be  smaller  at  low  altitudes  than  at  high  altitudes  but 
they  will  be  the  same  for  all  fighters  attacking  low  flying  bombers  and 
the  same  for  all  fighters  attacking  high  flying  bombers.  All  the 
independence  assumptions  already  discussed  still  apply.  Knowing  the 
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number  of  bombers  at  each  altitude  and  the  number  of  fighters  assigned 
to  each  group  of  bombers,  one  can  apply  the  above  analysis  to  each 
group  to  find'the  desired  distributions  of  the  numbers  of  high  and 
low  altitude  bombers  surviving  the  FAD. 

It  is  still  assumed  that  the  defense  assigns  fighters  uniformly 
to  all  the  bombers.  Hence,  if  the  number  of  fighters  is  an  integer 
multiple  of  the  number  of  bombers,  then  each  bomber  has  F/B  fighters 
assigned  to  it,  regardless  of  whether  it  is  flying  at  high  or  low 
altitude. 

The  only  question  that  arises  is  how  to  apportion  the  R  remaining 
fighters  if  the  ratio  of  fighters  to  bombers  is  not  an  integer.  For 
the  purposes  of  this  model,  it  is  assumed  that  the  remaining  fighters 
are  assigned,  one  each,  first  to  the  high  altitude  bombers  and  then, 
if  any  fighters  still  remain,  to  the  low  altitude  bombers.  This  is  a 
somewhat  arbitrary  assignment  but  as  justification  consider  that  bombers 
flying  at  high  altitude  will  most  likely  be  detected  earlier  than  those 
flying  at  low  altitudes.  If  the  defense  assigns  the  extra  fighters  to 
the  first  R  bombers  that  are  detected,  then  these  fighters  will 
probably  be  assigned  to  bombers  at  high  altitude. 

BSAM 

After  passing  through  the  forward  air  defense  zone,  the  bombers 
will  reach  the  continent.  The  last  line  of  defense  they  will  encounter 
before  beginning  penetration  is  the  barrier  SAMs. 

Before  approach  the  problem  of  the  barrier  of  SAMs  as  a  whole, 
consider  for  the  moment  a  single  SAM  site  and  a  bomber  penetrating 
through  the  site's  coverage.  The  problem  can  be  mathematically 


60 


described  using  two  variables.  Define  x  to  be  the  lateral  offset  dis¬ 
tance  between  the  path  of  the  bomber  and  the  center  of  the  SAM  site 
(the  distance  between  the  bomber  and  SAM  site  at4the  bomber's  point  of 
closest  approach).  Let  P^(x)  be  the  probability  that  a  bomber  flying 
at  a  specified  altitude  is  killed  by  the  SAM  site  if  the  lateral  off¬ 
set  distance  is  x.  A  sketch  of  the  function  might  look  something  like 

that  shown  in  Figure  4.1,  where  x  is  the  maximum  range  of  the  SAM. 

max 


FIGURE  4.2 

SAM  probability  of  kill  on  bomber  vs  offset  distance. 

For  the  purposes  of  this  model,  the  curve  in  Figure  4.1  is  reduced 
to  two  parameters,  a  probable  lethal  diameter  (PLD)  of  the  SAM  site  and 
a  constant  probability  of  kill,  Pk$,  within  the  PLD.  The  probability  of 
kill  outside  the  PLD  is  assumed  to  be  zero.  (See  Figure  4.3)  Since 
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the  effectiveness  of  the  SAM  may  vary  with  the  altitude  of  the  target. 


the  PLD  and  may  be  different  when  considering  bombers  at  two 
different  altitudes.  Both  the  PLD  and  P,  will  be  inputs  to  the  model, 
one  set  for  each  altitude. 


FIGURE  4.3 

PLD  and  P^  of  SAM  site. 

Now  consider  the  problem  of  a  barrier  of  Ng  SAM  sites,  all  having 
the  same  PLD  and  Pks.  The  sites  are  distributed  uniformly  along  the 
width  of  the  band,  W. 

To  find  the  probability  that  a  bomber  survives  the  barrier  SAMs, 
look  first  at  the  case  where  the  bomber  approaches  a  barrier  with  only 
one  site.  The  site  has  a  width  of  coverage  defined  by  its  PLD  and 
since  it  can  be  located  anywhere  along  the  barrier,  the  bomber’s 


I 


62 


i 


probability  of  crossing  through  the  site's  coverage  is  the  ratio  of 
the  site's  PLD  to  the  total  width  of  the  barrier.  If  the  bomber 
crosses  the  site's  coverage,  the  probability  that  it  is  killed  is  Pks . 
The  probability  that  the  bomber  survives  is  therefore 

Pss  =  l-(w/W)Pk$  4.29 

where  w  is  the  PLO  of  the  SAM  site. 

Now  if  there  are  Ng  identical  SAM  sites  in  the  barrier  and  the 
location  of  each  and  its  P^s  are  independent  of  all  the  other  sites, 
then  the  event  that  the  bomber  survives  any  svte  is  independent  of  its 
event  of  surviving  any  other  site.  The  probability  that  the  bomber 
safely  penetrates  through  the  barrier  is  the  probability  that  it  sur¬ 
vived  all  Ng  sites  which  is 

NB 

PBS  =  [l-(w/W)Pks]  .  4.30 

An  alternative  way  of  deriving  the  above  result  is  to  condition 
the  bomber's  probability  of  survival  on  the  number  of  SAM  sites  it 
encounters.  The  probability  that  the  bomber  encounters  an  individual 
SAM  site  is 

P  =  w/W  .  4.31 

If  there  are  Ng  sites  in  the  barrier,  each  located  independently  of  the 
others,  then  the  distribution  of  the  number  of  sites  encountered  by  the 
bomber  is  binomial 

N  N  -k 

P(encounter  k  sites)  =  (k8)Pk(l-p)  B  .  4.32 
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The  bomber's  probability  of  survival  if  it  encounters  k  sites  is 

k 

P(survives: encounters  k  sites)  =  (1-P^  )  4.33 

To  find  the  probability  that  the  bomber  survives  the  barrier  SAM  defense, 
condition  on  the  number  of  sites  it  encounters 


? 


BS 


Z  ?(survive  encounter  k  sites)P(encounter  k  sites) 
k 


nr  k  nr  k 

=  sB(l-Pks)k(kB)Pk(l-P) 

nr  nr  k 

=  s  (u  HPU-PJnl-p) 

0 

nr 

»  [P(l-Pks)  +  (1-p)] 
nr 

-  [1  -  PFW 

nb 

-  [1  -  (w/W)Pk) 


n-k 


N-k 


4.34 


as  above. 

Now  let  there  be  two  forces  of  bombers,  with  and  denoting 
the  numbers  of  bombers  at  high  and  low  altitude  respectively,  approaching 
the  barrier  of  SAMs.  Assume  that  (a)  within  each  force  each  bomber's 
probability  of  survival  is  the  same,  (b)  each  bomber's  survival  is 
independent  of  the  survival  of  the  other  bombers,  and  (c)  the  forces 
are  independent  of  each  other.  In  this  case,  the  distribution  of  the 
number  of  bombers  that  survive  from  each  force  is  binomial. 


P(n  survive  from  high  forcejB^  in  high  force)  = 

B.  _  Bh-n 

(  n^PBSh  U-W 


4.35a 
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P(m  survive  from  low  force' B  is  low  force)  = 


P3SI(1_PBSZ) 


Vm 


where 

PBSh  -  O-'VWksh1  *  a"d 

pbs;_  •  0-(»,/w)Pks()  8. 


4.35b 
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CHAPTER  V 


AREA  DEFENSES 

In  this  chapter,  the  model  representation  for  the  two  forms  of 
area  defense  will  be  developed.  Specifically,  these  defenses  are  ran¬ 
dom  area  SAMs  (RASAM)  and  fighter-interceptors  (FI). 

In  each  of  the  four  parts  of  the  model  developed  in  the  previous 
chapter,  it  was  only  necessary  to  know  for  each  bomber  whether  it 
"succeeded"  or  "failed"  that  aspect  of  its  mission.  In  the  case  of  FAD, 
for  example,  the  bomber  either  survived  the  entire  engagement  or  was 
killed;  if  it  was  killed,  it  did  not  matter  how  long  it  had  survived 
or  how  many  intercepts  it  had  survived  before  it  was  finally  shot  down. 

The  situation  is  different  in  the  case  of  area  defense.  Mow  the 
bomber  is  at  the  stage  where  it  will  begin  to  release  its  weapons. 

Each  bomber  has  been  assumed  to  carry  more  than  one  weapon  and  the 
bomber's  targets  are  located  at  different  depths  along  the  bomber's 
penetration  route  into  the  defended  area.  The  bomber  need  not  survive 
all  the  way  through  the  defended  area  in  order  to  release  some  or  maybe 
all  of  his  weapons.  Rather,  the  number  of  weapons  he  releases,  and 
hence  the  damage  he  extracts,  will  depend  on  how  far  into  the  defended 
area  the  bomber  successfully  penetrates. 

For  this  reason,  it  is  desired  to  find  a  probability  of  survival 
for  a  bomber  that  is  a  function  of  the  bomber’s  depth  of  penetration. 
This  penetration  probability  function  will  be  denoted  Pp(x)  and  will 
incorporate  the  bomber's  probabilities  of  successfully  penetrating  to 
at  least  a  depth  of  x  through  both  RASAM  and  FI  defenses. 


66 


In  this  model,  it  is  assumed  that  the  RASAM  and  FI  defenses  operate 
independently  of  each  other.  This  allows  each  defense  to  be  modelled 
separately  and  the  results  combined  to  give  the  desired  penetration 
probability  function.  The  remainder  of  this  chapter  describes  the 
modelling  of  the  RASAM  and  FI  engagements  individually  and  finally 
derives  the  penetration  probability  function. 

RASAM 

The  approach  used  to  model  RASAMs  is  very  similar  to  that  employed 
for  the  BSAMs.  The  major  difference  is  that  now  the  problem  has  depth 
and  it  is  desired  to  find  a  probability  of  survival  as  a  function  of 
depth  of  penetration. 

Consider  the  case  of  a  single  bomber  flying  a  course  due  south 
down  the  corridor.  As  in  the  case  of  the  BSAM  problem,  first  look  at 
the  case  where  the  entire  corridor  is  defended  by  just  one  SAM  site. 

The  SAM  site  will  be  modelled  using  a  probable  lethal  diameter  (PLD) 
and  a  probability  of  kill,  Pks>  within  the  PLD,  consistent  with  the 
description  in  the  BSAM  section.  The  PLD  and  Pks  need  not  be  the  same 
for  RASAMs  as  they  were  for  BSAMs. 

Three  possible  subcases  of  the  scenario  just  described  are  illus¬ 
trated  in  Figure  5.1,  showing  the  defended  territory,  the  bomber  course 
and  the  SAM  site.  In  Figure  5.1,  the  dashed  lines  represent  a  distance 
to  either  side  of  the  bomber  path  equal  to  half  the  PLD  of  the  SAM  site. 

The  bomber's  path  will  cross  through  the  coverage  of  the  SAM  site 
only  if  the  center  of  the  SAM  site  is  located  between  the  two  dashed 
lines  in  Figure  5.1.  In  other  words,  for  a  SAM  site  of  coverage  width 
w=PLD,  to  threaten  the  bomber,  the  site  must  be  centered  at  a  distance 
no  greater  than  ^w  from  either  side  of  the  bomber  path.  Furthermore, 
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no  sites  south  of  the  distance  x  which  the  bomber  has  penetrated  to 
will  have  yet  threatened  the  bomber. 


FIGURE  5.1 

Bomber  vs.  random  area  SAMs. 


Now  consider  the  probability  that  the  bomber  crossed  through  the 
coverage  of  the  SAM  site  while  penetrating  to  a  depth  x  into  the 
defended  territory.  The  bomber  will  have  passed  through  the  coverage 
of  the  SAM  site  if  the  center  of  the  SAM  site  is  located  in  the  area 
between  the  two  dashed  lines.  This  area  is  A=wx.  The  total  area  in 
which  the  SAM  site  may  be  found  is  just  the  total  defended  area  which 
is  Since  it  is  assumed  that  the  SAM  is  located  anywhere  in  Ay 

with  uniform  probability,  the  probability  the  SAM  engages  the  bomber, 
P,  is  the  ratio  of  the  two  areas  or 
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wx 

P  =  WD 


5.1 


If  the  bomber  crosses  through  the  coverage  of  the  SAM  site,  its 
probability  of  being  killed  is  Pks.  Therefore,  the  probability  that 
the  bomber  survives  to  at  least  a  depth  x  against  one  P.ASAM  (Pl(x))  is 


O'  Cv  )  =  i_pp  =  i-p  — 
, x,x)  i  ^Kks  l  Kks  WD 


5.2 


Note  that  the  probability  that  the  bomber  penetrates  all  the  way  through 
the  defended  area  is 


ps(D)  •  1_pks  if  5-3 

which  is  the  same  result  derived  for  the  BSAM  problem. 

Now  consider  the  case  where  there  are  NR  identical  RASAMs  defending 
the  country.  Assume  (a)  each  of  the  SAMs  is  located  independently  of 
the  others  and  (b)  the  engagements  between  the  SAMs  and  bombers  are 
independent.  Therefore,  if  the  bomber  is  to  survive  to  a  depth  x,  he 
must  survive  independent  threats  from  all  NR  SAMs.  His  probability  of 
surviving  to  a  depth  x,  against  SAM  defense  is  therefore 

ps<*>  ■  «-pks  ro>"R  5-4 

Two  implied  assumptions  should  be  noted  about  the  above  formulation 
of  the  RASAM  problem.  When  deriving  the  probability  that  a  bomber  would 
cross  through  the  coverage  of  a  given  SAM  site,  the  area,  wx,  within 
which  the  SAM  site  had  to  be  located  in  order  to  threaten  the  bomber 
was  implicitly  assumed  to  be  entirely  contained  within  the  total  defended 
area,  WO.  However,  if  the  bomber's  course  is  close  to  either  border  of 
the  defended  area,  this  may  not  be  the  case.  Any  bomber  course  less 
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than  a  distance  ^w  inside  the  boundary  would  place  part  of  the  area 
wx  outside  the  defended  territory.  For  such  a  bomber,  the  probability 
of  crossing  the  coverage  of  a  given  SAM  site  would  be  slightly  less  than 
the  probability  given  in  equation  5.1.  This  is  an  improbable  scenario; 
most  targets  will  be  more  centrally  located  and  the  bomber's  course  will 
therefore  not  be  that  close  to  the  boundary  of  the  defended  area. 

Also,  equation  5.1  implicitly  assumes  that  a  SAM  site  does  not 
threaten  a  bomber  until  the  bomber  has  reached  a  point  even  with  the 
center  of  the  SAM  site.  In  principle,  since  SAM  footprints  are  two- 
dimensional,  it  would  be  possible  for  the  bomber  to  be  killed  by  a  SAM 
centered  near  the  bomber's  projected  path  and  slightly  south  of  the  dis¬ 
tance  x,  roughly  between  x  and  x+4w.  To  correct  Equation  5.1  to  account 
for  this  wu  '' ^  require  making  further  assumptions  about  SAM  coverage. 

For  values  o.‘  ♦ ,  W,  and  D  much  larger  than  w,  the  practical  effect  of 
such  a  correction  would  be  negligible. 

FIGHTER-INTERCEPTOR 

At  first  glance,  it  might  appear  that  the  fighter-bomber  engage¬ 
ment  problem  over  land  might  be  modelled  similarly  to  the  FAD  problem. 
However,  several  things  prevent  such  an  approach.  First,  the  model  is 
concerned  with  finding  a  probability  of  survival  as  a  function  of  the 
depth  of  penetration.  There  was  no  problem  extending  the  BSAM  model 
to  include  he  depth  dimension,  thereby  arriving  at  a  model  for  RASAMs. 
However,  unlike  SAM  sites,  fighters  are  mobile  and  can  cover  much  larger 
areas  so  the  same  method  cannot  be  employed  to  extend  the  FAD  model  into 
an  area  FI  defense  model . 

Another  factor  preventing  any  easy  modification  of  the  FAD  approach 
is  the  problem  of  determining  the  number  of  intercepts  that  can  be  made 
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on  each  bomber.  It  is  no  longer  reasonable  to  assume  that  each  fighter 
will  have  only  one  chance  to  make  an  intercept.  If  the  bombers  are 
penetrating  to  depths  measured  in  thousands  of  kilometers,  the  engage¬ 
ment  will  last  long  enough  to  allow  fighters  to  make  an  intercept, 
return  to  the  base,  refuel,  and  then  repeat  the  process,  possibly 
several  times. 

For  these  reasons,  an  entirely  different  approach  is  used  to  model 
the  bomber-FI  engagement.  The  method  is  essentially  that  used  in  the 
SRI  CGPEM-1  model  and  is  described  below  (Ref  5:27-33). 

Before  beginning  the  actual  derivation,  some  variables  must  be 
defined . 

P(j(t'|td)  will  denote  the  probability  that  a  bomber  that  is 
detected  at  some  time  t^  after  beginning  penetration  survives  an  addi¬ 
tional  time  of  at  least  t' . 

P i ( t ‘ t  td )  will  be  defined  as  the  probability  that  exactly  i  of  the 
fighters  that  are  assigned  to  the  bomber  detected  at  time,  t^,  are  able 
to  intercept  the  bomber  within  a  time  t'  after  the  time  of  detection 
t^.  This  does  not  imply  that  all  i  fighters  will  successfully  detect 
and  convert  on  the  bomber,  only  that  the  i  fighters  are  close  enough 
to  the  projected  course  of  the  bomber  that  they  can  be  vectored  into 
the  vicinity  of  the  bomber  within  time  t'. 

The  probability  that  a  fighter  successfully  detects,  converts  on, 
and  kills  the  bomber  it  has  been  assigned  to  intercept  will  be  denoted 

her 

II 

Pp(t)  will  be  probability  that  the  bomber  survives  to  a  time  t 
after  initially  entering  the  corridor  and  beginning  penetration.  (The 
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primes  serve  only  to  distinguish  this  function  from  other  functions 
which  will  also  be  labelled  Pp) . 

The  probability  distribution  function  (pdf)  for  the  time  after 
entering  the  corridor  that  the  bomber  is  first  detected  by  the  radar 
net  is  Pd(td). 

The  pdf  for  the  time  at  which  a  bomber  first  enters  the  corridor 
wil 1  be  given  by  p  (t  ) . 

d  d 

Note  the  difference  between  t.  and  t  .  The  time  of  detection,  t,, 

da  d 

is  measured  for  each  bomber  from  the  time  the  bomber  enters  the  corridor 
and  does  not  specify  a  point  in  time  of  the  engagement.  The  arrival 
time  at  the  corridor  entrance  however  is  measured  from  some  standard 
time  0  for  all  bombers.  If  the  engagement  is  said  to  begin  at  time  0, 
when  the  first  bomber  begins  penetration,  then  a  subsequent  bomber  will 
be  said  to  begin  penetration  at  some  time,  ta  and  be  detected  at  some 

a 

time  td  later,  in  other  words  at  time  ta+td. 

Consider  first  P$( t ' | t^) ,  the  probability  that  the  bomber  survives 
at  least  an  additional  time  t'  after  being  detected  at  td.  This 
probability  can  be  conditioned  on  the  number  of  intercepts  that  occur 
within  the  time  t'  as  follows: 

PF(f|td>  =^oa“Pkf)iPi(t,'td)  5,5 

Assume  temporarily  that  the  bomber  begins  penetrating  at  time  t=0. 

•« 

Pp(t),  the  probability  that  a  bomber  survives  to  time  t,  can  be  found 
by  conditioning  on  the  time  of  detection,  td. 

PF(t>  *  't  Pstt-tdltd)Pd(td>dtd  *  5-6 
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The  first  integral  expresses  the  probability  that  the  bomber  is 
detected  at  some  time  t^  less  than  t  and  then  survives  an  additional 
time  t-tj.  The  second  integral  gives  the  probability  that  the  bomber 
is  not  detected  until  some  time  greater  than  t,  in  which  case  it  is 
assumed  that  he  is  guaranteed  of  surviving  to  time  t.  Equation  5.5 
can  be  substituted  into  5.6  to  give 


PF(t) 


.t 
■  o 


<1-h.f>’pi<t-td!td>Pd(t<jldVl-  \'o  »d(td)dtd 


1=0 


}  -r0  P,(t-td;td)Pd(td)dtdrt-Pd!q  5.7 


1=0 


where  Prf(t)  is  the  cumulative  distribution  function  (CDF)'  of  pd(td) 
and  the  functions  are  all  assumed  to  be  well  behaved  so  that  the  inte¬ 
gration  and  summation  can  be  exchanged. 

Equation  5.7  gives  a  bomber's  probability  of  survival  to  at  least 
a  time  t.  The  quantity  tyiat  is  of  interest  to  the  model  is  a  function 
Pp(x)  which  gives  the  probability  of  survival  for  the  bomber  as  a  func 
tion  of  penetration  depth.  This  can  be  found  by  assuming  that  the 
bombers  fly  at  a  constant  velocity  v.  In  a  time  t  then,  the  bomber 
will  travel  a  distance  x=vt.  Put  differently,  in  order  to  survive  to 
a  distance  x,  a  bomber  must  survive  a  time  t=x/v.  Thus, 


PI r(x) 


PFCx/v)  = 


z 

i=o 


-tdltd>pd(td)atc 


+  n-pd(|-)]  5.8 

The  last  thing  to  account  for  is  the  possibility  that  the  bomber 
enters  the  corridor  at  some  time  greater  than  0.  It  is  assumed  that 
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the  probability  of  detection  depends  only  on  the  depth  of  penetration 
and  not  on  the  time  of  arrival  at  the  corridor  entrance.  The  function 
P.j(t';t),  however,  may  in  general  depend  on  the  absolute  time  at  which 
detection  occurs.  For  example,  P . ( t ' ! t ) ,  the  probability  that  i  inter¬ 
cepts  could  take  place  within  time  t’  after  detection  of  the  bomber  at 
time  t,  may  depend  on  either  the  number  of  fighters  or  the  number  of 
bombers  in  the  air  at  time  t.  These  numbers  would  both  vary  with  abso¬ 
lute  time,  not  with  the  time  of  detection  t^  of  a  bomber.  Thus, 

Pi(t* | t)  =  p.( t ■ ; ta+td) 
instead  of  simply 

Pi  (t  • ;  t)  =  pi(f  |td) 


By  making  this  one  change  in  Equation  5.8,  one  can  find  a  function 
PF(xita)  which  is  the  probability  a  bomber  survives  to  penetrate  at 
least  a  distance  x  given  that  he  begins  penetration  at  time  t  ; 

a 


PF^X  1  ta)  =  .E  ^1_PkP  !o  *Vv  ^d'  ta+td'lpd^td^dtd 
i  =o 


+  U  -  PdC^-)3 


5.9 


Typically,  one  is  not  concerned  with  the  order  in  which  bombers 
penetrate.  The  unconditional  average  penetration  probability,  Pp(x),  is 

PpU)  =  fQ  PF^-Xita^pa^ta^dta 


=  ;o  {iE0^1'Pkf^  ;o/VpiC  VVpd(VdV[1*Pd(v^}pa(ta)dta 
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=  «-pkf>1  -i  -f V  Pi<7  -VVtd>‘,d<td>',a(t.>dt.dt. 

1  =0 


+  [1-Pd(v)]  5-10 

Once  again,  all  the  functions  have  been  assumed  to  be  well  behaved 
so  that  the  integral  and  the  summation  can  be  exchanged. 

This  is  the  general  form  of  the  penetration  probability  function 
derived  by  de  Sobrino,  Monahan,  Weinstein,  and  Yu.  (Ref  5:27-30).  The 
next  step  is  to  identify  choices  for  the  functions  P.(t';  t),  P^t.) 
and  pa(ta). 

Consider  first  the  function  P.(t'|t)  which  is  the  probability 
density  function  for  the  number  of  possible  intercepts  that  may  be  made 
on  a  bomber  within  a  time  t'  after  detection  at  time  t.  One  way  of 
looking  at  the  problem  is  to  consider  a  counting  function  N(t')  which 
counts  the  number  of  events,  in  this  case  intercepts,  that  occur  in 
the  time  interval  [0 , t ' ] . 

Since  it  has  been  assumed  that  bombers  are  detected  by  a  ground 
based  radar  net  and  not  by  chance  encounters  with  airborne  interceptors, 
it  is  reasonable  to  assume  that 

N ( 0)  =  0  5.11 

In  other  words  no  intercepts  will  have  been  made  at  the  time  of  detec¬ 
tion. 

In  addition,  assume  that  events  cannot  occur  simultaneously;  some 
time,  no  matter  how  small,  must  elapse  between  intercepts.  Mathemati¬ 
cally, 
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5.12 


,  P[N(t+h)-N(t)>2] 

For  any  t>0,  -  =  0 

h-o  P[N(t+h)-N(t)=l] 

Consider  now  a  third  assumption: 

For  any  t>0,  0<P[N{t)>0]<l  5.13 

This  assumption  implies  that  in  any  time  interval,  regardless  how 
small,  there  is  a  finite  probability  that  at  least  one  intercept  could 
occur  but  there  is  no  certainty  that  an  event  will  occur. 

Detection  of  a  bomber  does  not  guarantee  that  there  will  be  an 
unassigned  fighter  available  to  intercept  the  bomber.  There  is  a 
chance,  in  other  words,  that  even  though  the  bomber  is  detected,  there 
will  be  no  intercepts  made  on  it.  The  upper  inequality  in  Equation 
5.13  therefore  seems  reasonable. 

The  lower  inequality  poses  more  of  a  problem.  Physical  limitations 
of  the  command  and  control  system  will  cause  some  time  lag,  maybe  quite 
small,  between  the  actual  detection  of  the  bomber  and  assignment  of  an 
interceptor  to  it.  Even  if  instantaneous  communications  were  possible, 
some  time  would  still  have  to  elapse  to  allow  the  fighter  to  reach  its 
interception  point  with  the  bomber.  Therefore,  it  seems  likely  that 
P [N ( t ) >0] =0  for  some  small  values  of  t,0<t<e,  contrary  to  the  assump¬ 
tion  in  Equation  5.13.  Nevertheless,  allowing  extremely  fast  communi¬ 
cations  may  not  be  too  unreasonable  and  one  can  imagine  a  case  where 
an  interceptor  would  be  stationed  quite  close  to  the  place  where  the 
bomber  was  detected.  The  e,  which  would  describe  the  absolute  minimum 
time  that  would  elapse  before  the  first  intercept  was  possible,  might 
then  be  quite  small.  Assuming  that  e=0  would  give  the  defense  only  a 
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very  small  advantage  and  should  not  signi ficantly  affect  the  results 
of  the  model . 

Equations  5.11-13  are  three  of  the  four  assumptions  that  define 
a  Poisson  process.  The  fourth  assumption  states  that  the  process  has 
independent  increments.  Mathematically ,  this  says  that  for  all  choices 
of  ti  such  that  0<t^<t?. . .  .<tn>  the  n  random  variables  )-«%'( )  ] , 

[N(t2)-N(t1)J . .  [N(tn)-N(tn_j)j  are  independently  distributed. 

In  this  model  the  above  assumption  would  imply  that  the  distribu¬ 
tion  of  the  number  of  intercepts  that  could  occur  in  (t^,  t^)  would 
be  independent  of  the  number  that  had  occurred  in  (C,t^).  This  assump¬ 
tion  is  clearly  not  very  sound  for  small  numbers  of  interceptors. 
Imagine  the  case  where  there  is  only  one  fighter.  If  it  made  an  inter¬ 
cept  in  the  interval  (0,t^)  then  it  might  be  impossible  for  it  to  make 
another  intercept  in  (t^tg)  because  it  had  to  return  to  base.  This 
would  violate  the  assumption.  However,  as  the  number  of  interceptors 
increases,  this  becomes  less  of  a  problem.  Thus,  for  reasonably  large 
numbers  of  interceptors  the  assumption  is  probably  acceptable.  If  the 
number  of  interceptors  is  not  too  large,  which  might  be  the  case  if 
the  bomber  strike  were  preceded  by  a  nuclear  missile  strike,  then  this 
assumption  may  be  very  offense  conservative. 

If  all  four  of  these  assumptions  are  accepted,  keeping  in  mind 
their  limitations,  then  one  can  claim  that  the  intercept  probability, 
P.(t‘|t),  may  be  approximated  with  a  Poisson  distribution: 


P^t’ 


t)  * 


utt 


i  -xtr 
)’e  1  /i! 


0,  1,  2, 


5.14 


This  is  the  function  used  by  deSobrino,  Monahan,  Weinstein  and  Yu 
in  the  C0PEM-1  model  (Ref  5:30-33).  The  paramter  is  called  the 
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intercept  intensity  function  and  is  subscripted  to  indicated  that  it 


may  depend  on  time.  For  a  bomber  detected  at  time  t,  is  the 
average  number  of  intercepts  possible  against  the  bomber  per  unit 
time.  In  a  time  t‘  after  detection  of  the  bomber,  the  expected  number 
of  possible  intercepts  is 

Substituting  Equation  5.14  into  Equation  5.10,  one  has 


P  f  ( x ) 


i  =o 


u-hcf)1  -l  T  'ivy! 


td)'V1(ta*td)(rtd) 


i! 


Pd(td^a{ta)dtadta  +  [1‘Pd^^  5-15 

By  re-exchanging  the  summation  with  both  integrals  and  using  the 

identity 

co  i 

r  fr=ex  5.16 

i  =o  ' 

Equation  5.15  can  be  reduced  to 

PF(X)  m’l  !T  td)Pd<td>l’.(t.)dtddta 

*  5.17 

Leaving  aside  the  pdfs  of  t  and  td  temporarily,  Equation  5.17 
still  has  the  problem  of  a  choice  for  the  function  X^..  This  function 
will  in  general  depend  on  such  factors  as  the  number  and  location  of 
airborne,  unassigned  interceptors,  interceptor  characteristics  such  as 
velocity,  and  the  number  of  penetrators  in  the  defended  area.  At  any 
time,  the  function  will  also  depend  on  the  events  that  have  occurred 
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earlier  in  the  engagement,  de  Sobrino  et.al.  have  developed  an  itera¬ 
tive  routine  that  can  be  used  to  approximate  this  function.  The  con¬ 
cept  is  based  on  properties  of  the  Poisson  distribution.  For  a  Poisson 
distribution  with  parameter  A,  the  waiting  time  to  the  first  event  is 
distributed  exponentially  with  the  same  parameter  \ .  The  expected 
value  of  the  exponential  deistribution  is  i/v .  Therefore,  if  an  esti¬ 
mate  can  be  found  for  the  expected  waiting  time  to  the  first  intercept 
of  a  bomber  detected  at  time  t,  the  reciprocal  of  this  estimate  will 
in  turn  be  an  estimate  of  the  parameter  A  (Ref  5:33). 

A  description  of  the  calculation  to  find  the  expected  waiting 
time  to  the  first  intercept  can  be  found  in  references  4  and  8. 

This  calculation  can  be  performed  for  several  values  of  t  and  then 
appropriate  interpolation  procedures  can  be  employed  to  give  an 
estimate  of  At  for  other  times.  This  estimate  can  then  be  used  to 
numerically  integrate  Equation  5.17  (assuming  suitable  functions  or 
estimates  of  Pa(ta)  and  Pd(td)  have  been  specified). 

This  approach  probably  gives  the  best,  in  the  sense  of  most 
intuitively  appealing,  results.  However,  if  just  two  additional 
assumptions  are  made,  a  vast  simplification  results  on  Equation  5.17. 

First  consider  what  happens  if  the  parameter  A  is  assumed  to  be 
independent  of  time.  In  this  case,  one  has 

+  tl-Pd(7)]  5.18 

The  integrand  of  Equation  5.18  no  longer  contains  any  terms  which 
depend  on  t  except  the  pdf  of  t. .  Since  the  integral  of  a  pdf  over 

a  a 

its  full  range  is  1,  one  has 
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5.19 


PF(X)  •  :lh  exp[-Pkfi(i-  td)]pd(td)dta*[l-Pd(i)] 

With  x  assumed  to  be  a  (still  unspecified)  constant,  the  only  func¬ 
tion  of  time  remaining  in  the  integral  is  the  pdf  of  t^,  the  distri¬ 
bution  of  the  time  of  detection  of  the  bomber.  Consider  this  function 
now. 

It  has  already  been  assumed  that  the  bombers'  probability  of  detec¬ 
tion  is  a  function  only  of  its  altitude  and  depth  of  penetration. 
Temporarily  assume  that  all  the  bombers  are  flying  at  the  same  constant 
altitude.  Now  one  is  left  with  a  function  depending  only  on  the  depth 
of  penetration  which,  since  the  bombers  are  also  assumed  to  be  flying 
at  a  constant  velocity,  is  equivalent  to  being  a  function  of  penetration 
time. 

Assume  that  the  defensive  ground  based  detection  system  uniformly 
covers  the  defended  area.  If  this  is  the  case,  then  it  is  reasonable 
to  believe  that  an  undetected  bomber's  chance  of  being  detected  while 
flying  a  given  distance,  or  equivalently  a  given  time,  would  be  the 
same  for  all  distances  of  the  same  length  no  matter  where  they  were 
along  the  bomber's  path.  For  example,  consider  a  bomber  that  has 
penetrated  undetected  to  a  depth  of  1,000  kilometers.  His  probability 
of  escaping  detection  for  the  next  100  kilometers  will  be  the  same  as 
his  original  probability  of  penetrating  the  first  100  kilometers 
undetected . 

Mathematically,  using  time  instead  of  distance,  the  above  can  be 
expressed  as 
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p[(td>t1+t2) : td>t1]  -  p(td>t2) 


5.20 


where  td  is  the  time  of  detection  after  penetration  has  begun. 

Equation  5.20  is  a  property  of  the  exponential  distribution.  If 
one  is  willing  to  assume  that  the  distribution  of  the  time  to  detection 
of  a  bomber  is  characteri zed  by  the  above  assumptions,  then, 

-0 1 1  , 

Pd  (I'd)  =  ae  a>0,  td— ^  5.21 


The  parameter  a  is  a  measure  of  the  capability  of  the  radar  net  to 
detect  the  bomber.  Specifically,  1/a  is  the  expected  time  to  detec¬ 
tion  of  the  bomber.  An  estimate  of  the  expected  time  to  first  detection 
should  consider  the  intensity  of  the  defensive  radar  coverage  and  the 
effectiveness  of  bomber  tactics  such  as  low  level  penetration. 

Substituting  Equation  5.22  into  Equation  5.18,  one  has 


PFtx)  =  /*/v  exp[-Pdfx(|  -  td)]oe  atddtd+[l-  f*/v ae  5.22 

where  the  CDF  of  td  has  been  expressed  as  the  integral  of  the  pdf  in 
the  last  term. 

The  above  function  can  now  be  integrated  directly  to  give,  after 
combining  terms. 


PFCx)  *  Lf*- 


kf 


A -a 


■)e 


x 

-°y 


kPkfA-a 


)e 


-P  A— 
KkfAv 
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This  function  has  the  major  advantage  of  simplicity.  Mo  numerical 
integrations  are  required.  A  bomber's  probability  of  penetrating  to 
at  least  a  distance  x  Inside  the  defended  area  is  now  expressed  as  a 
simple  function  of  its  velocity  and  three  other  parameters.  Because 
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of  this  simplicity,  this  is  the  form  of  the  penetration  probability 
against  FI  defenses  that  is  used  in  this  model. 

As  an  input,  the  paramter  P^,  the  probability  that  a  fighter 
kills  the  bomber  it  is  assigned  to,  can  be  estimated  using  the  charac¬ 
teristics  of  the  aircraft  and  models  of  the  engagement  such  as  COLLIDE 
(Ref  1:8-5).  Reasonable  choices  for  the  parameters  a  and  X,  however, 
present  greater  difficulties.  One  approach  to  the  problem  of  esti¬ 
mating  x  has  already  been  referenced.  It  would  seem  reasonable  that  a 
similar  method  could  be  developed  to  estimate  a.  Explicit  techniques 
for  estimating  these  parameters  is  a  suggested  future  research  topic. 

For  the  purposes  of  this  thesis,  both  a  and  x,  as  well  as  P^p,  will  be 
taken  to  be  externally  provided  to  the  model. 

PENETRATION  PROBABILITY  FUNCTION 

The  previous  two  sections  have  derived  functions,  Ps(x)  and  Pp(x), 
giving  a  bomber's  probability  of  surviving  to  at  least  a  depth  x 
against  either  SAM  defenses  or  FI  defenses  individually.  The  only 
thing  that  remains  to  do  is  to  combine  these  results  into  a  single 
penetration  probability  function  Pp(x). 

In  order  for  a  bomber  to  sucessfully  penetrate  to  a  depth  x,  he 
must  survive  the  threats  from  both  forms  of  defense.  It  has  been 
assumed  that  the  two  defenses  effect  the  bomber's  survival  independently. 
Thus,  his  probability  of  surviving  both  defenses  is  just  the  product  of 
his  probability  of  surviving  each  defense  individually.  Therefore,  one 
has 


P p(x)  =  Ps(x)PF(x) 


5.24 
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Substituting  from  Equations  5.4  and  5.23  for  the  functions  P  (x) 


and  P.(x),  one  gets 


N-  P.  rA  -a- 

p  (x\  =  ri.p  wx,  Rf(_  kf__)e  v  _  , 
V  pks  WDJ  aP,  x-a;e  1 


V-* 


)e 


-P. 
kf  v. 
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This  is  the  final  form  of  the  penetration  probability  function. 

Up  to  this  point,  no  distinction  has  been  made  in  the  derivation 
between  bombers  flying  at  high  and  low  altitudes.  The  reasoning  will 
be  identical  in  both  cases  and  the  functional  form  of  the  penetration 
probability  function  will  thus  be  the  same  for  both.  Various  performance 
parameters  (Pksi  Pkf>  w,  v,  a  and  possibly  \)  may  have  different  values 
when  the  bombers  penetrate  at  high  altitudes  than  when  bombers  penetrate 
at  low  altitudes.  Therefore,  the  model  will  use  two  different  probab¬ 
ility  penetration  functions  corresponding  to  the  two  sets  of  values 
for  the  parameters. 
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CHAPTER  VI 


MISSION  COMPLETION  AND  CRUISE  MISSILES 

In  this  chapter,  the  model  descriptions  for  the  last  two  stages 
of  the  bomber  mission  are  developed.  Specifically,  these  stages  are 
weapons  delivered,  including  the  terminal  SAM  defenses,  and  recovery. 

Also  in  this  chapter  the  model  representations  for  cruise  missiles 
(CM)  are  described.  Three  different  cases  will  be  covered  corresponding 
to  CMs  of  three  different  ranges. 

TERMINAL  SAMs  AND  WEAPONS  DELIVERED 

The  next  problem  to  be  addressed  is  the  modelling  of  terminal  SAM 
engagements  and  the  number  of  weapons  delivered. 

Consider  first  the  case  of  one  bomber  and  its  load  of  N  weapons. 

It  has  already  been  assumed  that  the  weapons  are  homogeneous  and  that 
only  one  weapon  will  be  expended  on  each  target.  The  bomber's  targets 
are  assumed  to  have  been  preselected  as  part  of  his  mission  plan  and 
are  located  at  various  distances  along  his  course  into  the  defended 
area.  The  targets  are  all  assumed  to  be  uniformly  defended  by  some 
number,  possibly  zero,  of  terminal  SAMs. 

Label  the  points  at  which  the  bomber's  targets  are  located  Xj, 

X£,  .  xn,  where  o  <  x^  <  Xg  <  .  <  xn  <  D.  The  distances  are 

measured  from  where  the  bomber  began  penetration  and  D  is  the  total 
depth  of  the  defended  area. 

Let  n  be  the  number  of  weapons  released.  The  probability  that  at 
least  one  weapon  is  released  may  be  written  P(n  >_  1).  For  one  weapon 
to  be  released,  the  bomber  must  have  penetrated  to  at  least  the  depth 
of  his  first  target  x^.  This  probability  that  the  bomber  survives  as  a 
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function  of  the  depth  of  penetration  has  been  derived  in  Chapter  V  and 
is  given  by  the  penetration  probability  function  P  (x^). 

Besides  reaching  the  target  through  the  area  defenses,  it  is 
assumed  that  the  bomber  must  survive  the  terminal  SAM  defenses  before 
he  can  release  a  weapon  at  that  target.  Denote  this  probability  of 
survival  as  Py.  Since  the  targets  are  assumed  to  be  uniformly  defended 
Py  will  be  a  constant  in  the  model  for  bombers  flying  at  a  given  alti¬ 
tude.  Py  need  not  be  the  same,  however,  for  bombers,  at  different 
altitudes.  Note  also  that  Py  does  not  have  to  be  the  same  as  the 
bomber's  chance  of  surviving  a  SAM  attack  in  general.  Several  different 
assumptions  can  be  made  about  terminal  SAM  defenses  and  the  offensive 
tactics  used  against  them.  As  just  one  example,  consider  targets 
defended  by  fixed  SAM  sites.  The  bombers  may  attack  the  SAM  sites  with 
SRAMs  before  moving  in  closer  to  release  a  weapon  at  the  target  itself. 
In  this  scenario,  Py  would  be  the  same  as  the  probability  that  the 
bomber  successfully  destroys  the  SAM  site  or  fails  to  destroy  it  but 
survives  the  SAM  attack.  In  the  case  of  undefended  targets  Py  is  set 
to  1 . 

From  the  above,  the  probability  that  at  least  one  weapon  is 
released  is 

p(n  >  1)  *  Ps(Xl)Py  .  6.1 

The  probability  that  at  least  two  weapons  are  released  can  be 
found  in  the  same  way.  To  have  released  two  weapons,  the  bomber  must 
penetrate  to  at  least  a  depth  x^  and  must  have  survived  two  terminal 
SAM  engagements,  one  at  each  target.  Hence, 
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5.2 


P(n  >  2)  =  PjxJP; 

—  S  c  i 

In  general ,  one  fi nds 

P(n  >  i)  =  Ps(x.)pj  i  =  0,  M;  xq  h  o  6.3 

Mow  consider  the  probability  that  exactly  one  weapon  is  released. 
Since  the  number  of  weapons  released  is  a  discrete  distribution  one 
has 

P(n  =  1)  =  P(n  >  1)  -  P(n  >  1) 

=  P(n  >_  1)  -  P(n  >  2)  6  A 

Therefore  the  probability  that  exactly  one  weapon  is  released 
can  be  expressed 

P(n  =  1)  =  -  Ps(x2)PT  6-5 

and  in  general  for  exactly  i  weapons,  one  has 

P(n  =  i)  =  Ps(Xi -  Ps(xi+1)PT+1  i  =  °*  N_1  6-6a 

PCn  =  N)  =  Ps(xN)Py  6.6b 

Since  the  bomber  only  carries  N  weapons  the  probability  of 
releasing  more  than  N  is  zero  or,  to  put  it  differently,  the  probability 
of  releasing  exactly  N  is  the  same  as  the  probability  of  releasing  at 
least  N.  Thus  the  difference  between  equations  6.6a  and  6.6b. 

Equations  6.6  give  the  distribution  for  the  number  of  weapons 
released  by  each  bomber.  For  a  given  number  of  bombers  that  survived 
the  BSAM  defense  to  begin  penetration,  it  would  in  principle  be 


86 


possible  to  find  a  distribution  for  the  total  number  of  weapons  released 
by  all  (say  B)  bombers.  This  would  be  given  by 
N  N  N 

?(nt=k)  =.  Z  z  Z  P_(n1=i 1)P(n2=i2) - p(nB.i-iB_1) 

ir°  i2=oiB-r0 

Where  nfc  is  the  total  number  of  weapons  releases  and  n.  is  the  number 

of  weapons  released  from  the  ith  bomber  and  P(nt  <  0)  =  ?(nt  >  11)  =  0. 

For  large  numbers  of  surviving  bombers,  this  would  be  a  very  unweildy 
sum.  An  approximate  alternative  method  is  outlined  below. 

Consider  a  single  arbitrary  weapon  out  of  the  N  weapons  on  a 
particular  penetrating  bomber.  Since  this  weapon  is  identical  to  all 
the  others  on  the  bomber,  it  will  not  be  assigned  to  a  particular 
target,  rather  it  will  have  equal  probability  of  being  released  over 
any  of  the  targets.  Another  way  of  looking  at  that  is  that  as  the 
bomber  reaches  any  given  target,  that  weapon  is  as  likely  to  be 
released  as  any  other  weapon. 

If  the  bomber  survives  to  release  a  weapon  over  exactly  one  tar¬ 
get,  the  weapon  will  have  a  chance  P1  =  1/N  of  being  the  one  that  is 
released.  Similarly,  if  the  bomber  survives  to  attack  exactly  two 
targets,  the  weapon  will  have  a  probability  P^  =  2/N  of  being  one  of 
the  two  weapons  released  and  in  general  the  weapon’s  probability  of 
being  released  if  the  bomber  succeeds  in  attacking  i  targets  is 
Pi  =  i/N. 

By  conditioning  on  the  number  of  weapons  released  by  the  bomber, 
the  probability  that  the  given  weapon  is  released,  denoted  P  ,  is 
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p 

r 


P(releasei  weapons  released)P(i  weapons  released) 


N 

=  Z  (i/N)P(n=i) 
i=o 

N 

=  (1/M)  z  i  P(n=i ) 
l  =o 


=  (l/N)E(n) 


5.3 


where  E(n)  is  the  expected  number  of  weapons  released  by  the  bomber. 

The  number  of  bombers  that  survive  the  BSAM  defense  is  known  (Bp) 
and  hence  the  total  number  of  weapons  that  may  be  released  is  known 
(NBp)  (since  the  bombers  are  assumed  to  have  equal  weapon  loads). 

Assume  that  each  bomber's  set  of  targets  is  located  at  the  same  dis¬ 
tances  xi  along  its  penetration  path.  Then  the  probability  of  being 
released  over  a  target  as  given  by  equation  6.3  is  the  same  for  all 
weapons.  If  it  is  further  assumed  that  each  weapon's  release  event  is 
independent  of  what  has  happened  to  the  other  weapons  (discussed  later), 
then  the  distribution  of  the  total  number  of  weapons  released  is 
binomial : 

BN  u  MB  — K 

Pr(klBp)  =  (  KP)  Pf  1-Fr)  P  o*k<BpN  6.9 

where  B  is  the  number  of  penetrating  bombers  that  have  survived  the 
BSAMs  each  carrying  N  weapons. 

The  distribution  in  Equation  6.9  just  gives  a  distribution  for  the 
number  of  weapons  released.  The  distribution  of  the  number  of  targets 
actually  destroyed  will  not  be  the  same  unless  each  weapon  that  is 
released  has  a  probability  1  of  hitting  and  destroying  its  target. 

It  is  possible  for  the  weapon  to  malfunction,  or,  if  the  target  is 
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extremely  hard,  the  weapon  may  not  impact  close  enough  to  th<=  target  to 
destroy  it.  If  q  is  the  probability  that  a  weapon  that  is  released 
will  destroy  its  target,  then  the  probability  that  a  weapon  is  released 
and  destroys  its  target,  P^,  is 


P 


d 


6.10 


and  the  distribution  of  the  number  of  targets  destroyed  is  given  by 


MB  ,  MB  -k 

vy  ■  <  yyi  -  pd!  p 


6.11 


Two  of  the  assumptions  that  have  been  made  require  further  dis¬ 
cussion.  The  assumption  that  each  weapon's  probability  of  being 
released  is  independent  of  what  has  happened  t  all  the  other  weapons 
is  not  strictly  accurate.  If  it  is  known  that  a  specific  bomber  has 
already  released  half  its  weapons,  then  one  would  expect  to  find  an 
increased  probabiltiy  of  being  released  for  any  weapon  still  aboard. 
However,  it  may  be  reasonable  to  assume  that  the  probability  of  being 
released  for  any  weapon  on  one  bomber  would  be  independent  of  the 
probability  of  being  released  for  any  weaon  on  any  other  bomber.  In 
the  case  where  each  bomber  carries  only  one  weapon,  all  the  weapons 
will  be  independent  and  the  binomial  distribution  given  above  will 
accurately  describe  the  total  number  of  weapons  released.  When  there 
is  more  than  one  weapon  on  each  bomber,  then  for  reasonably  large 
numbers  of  bombers  most  of  the  weapons  will  still  be  independent  of 
each  other  and  the  binomial  distribution  given  above  should  be  a  good 
approximation  to  the  total  number  of  weapons  released. 

The  other  assumption  that  needs  clarification  is  the  assumption 
that  each  bomber's  set  of  targets  is  located  at  the  same  depth 
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along  the  bomber's  penetration  path.  Strictly  interpreted,  this 
assumption  implies  that  all  the  targets  in  the  defended  area  are  located 
on  N  distinct  lines  of  latitude  and  that  there  are  the  same  number  of 
targets  on  each  line.  In  reality,  of  course,  the  targets  will  be  dis¬ 
tributed  over  all  depths  of  penetration.  The  problem  is  to  choose 
values  of  x.  so  that  the  assumed  target  distribution  is  a  reasonably 
close  model  of  the  real  one.  A  method  to  do  this  is  outlined  below. 

One  possible  target  distribution  is  shown  in  Figure  6.1.  Targets 
are  denoted  by  T's.  The  targets  are  not  distributed  uniformly;  more 
targets  are  shown  in  the  middle  parts  of  the  defended  area  than  in  the 
northern  or  southern  extremes.  This  seems  reasonable  but  is  not  a 
necessary  assumption;  any  distribution  would  work. 

The  target  map  in  Figure  6.1  shows  forty-eight  targets.  Assume 
those  targets  are  to  be  attacked  by  a  force  of  eight  bombers  each 
carrying  six  weapons.  The  map  can  be  partitioned  into  horizontal 
bands  each  containing  eight  targets.  This  is  done  in  Figure  6.2.  Now 
each  bomber  can  be  assumed  to  release  one  of  its  weapons  over  one 
target  in  each  of  the  six  bands.  The  only  thing  that  needs  to  be  done 
is  to  choose  a  value  of  x^  for  each  band  that  will  approximate  the 
depth  of  all  eight  targets  within  the  band.  One  choice  would  be  to 
assume  that  all  targets  within  a  band  are  on  the  midline  of  the  band. 

A  more  offense  conservative  choice  is  to  assume  that  all  targets  in 
a  band  are  on  the  southern  boundary  of  the  band. 

An  alternative  method  is  available  if  it  is  assumed  that  the 
analyst  knows  the  distribution  of  targets  as  a  function  of  penetration 
depth.  Consider  a  pdf  that  gives  a  random  target's  probability  density 


FIGURE  6.1 

Distribution  of  targets  within  defended  area. 


FIGURE  6.2 

Target  set  partitioned  Into  six  bands. 
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of  being  located  at  a  depth  x.  This  probability  might  follow  a  normal, 
uniform,  or  any  other  distribution  one  wished  to  assume.  A  possible 
pdf  is  given  in  Figure  6.3. 


FIGURE  6.3 

Probability  distribution  of  targets  as  a  function  of  depth. 

This  pdf  can  be  partitioned  into  six  intervals  of  equal  area, 
corresponding  to  bands  with  equal  probabilities,  1/6,  that  the  target 
would  be  within  each  band.  If  the  intervals  are  (C,a),  (a,b),  ••*, 
(e,D),  then  the  values  of  a,  ...,  e  will  be  those  values  such  that 

/*ptx)dx  *  /gP(x)dx  =  ...  =  /gP(x)dx 

The  resulting  partition  Is  sketched  in  Figure  6.4.  Since  each 
of  the  intervals  defines  a  band  within  which  a  target  has  equal 
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FIGURE  6.4 

Target  pdf  partitioned  into  eight  equal  parts. 

probability  of  being  located,  then  for  large  numbers  of  targets  that 
follow  the  assumed  distribution,  one  could  expect  to  find  equal 
numbers  of  targets  within  each  band.  As  before,  it  still  remains  to 
choose  a  value  of  x.  for  each  band.  Various  procedures  exist  that 
could  be  used  to  generate  values  that  are  mathematically  best  under 
various  criteria;  these  will  not  be  discussed  here.  For  many  applica¬ 
tions,  either  of  the  simple  choices  mentioned  previously  would  be 
adequate. 

In  this  computer  model,  the  x^  are  assumed  to  be  Input,  the’model 
calling  for  as  many  x^  as  there  are  weapons  on  each  bomber. 
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It  was  pointed  out  earlier  in  this  section  that  Py  need  not  be 
the  same  for  bombers  flying  at  high  and  low  altitudes.  The  same 
applies  to  P  (x).  Therefore,  there  will  be  two  distributions  for  the 
number  of  weapons  that  destroy  targets,  corresponding  to  each  of  the 
two  groups  of  bombers.  These  distributions  can  be  combined  to  give 
the  distribution  for  the  total  number  of  weapons  delivered  (see 
Chapter  VII). 

RECOVERY 

The  only  aspect  of  the  recovery  problem  to  be  addressed  by  this 
model  is  the  bomber's  problem  of  leaving  the  defended  area  after 
releasing  all  its  weapons. 

In  this  model,  it  is  being  assumed  that  after  attacking  its  last 
assigned  target,  a  bomber  will  continue  to  penetrate  the  rest  of  the 
way  through  the  corridor  and  exit  at  the  bottom  or  southern  end. 

Using  the  penetration  probability  function  Pp(x),  the  bomber's  chance 
of  surviving  the  area  defenses  to  reach  the  end  of  the  defended  area 
is  PS=PS(D)  where  D  is  the  total  depth  of  the  defended  area. 

It  is  assumed  that  if  the  bomber  has  reached  the  end  of  the 
defended  area,  he  has  already  attacked  all  his  assigned  targets  and 
survived  the  terminal  defenses  of  each.  The  probability  that  a 
bomber  that  has  survived  the  BSAM  defense  to  begin  penetration  will 
survive  to  exit  the  defended  area,  Pg,  is,  therefore, 

Pe  =  Pp(0)p!J  6.13 

where  N  is  the  number  of  targets  he  attacked  during  penetration. 
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Assume  that  each  bomber's  recovery  event  is  independent  of  every 

other  bomber's  recovery  event.  Then  the  distribution  of  the  number  of 

bombers  that  exit  the  defended  area  for  a  given  number  B  of  bombers 

P 

that  survive  the  BSAM  defense  is  binomial. 


pR(k;6p5 


=  (  p 
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P^(i  - 


B  -k 


6.14 


Since  both  P  (x)  and  Py  may  have  different  values  for  bombers 
flying  at  high  and  low  altitudes,  there  will  be  two  recovery  distri¬ 
butions,  one  for  each  force  of  bombers. 


CRUISE  MISSILES 

This  section  will  develop  the  representation  for  cruise  missiles 
(CM)  under  each  of  the  three  cases  discussed  in  Chapter  III.  Specifi¬ 
cally,  these  cases  are  short  range  cruise  missiles  (SRCM)  that  are 
released  at  various  points  during  penetration  after  the  bombers  have 
passed  through  the  BSAM  band,  intermediate  range  cruise  missiles  (IRCM) 
which  are  released  after  the  FAD  zone  but  before  the  bombers  encounter 
the  BSAM  defense,  and  long  range  cruise  missiles  (LRCM)  which  can  be 
released  before  the  bombers  enter  the  FAD  zone. 

Before  discussing  each  of  the  three  cases,  some  general  assumptions 
will  be  made  about  the  effect  CMs  would  have  on  a  bomber's  penetration 
probabilities  that  have  been  discussed  in  the  last  three  chapters. 

The  first  assumption  to  be  made  is  that  introduction  of  CMs  will 
not  effect  a  bomber's  probability  of  penetrating  through  BSAM  or  RASAM 
defenses.  This  assumption  is  reasonable  in  cases  where  there  are 
Insufficient  penetrators  to  saturate  the  SAM  defenses.  As  the  number 
of  penetrators  becomes  larger,  it  would  be  possible  to  degrade  the 
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defenses  simply  by  exhausting  the  supply  of  missiles  at  the  SAM  sites. 
For  large  total  numbers  of  penetrators,  an  individual  penetrator's 
chance  of  surviving  the  SAM  defenses  would  tend  to  increase  and  the 
assumption  would  become  quite  offensive  conservative. 

The  next  assumptions  deal  with  the  effect  of  CMs  on  a  bomber's 
probability  of  surviving  the  area  fighter-interceptor  defenses.  Recall 
that  this  probability  depends  on  a  parameter  A  which  gives  the  average 
number  of  intercepts  possible  on  a  penetrator  per  unit  time.  It  has 
been  assumed  in  this  model  that  X  is  a  constant.  It  is  now  further 
assumed  that  x  remains  a  constant  when  CMs  are  added  to  the  problem 
and  that  x  is  the  same  for  both  bombers  and  cruise  missiles.  This 
does  not  imply  that  X  must  be  the  same  constant  for  a  model  with  CMs 
added  as  it  was  for  a  model  without  CMs,  only  that  it  is  constant  within 
each  model  and  that  it  is  the  same  for  both  CMs  and  bombers. 

The  last  area  in  which  CMs  could  effect  a  bomber's  probability  of 
survival  is  the  FAD.  If  LRCMs  have  been  released  and  are  entering  the 
FAD  zone  with  the  bombers,  it  is  assumed  that  the  defense  will  try  to 
attack  them  as  well  as  the  bombers.  Since  there  are  a  fixed  number  of 
defensive  interceptors,  each  capable  of  attacking  only  one  penetrator, 
an  increase  in  the  number  of  penetrators  will  decrease  the  number  of 
interceptors  assigned  to  each  penetrator.  This  will  increase  each 
individual  penetrator's  probability  of  surviving  the  FAD.  This  is 
accounted  for  in  this  model. 

Short  Range  Cruise  Missiles.  Consider  a  cruise  missile  of  range  R 
destined  for  a  target  at  some  distance  x+R  inside  the  defended  area. 
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bomber  will  release  the  CM  at  its  maximum  range,  the  CM  will  -be 
released  when  the  bomber  has  penetrated  to  a  dpeth  x  inside  the 
defended  area.  In  order  for  the  CM  to  be  released,  the  bomber  carrying 
it  must  survive  to  at  least  the  depth  x.  The  probability  of  this  is 
given  by  the  penetration  probability  function  derived  in  the  last 
chapter . 
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Recall  that  the  term  inside  the  second  set  of  square  brackets  in 
Equation  6.15  is  Pp(x),  the  probability  that  a  bomber  survived  to  at 
least  a  depth  x  through  FI  defenses.  The  last  step  of  the  derivation 
of  Pp(x)  in  Chapter  V  was  the  integration  of  Equation  5.23. 
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If  the  terms  in  equation  6.16  are  combined,  one  will  find  that 
the  expression  is  equivalent  to  the  final  form  of  Pp(x)  given  in 
equation  5.24  and  used  in  equation  6.15.  However,  in  the  expression 
for  PF(x)  in  equation  6.15,  each  term  has  an  interpretation  in  terms 
of  the  model.  The  first  set  of  square  brackets  in  equation  6.15  gives 
the  probability  that  the  bomber  is  detected  by  the  defense  but  still 
survives  to  at  least  a  depth  x.  The  term  in  the  second  set  of  square 
backets  is  the  probability  that  the  bomber  has  reached  a  depth  x 
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undetected.  (See  equation  5.6  and  the  accompanying  text).  These 
interpretations  will  be  useful  in  the  development  of  the  SRCM  model 
and  for  this  reason  equation  6.14  is  rewritten 
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where  equation  6.15  and  the  distributive  property  have  been  used. 

Return  now  to  the  cruise  missile.  If  it  is  successfully  released 
at  a  depth  x  inside  the  defended  area,  it  must  then  penetrate  a  further 
distance  R  on  its  own.  Since  it  is  assumed  that  the  CM  is  threatened 
by  the  same  defenses  that  threaten  the  bombers,  it  is  reasonable  to  use 
the  same  general  analysis  that  was  developed  in  Chapter  V  for  bomber 
penetration  through  area  defenses.  The  only  modifications  will  be  to 
account  for  the  fact  that  at  the  time  of  its  release,  the  CM  is  already 
some  distance  into  the  defended  area. 


As  has  already  been  suggested  by  equation  6.16,  there  a» a  two 
possible  situations  at  the  time  the  bomber  releases  the  CM.  First, 
the  bomber  may  have  penetrated  to  the  release  point  undetected.  The 
probability  of  this  event,  call  it  Px(u)  has  already  been  shown  to  be 
the  second  term  on  the  RHS  of  Equation  6.17. 


Px(u) 
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It  seems  reasonable  to  believe  that  if  the  bomber  has  not  been 
detected  at  the  time  of  release  of  the  CM,  then  the  CM  itself  will 
begin  its  flight  undetected.  In  this  case,  the  analysis  that  the  CM 
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will  penetrate  an  additional  depth  R  on  its  own  will  parallel  exactly 
the  analysis  given  in  Chapter  V  for  bombers. 

The  analysis  in  Chapter  V  of  bomber  penetration  through  RASAMs 
gave  a  probability  of  bomber  survival  P$U)  that  depended  on  the  depth 
of  penetration  x  into  the  defended  area.  The  analysis  assumed  that  x 
was  measured  from  the  northern  boundary  of  the  defended  area  since  there 
was  no  way  a  bomber  could  begin  penetration  at  any  point  south  of  the 
boundary.  However,  a  cruise  missile,  because  it  is  being  released 
from  a  bomber  which  has  already  penetrated  to  a  depth  x,  begins  its 
own  flight  at  a  distance  x  inside  the  boundary.  Nonetheless,  using 
the  same  assumptions  and  analysis,  it  can  be  shown  that  the  CMs 
probability  of  penetrating  a  distance  R  from  its  release  point  x  is  a 
function  only  of  the  distance  R,  not  x.  This  is  because  the  analysis 
used  a  ratio  of  two  areas  to  find  the  probability  of  encountering  a 
particular  SAM  site.  The  area  in  the  denominator  was  just  the  area  of 
the  entire  defended  area.  The  area  in  the  numerator  depended  on  the 
length  of  the  penetration  path  and  the  PLD  of  a  SAM  site;  this  would 
be  the  same  for  a  path  of  given  length  no  matter  where  the  path  began. 
See  Figure  6.5. 

The  probability  that  a  cruise  missile  penetrates  a  distance  R 
through  RASAM  defenses  is  therefore 

pscM  ■  n-'kscTff  )Nr  6-19 

Now  consider  the  FI  defenses.  Since  the  parameter  x  is  assumed 
to  be  a  constant,  the  number  of  intercepts  that  are  possible  against 
the  cruise  missile  depends  only  on  the  amount  of  time  from  first 
detection  of  the  CM  until  it  reaches  its  target.  The  time  until 
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FIGURE  6 ;5 

Cruise  missile  vs.  random  area  SAMs. 


detection  is  distributed  exponentially  and  therefore  will  not  depend 
on  where  the  CM  begins  its  flight.  The  CMs  probability  of  penetrating 
a  distance  R  through  FI  defenses  will  not  depend  on  where  it  begins 
penetration  and  can  be  found  using  the  analysis  for  FI  defenses  in 
Chapter  V. 
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Combining  equations  6.19  and  6.20,  one  finds  the  CM's  probability 
of  penetrating  a  distance  R  through  the  area  defenses,  given  that  it 
is  undetected  at  the  time  it  is  released,  to  be 
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The  other  possible  situation  at  the  time  the  bomber  releases  the 
CM  is  that  the  bomber  has  already  been  detected  but  has  still  survived 
to  reach  the  release  point.  The  probability  of  this  event,  P  (V),  is 
given  by  the  first  term  on  the  RHS  of  Equation  6.15 
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If  the  bomber  is  being  tracked  at  the  time  it  releases  the  CM, 
it  seems  reasonable  to  assume  that  the  CM  will  also  be  detected  at  the 
time  of  its  release.  This  will  not  affect  the  CM’s  probability  of 
surviving  to  travel  a  distance  R  through  RASAM  defenses.  This  will 
remain  the  same  as  given  in  equation  6.19.  However,  if  the  CM  is 
detected  as  it  is  released,  then  the  defense  can  immediately  begin 
assigning  interceptors  to  it.  The  total  number  of  intercepts  that  can 
occur  in  the  time  it  takes  the  CM  to  travel  a  distance  R  will  be,  by 
the  assumptions  made  in  Chapter  V  about  FI  defenses,  a  Poisson  random 
variable  with  parameter  A.  From  equation  5.16,  P^R),  the  probability 
of  i  intercepts  in  a  distance  R  is 

i  -xR/vc 

P^R)  -  (AR/vc)’e  c/i!  6.23 

To  travel  the  total  distance  R  through  the  FI  defenses,  the  CM 
must  survive  all  such  attacks,  each  of  which  may  kill  the  CM  with 
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Combining  equations  6.19  and  6.24  gives  the  probability,  PR  (R!V) 
that  a  CM  penetrates  a  distance  R  on  its  own  through  the  area  defenses, 
given  that  it  is  detected  as  soon  as  it  is  released,  to  be 


w.R  Nr  -Pkf  *R/v 


Using  the  above  analysis,  one  can  now  find  the  probability  that  a 
CM  is  released  at  a  depth  x  and  penetrates  on  its  own  an  additional 
distance  R,  thus  reaching  its  target.  By  conditioning  this  probability, 
denoted  P^Cx),  on  whether  the  bomber  has  reached  the  CM  release  point 
x,  detected  (V)  or  undetected  (U),  one  finds 


PCM(x)  *  pPc(R!V)Px(V)  +  PPc(R[U)Px(U) 
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Equations  6.19,  6.20,  6.22,  and  6.25  can  be  substituted  for  the 
quantities  on  the  RHS  of  Equation  6.26,  giving,  after  some 
rearranging. 
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Equation  6.27  gives  the  probability  that  a  SRCM  reaches  its 
target  through  the  area  defenses  as  a  function  of  its  release  point, 
x.  When  the  CM  reaches  its  target,  it  win  encounter  one  last  defense, 
namely  a  ISAM.  The  probability  that  the  CM  survives  the  TSAM  will  be 
denoted  by  Py  . 

Equation  6.27  also  fails  to  account  for  the  possibility  that  a  CM 
doesn't  reach  its  target  because  of  some  kind  of  malfunction.  Another 
possibility  is  terrain  clobber.  Since  CMs  may  fly  at  very  low  alti¬ 
tudes,  there  is  a  chance  that  they  would  run  into  a  natural  or  manmade 
object  that  rose  sharply  above  the  surrounding  terrain.  These  non¬ 
defense  related  obstacles  to  the  CM  will  be  represented  by  qc>  where 
qc  is  the  probability  that  a  CM  fails  to  reach  and  destroy  its  target 
for  reasons  not  related  to  the  defenses. 


All  this  can  now  be  put  together  to  give  the  probability  that  a 
SRCM  with  a  target  at  x+R  will  be  released  by  its  bomber  and  will 
reach  and  destroy  its  target.  This  probability  is  a  function  of  x  and 
is  denoted  PCM(x). 


PCM(x)  s  pCM^PTc^"qc^ 


6.28 


To  find  an  average  probability  that  a  CM  reaches  and  destroys  its 
target,  it  is  necessary  to  know  something  about  the  distribution  of 
the  targets  within  the  defended  area.  Consider  a  pdf  for  the  location 
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(in  terms  of  penetration  depth)  of  a  random  target,  p  (y)  .  Since 
these  targets  are  all  assumed  to  be  at  least  a  distance  R  inside  the 
defended  area,  the  pdf  might  look  something  like  the  one  shown  in 
Figure  6.6. 


FIGURE  6.6 

Density  function  for  CM  targets  as  a  function  of  depth. 

The  probability  that  a  CM  reaches  and  destroys  its  target  given 
in  equation  6.28  is  a  function  of  its  release  point  x,  not  the  depth 
of  the  target  y.  However,  note  that 

x  =  y  -  R.  6.29 
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This  is  because  it  is  assumed  that  the  CM  will  always  be  released  at 
its  maximum  range  R  away  from  its  target. 

The  average  probability  that  a  CM  reaches  and  destroys  its  target, 
denoted  P^,  can  now  be  found  to  be  conditioning  on  the  location  of  its 
target 

PCM  =  ’"r  PCM^y‘R'pc^y'dy  ‘  6-30 

In  practice,  this  integral  would  be  difficult  to  impossible  to 
evaluate  analytically.  However,  a  good  approximation  could  always  be 
found  using  any  of  several  different  computer  integration  techniques. 

Once  p£w  has  been  found,  the  distribution  of  the  number  of  targets 
destroyed  will  be  binomial 

V"!  ■  <Bk>PCH<1-PCM>BC'k  6-31 

where  B  is  the  number  of  bombers  that  have  survived  the  BSAM  defense 
and  C  is  the  number  of  CMs  on  each  bomber.  As  in  the  case  of  weapons 
delivered  discussed  in  the  previous  section,  the  use  of  the  binomial 
distribution  implies  the  assumption  that  each  CM's  probability  of 
success  is  independent  of  what  has  happened  to  all  the  others.  The 
comments  made  on  this  assumption  in  the  section  on  weapons  delivered 
are  equally  appropriate  here. 

Intermediate  Range  Cruise  Missiles.  Consider  next  a  CM  that  is 
released  after  the  FAD  zone  but  before  the  bomber  reaches  the  BSAM 
band.  The  CM  is  destined  for  a  target  at  some  distance  y  inside  the 
defended  area. 
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To  reach  its  target,  the  CM  must  first  survive  the  BSAM  defense 
and  then  penetrate  on  its  own  to  a  depth  y  inside  the  defended  area. 
Each  part  of  the  CM's  problem  can  be  treated  exactly  the  same  as  for 
a  bomber. 

In  Chapter  IV,  the  probability  that  a  particular  bomber  survived 
the  3SAM  defense  was  derived  (Equation  4.22).  For  a  cruise  missile, 
the  analysis  procedes  in  the  same  way  and  one  finds 
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6.32 


If  the  CM  survives  the  BSAM  defense,  then,  like  a  bomber,  it  will 
begin  penetration  undetected.  Its  probability  of  penetrating  to  a 
depth  y  within  the  defended  area  is  given  by  the  penetration  probability 
function  Ppc(y) 
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The  CM's  total  probability  of  reaching  its  target  through  the 
defenses,  P^(y),  can  now  be  expressed  as  the  product  of  its  probability 
of  surviving  the  BSAM  defense  and  its  probability  of  penetrating  the 
area  defenses  to  a  depth  y  within  the  country. 

PCM^  =  PBScPPc^ 
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6.34 
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As  already  described,  the  CM  has  a  probability  of  failing  to  reach 
or  destroy  its  target  that  is  unrelated  to  the  defense,  q  .  And  if  it 
reaches  its  target,  it  still  faces  the  TSAM  defense  (°jc)-  These 
factors  combine  as  in  equation  6.29  to  give 

PCM(y)  »  PcM(y)  PTc  (1'qc}  6-35 

where  P^(y)  1s  the  probability  a  CM  reaches  its  target  if  the  target 
is  at  a  depth  y. 

If  the  distribution  of  the  CM's  targets  is  known,  the  average 
probability  that  an  IRCM  reaches  and  destroys  its  target  can  be  found 
by  conditioning  on  the  target's  location. 

PCM  =  ;0  PCM(y)pc(y)dy  6-36 

Notice  that  since  the  CMs  are  released  outside  the  defended  area  any¬ 
way,  it  is  no  longer  assumed  that  the  targets  are  at  a  distance  greater 
than  the  range  of  the  CM  inside  the  defended  area. 

Given  PCM,  the  distribution  of  the  number  of  CMs  that  reach  and 
destroy  targets  is  given  by  equation  6.31  where  is  the  total  number 
of  bombers  that  survive  the  FAD  and  C  is  the  number  of  CMs  on  each 
bomber. 

BcC  .  BcC-k 

’dcWM  =  <•  [  > <'-'W  6-37 
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Long  Range  Cruise  Missiles.  Finally,  consider  the  case  of  a 
cruise  missile  of  sufficiently  long  range  that  it  can  be  released 
before  the  bomber  reaches  the  FAD  zone.  The  CM  is  destined  for  a  tar¬ 
get  at  depth  y  inside  the  defended  area. 

To  reach  its  target,  the  CM  must  now  survive  both  the  FAD  zone 
and  the  5SAM  and  then  penetrate  to  its  target  at  a  distance  y  inside 
the  defended  area.  The  analysis  for  the  BSAM  defense  and  for  pene¬ 
tration  against  the  area  defenses  will  be  the  same  as  for  IRCMs.  The 
only  thing  that  has  to  be  added  is  the  probability  that  the  CM  survives 
the  FAD. 

As  in  the  case  of  both  the  BSAM  and  area  defenses,  the  analysis 
of  a  CM's  probability  of  surviving  the  FAD  will  exactly  parallel  the 
analysis  done  for  bombers.  For  a  given  number  nc  of  fighters  assigned 
to  intercept  the  CM,  the  probability  that  the  CM  survives  all  of  the 
attacks  and  so  survives  the  FAD  is  given  by  Equation  4.17, 

Pc  ■  (1  -  Pkc)"C  6.38 

where  Pkc  is  the  probability  that  the  CM  is  killed  by  a  single  FAD 
intercepter  assigned  assigned  to  it. 

It  remains  to  determine  the  number  of  fighters,  nc>  assigned  to 
attack  the  CM.  In  the  section  in  Chapter  IV  on  bombers  against  the 
FAD,  it  was  assumed  that  the  defense  would  try  to  assign  fighters 
uniformly  to  all  approaching  bombers.  In  the  present  case  where  LRCMs 
are  added  to  the  attack,  the  assumption  will  be  broadened  to  assign 
the  fighters  uniformly  to  all  approaching  penetrators,  including  both 
bombers  and  CMs .  This  assumption  might  imply  that  the  defense  feels 
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that  individual  bombers  and  CMs  pose  a  roughly  equal  threat  or  simply 
that  the  defense  cannot  distinguish  between  a  bomber  and  a  CM  on 
radar  and  so  assigns  fighters  uniformly  to  all  penetrators. 

If  there  are  SgC  bombers  and  C  cruise  missiles  entering  a  FAD 
zone  defended  by  F  fighters,  then  each  penetrator  will  have  at  least 
nc  fighters  assigned  to  it  where 


nc  •  <  F/Sbe  +  C)>  6.39 

If  F/(BgE+C)  is  a  noninteger,  the  model  assigns  the  remaining 
fighters  one  each  to  the  bombers  penetrating  at  high  altitude,  consistent 
with  the  analysis  in  Chapter  IV.  If  there  are  still  fighters  remaining, 
it  is  arbitrarily  assumed  that  they  are  assigned  next  to  the  low  alti¬ 
tude  bombers  and  finally  to  the  CMs. 

Equations  6.34  and  6.38  can  now  be  combined  to  give  the  probability 
that  a  CM  survives  the  FAO  and  8SAM  defenses  and  penetrates  the  area 
defenses  to  its  target  at  a  depth  y  inside  the  defended  area. 
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Once  again,  account  must  be  made  of  the  CM's  probability  of  sur¬ 
viving  the  TSAM  defense  and  of  the  probability  that  it  fails  to  reach 
or  destroy  its  target  for  reasons  unrelated  to  the  defenses.  This  gives 

PCM^  *  pCM^PTc^‘qc^  6,41 
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As  before,  this  equation  can  be  multiplied  by  the  pdf  of  the 


location  of  the  target  and  integrated  to  give 


CM  =  ;0  PCM^Pc^dy 


6.42 


Again,  given  P^,  the  distribution  of  the  number  of  targets  de¬ 
stroyed  is  binomial  and  is  given  in  Equation  6.30  or  6.36.  is 

the  number  of  bombers  that  survive  the  base  escape. 
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High  vs.  Low  Altitude.  As  has  been  the  case  in  all  parts  of  this 
model  since  aerial  refueling,  it  is  necessary  to  consider  the  effect 
due  to  two  groups  of  bombers  flying  at  different  altitudes. 

It  is  assumed  that  all  CMs,  regardless  of  what  altitude  they  are 
released  at,  will  penetrate  at  some  low  altitude.  With  this  assumption, 
the  altitudes  of  the  releasing  bombers  will  not  effect  the  distribution 
of  the  number  of  LRCMs  that  reach  their  targets.  All  the  bombers  that 
survive  the  base  escape  will  release  their  CMs  whether  they  are  refueled 
or  not  and  the  CMs  will  all  penetrate  at  the  same  altitude;  thus  they 
will  all  have  the  same  probability  of  reaching  and  destroying  their 
targets. 

In  the  case  of  IRCMs,  all  the  CMs,  regardless  of  the  altitude  at 
which  they  are  released,  will  go  to  low  altitude  and  thus  all  have 
equal  probabilities  of  reaching  their  targets.  However,  the  total 
number  of  CMs  released  will  depend  on  the  total  number  of  bombers  that 
survive  the  FAD  which  will  depend  on  the  numbers  of  bombers  penetrating 
at  high  and  low  altitudes. 
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In  the  case  of  SRCMs,  each  CM  will  have  the  same  probability  of 
reaching  its  target  given  that  it  is  released,  since  they  will  all 
fly  at  the  same  altitudes.  However,  a  CM's  probability  of  being 
released  depends  on  its  bomber's  penetration  probability  function  which 
in  turn  depends  on  the  altitude  at  which  the  bomber  is  penetrating. 
Therefore,  in  this  case,  two  separate  distributions  are  used,  giving 
the  number  of  targets  destroyed  by  CMs  released  by  the  two  different 
groups  of  bombers.  These  distributions  can  be  combined  to  give  the 
total  number  of  targets  destroyed  in  the  case  of  SRCMs. 
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CHAPTER  VII 


INTEGRATION  AND  COMPUTERIZATION 

The  preceding  chapters  have  described  the  model  representations 
for  each  of  the  stages  of  the  strategic  bomber  mission.  This  chapter 
will  briefly  describe  how  the  various  stages  of  the  model  are  integrated 
for  computerization. 

In  the  program  written  for  this  paper,  there  are  six  basic  outputs 
desired.  These  are  a)  a  mean  and  a  variance  for  the  total  number  of 
weapons  delivered,  b)  a  mean  and  a  variance  for  the  total  number  of 
bombers  that  survive  the  mission,  and  c)  a  mean  and  a  variance  for  the 
total  number  of  cruise  missiles  that  reach  and  destroy  their  targets. 
These  are  not  the  only  possible  outputs  that  might  be  obtained  from  the 
model;  with  relatively  minor  modifications,  the  program  could  give  inter¬ 
mediate  results  such  as  the  number  of  bombers  that  survive  any  parti¬ 
cular  layer  of  defense,  for  example  BSAMs.  The  program  could  also  be 
expanded  to  give  plots  of  the  various  distributions  instead  of  just 
means  and  variances. 

Conditioning 

To  find  the  mean  and  variance  of  the  quantities  listed  above,  it 
is  generally  necessary  to  know  their  distributions.  Since  the  number 
of  weapons,  bombers,  and  cruise  missiles  will  always  be  integers,  the 
distributors  will  all  be  discrete.  Given  a  discrete  distribution 
P(K* k),  the  mean,  E{K),  and  variance,  Var(K),  are  expressed  by  the 
following  formulas: 

E ( K)  *  Z  1P(K-1)  7.1 

i 
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7.2 


Var(K)  »  E  i2P(K=i)  -  E(K)2 
i 

Thus,  to  obtain  the  six  desired  outputs,  it  is  necessary  to  find 
three  distributions,  one  for  the  total  number  of  weapons  delivered, 
another  for  the  total  number  of  bombers  that  survive  the  mission,  and 
the  third  for  the  total  number  of  cruise  missiles  that  reach  and 
destroy  their  targets. 

Obtaining  any  one  of  these  distributions  is  a  problem  in  condi¬ 
tional  probability.  The  reasoning  process  for  all  three  distributions 
will  be  the  same.  Because  it  is  convenient  to  deal  with  a  specific 
example,  the  distribution  for  the  total  number  of  bombers  that  survive 
the  mission  will  be  used. 

It  is  desired  to  find  the  probability  that  K  bombers  survive  the 
mission  where  the  initial  number  of  bombers  that  begin  the  mission  is 
a  input  to  the  problem.  To  find  this  probability,  work  backwards  from 
the  last  stage  of  the  mission,  recovery,  to  the  first,  base  escape. 

Recall  that  the  model  is  dealing  with  two  groups  of  bombers,  one 
flying  at  high  altitude  and  the  other  flying  at  low  altitude.  Let  K 
be  the  total  number  of  bombers  that  survive  and  let  KH  and  be  the 
numbers  that  survive  at  high  and  low  altitudes  respectively.  Then  one 
has 


P(K*k)  =  P(KH+KL=k)  7.3 

By  conditioning  on  the  number  of  bombers  that  survive  at  high 
altitude,  equation  7.3  can  be  expressed  as 
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7.4 


P( K=k)  =  EP(KH+KL-k; KH*h)P(KH-h) 
h 

=  :P(  K.  =k-h)P(KM=h) 
h  L  H 

Now  let  5p^  and  Bp,  represent  the  number  of  bombers  that  have 
survived  the  BSAMs  at  each  altitude.  Then  the  two  probabilities 
?('<^=k-h)  and  P ( K,_j = h )  can  be  conditioned  on  Bp?  and  3p^  respectively 
to  give 

P(KH=h)  =  :P(KH=hj8ph=i)P(Bph=i)  7.5a 

P(Kl=0  =  rP(KL-a;Bpi-i)P(Bpjl»i)  7.5b 

where  for  convenience  k-h  has  been  replaced  by  l. 

An  expression  for  the  quantities  P( -H=h | Bph=i)  and  P( j Bpi=i ) 
was  derived  in  Chapter  VI  in  the  section  on  recovery  and  is  given  in 
equation  6.14. 

The  terms  P(Bp^=i )  and  P(Bpi=i)  can  each  be  expanded  by  conditioning 
on  the  numbers  of  bombers  that  survive  the  FAD  at  high  and  low  altitudes 
respectively.  Denote  these  numbers  Bph  and  Bp^.  Then 

P(Bph=i)  -  zP(Bph-i|Bph=j)P(Bph«j)  7.6a 

J 

P(Bpt  =  iJ  ■£PlBM-1|BFt-J)P(BFt-J)  7.6b 

j 

Expressions  for  p(Bph=i i Bph= j )  and  P(Bpi=i |Bpt=j)  were  derived  in 
Chapter  IV  in  the  section  on  BSAM  defense.  They  are  given  in  Equations 
4.35a  and  b. 

The  terms  P(Bp^=j)  and  P(Bpj=j)  can  each  be  conditioned  on  the 
numbers  of  bombers  entering  the  FAD  zone  at  each  altitude.  If  these 
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numbers  are  denoted  and  8,,  then 


P(5Fti’J'  ‘ 

T  P(B  =j;Bh=m)P(Sh=m) 
m 

7.7a 

r  P(BFz=j!8.=n)P(B,=n) 

7.7b 

n 


In  Chapter  IV,  an  expression  was  found  for  the  terms  P(Br.  =j  8,  =rrl 

Fh  n 

and  P(Bp. =j ;B. =n) .  This  is  given  oy  Equation  4.23. 

The  numbers  of  bombers  flying  at  high  and  low  altitudes  is  deter¬ 
mined  by  the  number  of  bombers  that  successfully  refuel.  Specifically, 
one  has 


P(Bn=m)  »  P(m  bombers  fail  to  refuel)  and  7.3a 

PlB^n)  =  P(n  bombers  successfully  refuel)  .  7.3b 


Each  of  these  probabilities  can  be  conditioned  on  the  number  of 
bombers  that  survive  the  base  escape,  BBE. 


P(Bh=m)  =  SP(Bh=mJBBE=b)P(B8£=b) 
b 


7.9a 


P(B^n)  =  ZP(Bjl=niBBE=b)P(BBE=b) 


7.9b 


Now  consider  that  the  sum  of  the  numbers  of  bombers  flying  at  high 
altitude  and  at  low  altitude  must  be  the  total  mimber  of  bombers  that 
survived  the  base  escape,  B0E. 


B 


h 


+  Bi 


7.10 


Using  7.10,  the  conditional  probability  P(8^=m 
expressed  in  terms  of  as  follows 


BBE=b)  can  be 
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7.11 


P(B^=m; Bg£=b)  =  P(3.,=b-m;  Bgg=b) 

Equation  7.11  can  be  substituted  into  7.9a  to  give  (with  equation 
7.9b  unchanged) 

KV*: '  i  7-12a 

0 

°( 3 . =n )  =  :  P(B,=n'3nr=b)P(2n^=b)  7.12b 

Vi  1_  c  D  fc,  DL 

b 

The  quantity  P (B^  =  i i 8gg=b)  has  been  derived  in  Chapter  IV  in  the 
section  on  refueling.  It  is  given  in  either  Equation  4.18  or  4.21  as 
appropriate.  The  expression  for  P(Bgg=b)  was  derived  on  the  base 
escape  section  of  Chapter  IV  and  is  given  in  equation  4.16  when  the 
initial  numbers  of  bombers  on  the  gound  and  in  the  air  prior  to  the 
attack  are  given. 

This  completes  the  integration  of  the  model  to  find  the  distribu¬ 
tion  of  the  number  of  bombers  that  survive  the  mission.  Parallel 
derivations  can  be  done  to  give  the  distributions  of  the  numbers  of 
w eaons  and  cruise  missiles  delivered  on  targets. 

COMPUTERIZATION 

In  writing  the  computer  program  for  this  paper,  the  intent  was  to 
demonstrate  the  validity  and  practicality  of  the  model  that  has  been 
developed.  The  resulting  computer  model  is  therefore  not  intended  or 
expected  to  be  optimum  in  terms  of  flexibility  or  efficiency.  However, 
it  is  hoped  that  the  program  might  provide  a  starting  point  for  someone 
who  desired  to  write  his  own. 
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For  convenience,  one  assumption  and  one  approximation  have  been 
built  into  the  program.  These  are  described  below. 

The  assumption  involves  the  base  escape  part  of  the  model.  Recall 
that  when  this  part  was  developed  in  Chapter  IV,  the  distribution  for 
the  time  T^  which  the  bombers  have  to  take  off  in  was  left  as  an 
unspecified  input.  For  convenience,  in  the  computer  program  this  dis¬ 
tribution  has  been  built  in  as  a  Beta  distribution  with  parameters 
and  3=5.  The  Beta  distribution  is  given  by 

f(x)  »  ( -  x'Mi-x)5  o<x<i  7.13 

when  a  and  3  are  integers.  In  the  model,  the  maximum  and  minimum  pos¬ 
sible  values  of  T^  are  input  and  the  computer  fits  the  distribution 
between  them  using  the  substitution  x=  (t-T  .  )/(T  -T  .  ). 

For  illustrative  purposes,  this  program  contains  the  code  for 
modelling  long  range  cruise  missiles.  The  approximation  made  in  the 
program  involves  both  the  CMs  and  the  bombers  in  the  FAD  part  of  the 
mission.  Recall  that  because  the  numbers  of  fighters,  bombers  and 
CMs  are  all  integers,  it  is  possible  to  have  some  bombers  with  one 
more  interceptor  assigned  to  them  than  the  other  bombers  or  some  CMs 
with  one  more  inteceptor  assigned  to  them  than  the  other  CMs.  In  other 
words,  within  any  one  of  three  distinct  groups  -  high  altitude  bombers, 
low  altitude  bombers,  or  CMs,  -  it  will  be  possible  to  have  two  sub¬ 
groups  of  penetrators  having  different  probabilities  of  surviving  the 
FAD.  The  interceptor  assignment  policy  which  assigns  remaining  inter¬ 
ceptors  first  to  high  altitude  bombers,  then  to  low  altitude  bombers. 
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and  finally  to  CMs ,  prevents  this  subgrouping  form  from  occurring  in 
more  than  one  of  the  three  groups  of  penetrators. 

The  model  treats  these  subgroups  separately  and  then  convolutes 
the  distributions  of  the  number  of  survivors  from  each  subgroup  to 
find  the  total  number  of  survivors  from  the  particular  group  of  pene¬ 
trators.  This  is  done  in  equations  d.27  and  4.28.  These  equations 
are  repeated  nere  with  the  assumption,  for  the  sake  of  example,  that  it 
is  the  CMs  that  are  divided  into  subgroups. 

PS1  =  (1-Pk)  <F/8+C>  +  1  7.14a 

ps2  =  (1*pk)  <p/B+c>  7.14  b 

max( k  ,R)  _  _  .  .  „  _  ,  . 

P(k  survive  FAD)  =  I  (?  ?)Pk  1  (l-p_?)C'R  k+1 

i-max(o,k-C+R)  k 

'(i)Psl{1-psl)R_1'  7AS 

where  F  is  the  total  number  of  fighters,  B  is  the  total  number  of 
bombers,  c  is  the  total  number  of  CMs,  and  R  is  the  number  of  CMs  with 
one  extra  fighter  assigned  to  them. 

Since  the  FAD  routine  is  called  numerous  times  by  the  main  computer 
program,  time  can  be  saved  by  eliminating  the  summation  in  equation 
7.15.  To  do  this,  the  approximation  is  made  that  all  CMs  have  the 
same  probability  of  survival  which  is  given  by  the  average  for  the 
survival  probabilities  of  all  the  CM's 

ps=  S^sl  +  IT  ps2 

Using  this,  the  probability  that  k  CMs  survive  the  FAD  is  now  given  by 

p(k  survive  FAD)  =  ( jj) pk(l-ps)C_k  7.16 
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Some  care  should  be  exercised  using  this  approximation  when  the 
number  of  fighter-interceptors  is  small  compared  to  the  number  of 
penetrators.  Since  this  program  was  written  largely  to  demonstrate 
the  validity  and  practicality  of  the  model,  the  approximation  is 
used  for  both  bombers  and  cruise  missiles  to  save  computer  time. 

A  flow  chart  of  the  program  is  shown  in  Figure  7.1.  The  program 
will  not  be  described  in  detail  in  this  paper.  A  listing  of  the  pro¬ 
gram  and  the  inputs  to  it  is  given  in  Appendix  A.  Some  verification 
and  results  of  the  program  are  given  in  Chapter  VIII. 
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FIGURE  7.1 

ow  chart  of  computar  program 


chao'ep  v : i : 

VERIFICATION  AND  RESULTS 

This  chapter  will  briefly  outline  some  of  the  steps  taken  to  verify 
the  computer  model.  Also  in  this  chanter  some  initial  results  (r 
the  model  are  given . 

VERIFICATION 

Any  computer  model  of  this  size  and  comp'exit;  -eguires  care^u' 
verification  to  insure  that  the  program  actually  joes  w^at  t*e  analyst 
wants  it  to  do.  Due  to  the  overall  time  constraint  placed  on  this 
project  and  various  problems  that  arose  with  the  computer  system  that 
was  used  the  program  may  not  have  received  the  extensive  verification 
it  deserves.  However,  some  initial  steps  have  been  taken  to  verify 
the  program. 

The  first  step  taken  to  verify  the  program  was  to  run  the  program 
for  an  extremely  scaled  down  version  of  the  problem.  This  same  problem 
was  also  done  manually  and  the  results  were  compared  to  those  from  the 
computer  version.  Specifically,  the  problem  involved  the  following 
basic  inputs: 

2  bombers,  initially  both  on  the  ground,  each  armed  with  one 
weapon; 

2  tankers,  both  initially  airborne,  and  a  bomber  probability  of 
successfully  refueling  given  a  tanker  of  1; 

2  FAD  interceptors  with  a  combined  probability  of  detecting,  con¬ 
verting  on,  and  killing  each  bomber  of  0.5; 

1  BSAM  site  with  a  PLD  of  half  the  width  of  the  barrier  (PLD*1) 
and  a  of  1  on  bombers  within  the  PLD;  and 

0  RASAM  sites. 
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The  intercept  intensity  parameter  •  was  set  equal  to  1  with  area 
FIs  having  a  p^  of  1.  The  detection  parameter  *  was  set  equal  to  0.5. 
The  velocity  of  a  bomber  was  1,  the  size  of  the  defended  area  was  2  by 
2  and  the  bombers  targets  were  located  at  a  depth  1  inside  the  defended 
area.  The  other  parameters  of  the  model  were  cnosen  to  give  the  bombers 
a  probability  of  0.5  of  successfully  escaping  their  base.  Note  that 
the  bombers  have  a  probability  1  of  successfully  refueling  in  this 
probl em. 

The  computer  program  and  the  hand  calculations  gave  the  same 
results  for  the  mean  and  variances  of  the  number  of  weapons  delivered 
and  the  number  of  bombers  that  survived  the  mission.  At  the  time 
this  initial  verification  was  done,  the  code  for  cruise  missiles  was 
not  yet  contained  in  the  program. 

It  should  be  noted  that  the  program  contains  several  different 
IF/THEN/ELSE  blocks.  In  any  simple  run,  such  as  they  are  described 
above,  many  small  parts  of  the  program  may  not  be  tested.  To  test 
every  part  of  the  program  this  way,  it  would  be  necessary  to  select 
values  of  the  inputs  that  will  insure  that  the  parts  of  the  program  in 
question  are  called.  The  time  constraints  placed  on  this  project  pre¬ 
cluded  such  in  depth  verification. 

The  second  major  form  of  verification  undertaken  for  this  program 
was  to  vary  some  of  the  major  inputs  of  the  model  and  verify  that  the 
changes  in  the  outputs  were  what  would  be  predicted. 

The  first  parameter  to  be  varied  was  the  number  of  bombers  in  the 
problem.  In  general,  as  the  number  of  bombers  is  increased,  all 
other  parameters  remaining  the  same,  the  numbers  of  weapons  delivered 


122 


and  bombers  surviving  both  increase.  However,  when  all  the  bombers 
are  initially  on  the  ground,  there  is  an  upper  limit  to  this  increase 
due  to  the  combination  of  the  rate  at  which  bombers  can  take  off  and 
the  time  in  which  they  can  take  off  before  the  SLBMs  arrive  at  their 
bases.  For  example,  if  the  missile  flight  time  is  twelve  minutes  and 
bombers  take  off  at  a  rate  of  one  per  minute,  there  is  nothing  to  be 
gained  by  increasing  the  number  of  bombers  initially  on  the  ground 
beyond  twelve. 

The  numbers  of  defenses  have  also  been  varied  with  the  program 
yielding  decreases  in  weapons  delivered  and  bombers  surviving  as  the 
numbers  of  defenses,  such  as  FAD  intercepters  and  SAMs,  are  increased. 
The  same  holds  true  when  the  intercept  intensity  paramter  X  is  in¬ 
creased,  which  is  the  same  as  increasing  the  level  of  area  FI  defenses 
in  the  problem. 

The  last  major  category  of  parameters  to  be  varied  included  the 
probabilities  of  detection  and  conversion  and  the  kill  probabilities 
for  the  different  defenses.  In  general,  two  results  are  expected. 

First,  as  the  probabilities  are  increased,  corresponding  to  improvements 
in  defensive  capabilities,  the  numbers  of  weapons  delivered  and  bombers 
surviving  should  both  decrease.  Second,  it  was  to  be  expected  that 
probabilities  near  0.5  will  give  the  largest  variances  for  the  results 
while  probabilities  near  0.1  and  0.9  will  give  smaller  variances.  Both 
of  these  expectations  are  born  out  by  the  program. 

Because  the  code  for  the  cruise  missiles  was  added  after  the  other 
code  had  been  written,  verification  of  the  outputs  for  CMs  was  less 
extensive  than  for  the  weapons  delivered  and  bombers  surviving. 
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However,  there  are  at  least  two  indications  that  the  CM  part  of  the 
program  works  as  it  should.  First,  as  the  number  of  CMs  in  the  pro¬ 
blem  is  increased,  the  expected  number  that  reach  and  destroy  targets 
also  increases.  Second  as  the  number  of  CMs  are  increased  the  expected 
numbers  of  weapons  delivered  and  bombers  surviving  also  initially 
increases.  Since  the  program  is  written  for  long  range  CMs,  increases 
in  the  number  of  CMs  increase  the  number  of  penetraters  entering  the 
FAD  zone,  initially  giving  each  one  a  better  chance  of  surviving.  After 
the  number  of  CMs  reaches  a  certain  level,  there  will  be  only  one  inter- 
cepter  assigned  to  each  bomber  and  all  the  rest  of  the  intercepters  will 
be  assigned  to  CMs.  Further  increases  in  the  number  of  CMs  will  not 
improve  the  number  of  bombers  surviving  or  weapons  delivered. 

RESULTS 

The  following  is  a  brief  overview  of  some  of  the  results  obtained 
from  the  model.  These  results  are  intended  to  show  that  the  model  works 
and  give  some  feel  for  the  output  that  can  be  obtained.  The  inputs 
have  been  selected  solely  for  this  purpose  and  are  not  necessarily 
meant  to  describe  any  specific  realistic  situation. 

The  base  case  for  all  of  the  results  to  be  presented  below  is  as 
follows: 

15  bombers,  initially  all  on  the  ground,  each  carrying  3  weapons; 

15  tankers,  initially  all  on  the  ground; 

30  FAO  fighter- intercepters; 

50  BSAMs;  and 

100  RASAMs. 
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The  intercept  intensity  paramter  >.  was  set  to  0.5.  The  size  of 
the  defended  area  was  4000  by  4000.  SAM  sites  had  a  PLD  of  3.  A 
bomber's  probability  of  refueling,  given  a  tanker  was  set  to  0.95.  All 
other  probabilities  were  set  to  0.5.  The  detection  parameter  a  was 
set  to  0.5  for  low  altitude  bombers  and  1.0  for  high  altitude  bombers. 
This  was  the  only  distinction  made  between  the  two  groups  of  bombers. 

To  reduce  the  amount  of  data  to  be  displayed,  and  because  the 
cruise  missile  code  had  received  less  verification  than  that  for 
bombers  and  weapons  released,  cruise  missiles  were  left  out  of  the 
probl em. 

The  results  of  the  test  case  were  as  follows: 

Expected  number  of  weapons  released  -  2.31 
Variance  of  number  of  weapons  released  -  9.10 
Expected  number  of  bombers  surviving  -  0.29 
Variance  of  number  of  bombers  surviving  -  0.33 
Four  different  parametric  variations  were  done  on  the  test  case. 
These  are  described  below. 

In  the  first  example,  the  number  of  weapons  on  each  bomber  was 
varied  between  1  and  5.  The  results  of  the  expected  number  of  weapons 
delivered  and  the  expected  number  of  bombers  that  survive  the  mission 
are  shown  in  Figure  8.1.  Two  things  can  be  noted.  First,  the  expected 
number  of  bombers  that  survive  the  mission  decreases  as  the  number  of 
weapons  on  each  bomber  increases.  This  is  because  the  bombers  must 
face  terminal  SAM  defenses  at  each  target  they  attack.  Also,  as  the 
number  of  weapons  on  each  bomber  increases,  the  rate  of  increase  of  the 
expected  number  of  weapons  delivered  decreases.  This  is  also  a  result 
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FIGURE  8.1 

Bombers  surviving  and  weapons  delivered  as  a  function  of 
the  number  of  weapons  on  each  bomber. 


In  the  second  example,  the  bombers'  probability  of  surviving  the 
TSAM  defenses  at  each  target  was  varied  from  0.1  to  0.9.  The  results 
are  shown  in  Figure  8.2.  As  expected,  more  weapons  are  delivered  and 
more  bombers  survive  the  mission  as  the  probability  of  surviving  the 
TSAMs  increase. 
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Bombers  surviving  and  weapons  delivered  as  a  function 
of  the  probability  of  surviving  TSAM  defenses. 


In  the  third  example,  the  intercept  intensity  parameter  A  is  varied 
from  0.1  to  0.9.  The  results  are  shown  in  Figure  8.3.  As  the  intercept 
intensity  increases,  more  FI  attacks  are  expected  on  a  bomber  per  unit 
time.  As  a  result,  the  expected  numbers  of  weapons  delivered  and 
bombers  surviving  decreases  with  increases  in  A. 

In  the  last  example,  all  of  the  probabilities  input  to  the  model 
(except  refueling)  are  varied  between  0.1  and  0.9.  Increases  in  the 
defensive  probabilities  of  kill  on  a  bomber  and  increases  in  a  bombers 
probability  of  surviving  the  TSAM  defenses  tend  to  have  offsetting  effects 
on  the  expected  numbers  of  bombers  surviving  and  weapons  released. 
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FIGURE  8.3 

Bombers  surviving  and  weapons  delivered  as  a 
function  of  intercept  intensity. 

;  example  was  chosen  to  show  the  effects  these  changes  will  have  on 
variances  of  the  model  outputs.  The  results  are  shown  in  Figure 


FIGURE  8.4 

Variances  of  model  outputs  asa 
function  of  the  input  probabilities. 
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This  has  just  been  a  small  sample  of  the  various  results  that 
can  be  obtained  with  this  model.  However,  they  should  be  sufficient 
to  indicate  that  the  model  works  and  to  demonstrate  some  of  its 
potential  power. 
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CHAPTER  IX 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  chapter  is  to  recap  the  model  that  has  been 
presented  in  this  paper  and  to  suggest  some  possible  areas  where 
further  work  could  produce  significant  improvment  or  extensions  to  the 
model.  At  the  end  of  this  chapter,  some  comments  and  conclusions 
about  this  project  will  be  made. 

SUMMARY 

This  paper  set  out  to  develop  an  analytic  model  of  strategic 
bomber  penetration  from  which  both  expected  value  and  variance  calcu¬ 
lations  on  the  results  could  be  made.  A  model  of  the  strategic  bomber 
mission  has  been  constructed  which  covers  the  entire  mission  from  the 
time  bombers  take  off  from  their  bases  until  they  leave  the  defended 
area  after  releasing  their  weapons.  It  was  recognized  from  the  out¬ 
set  that  the  scope  of  the  problem  would  ultimately  limit  the  depth- 
to  which  each  part  of  the  mission  could  be  modelled;  however,  it  is 
felt  that  the  current  model  contains  sufficient  detail  to  adequately 
represent  the  bomber  mission  for  many  purposes.  A  computer  program 
has  been  written  that  demonstrates  some  of  the  various  possible  out¬ 
puts  that  can  be  obtained  from  the  model.  In  particular,  the  program 
shows  that  the  model  can  give  variance  calculations  on  its  results, 
as  desired. 

There  probably  exists  no  model  of  this  scope  which  can  claim  to 
be  perfect  in  all  respects  and  the  same  applies  to  the  model  presented 
in  this  paper.  In  the  following  sections,  several  areas  are  suggested 
where  additional  work  could  yield  substantial  improvements  to  the 
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model.  In  some  cases  possible  approaches  or  potential  directions  for 
work  are  indicated;  the  other  areas  are  open  to  suggestions. 

COMPUTERIZATION 

The  computer  program  written  for  this  paper  was  designed  to  demon 
strate  the  validity  and  practicality  of  this  model.  The  orogram  will 
give  the  six  basic  results  described  in  Chapter  VII  but  as  it  current! 
stands,  it  will  not  output  any  of  the  numerous  intermediate  results  of 
the  model  that  may  be  of  interest  to  users. 

In  the  current  program,  it  is  not  easy  to  model  only  parts  of 
the  problem.  For  example,  by  inputting  an  initial  number  of  tankers 
in  the  air  that  is  greater  than  the  total  number  of  bombers  in  the 
problem  and  by  setting  the  probability  of  successfully  refueling  given 
a  tanker  to  1,  it  is  possible  to  insure  that  all  bombers  get  refueled. 
This  will  effectively  leave  the  refueling  stage  of  the  mission  out  of 
the  problem  as  far  as  the  final  results  are  concerned.  However,  the 
computer  will  still  run  the  refueling  routine  and  spend  the  time  cal¬ 
culating  all  the  various  probabl ities .  This  is  not  efficient.  In 
short,  the  computer  model  lacks  flexibility. 

Therefore,  for  this  model  to  achieve  its  full  potential  value, 
the  most  important  next  line  of  work  is  to  develop  a  flexible, 
efficient,  well  documented  program  for  general  use. 

INTERCEPT  DENSITY  FUNCTION 

A  major  part  of  this  model  rests  on  the  assumption  that  a  time 
independent  paramter  (\)  can  be  found  that  gives  the  average  number 
of  intercepts  possible  on  a  penetrating  bomber  per  unit  time.  In 
this  model,  this  parameter  has  been  assumed  to  be  input  and  minimum 
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consideration  has  been  given  to  the  derivation  or  selection  of  a 
value  for  the  parameter.  This  opens  several  possible  lines  of 
further  research. 

The  first  thing  that  should  be  done  is  to  find  a  derivation  for 
the  parameter  that  incorporates  as  much  information  about  the  defen¬ 
sive  fighter- i nterceptor  force  as  possible.  An  iterative  procedure  to 
find  a  time-dependent  function  for  the  parameter  has  already  been 
developed  by  D.  H.  Monahan  for  CCPEM-I  (Ref  12:3).  This  procedure 
is  described  in  Appendix  B.  It  should  be  possible  by  making  some 
appropriate  assumptions  to  modify  the  method  to  give  a  constant 
instead  of  time  dependent  parameter. 

There  is  an  obvious  alternative  approach  to  the  problem.  Since 
a  derivation  for  a  time  dependent  form  of  the  parameter  already 
exists,  and  since  there  are  good  reasons  to  believe  that  such  a  para¬ 
meter  would  in  fact  be  time  dependent,  it  might  be  worth  some  study 
to  determine  if  it  would  be  possible  to  modify  the  model  to  allow  for 
the  use  of  a  time  dependent  parameter.  The  major  disadvantage  to 
such  an  approach  lies  in  the  fact  that  the  penetration  probability 
function  derived  in  Chapter  V  would  no  longer  have  a  closed  analytic 
solution.  The  computer  model  would  have  to  be  expanded  first  to 
calculate  the  parameter  and  second  to  do  a  numerical  integration  of 
the  penetration  probability  function.  This  problem  could  only  be 
avoided  if  some  form  of  analytic  function  could  be  found  to  describe 
the  parameter  that  would  still  allow  the  integral  to  be  calculated 
analytical ly. 
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If  the  intercept  density  parameter  is  not  a  function  of  time, 
one  would  still  expect  it  to  depend  on  the  relative  numbers  of  pene- 
trators  and  fighters  in  the  air-to-air  battle  over  the  defeneded  area. 

In  this  paper,  the  parameter  has  been  implicitly  assumed  to  be  the 
same  for  all  possible  numbers  of  penetrators. 

Two  possible  approaches  to  this  problem  exist.  First,  an  expected 
initial  number  of  penetrators  in  the  air-to-air  battle  could  be  used 
to  calculate  the  parameter.  This  would  most  likely  give  offense  con¬ 
servative  results  for  the  total  expected  number  of  weapons  delivered 
in  the  attack.  However,  it  would  also  tend  to  understate  the  variance 
of  the  attack. 

The  second  possible  solution  to  this  problem  is  to  modify  the  FI 
part  of  the  model  to  permit  the  parameter  to  be  a  function  of  the 
number  of  penetrators  that  survive  the  BSAM  defense.  The  way  the 
current  model  is  written,  this  would  not  be  a  difficult  change  to  make. 
An  array  giving  the  parameter  as  a  function  of  the  number  of  penetrators 
could  be  input,  or,  if  an  analytic  approximation  to  the  function  could 
be  found,  the  computer  could  perform  the  calculations  during  each  run 
of  the  model  . 

WEAPONS  DELIVERED 

In  order  to  keep  the  scope  of  this  problem  down  to  something  that 
could  be  handled  in  the  time  frame  of  this  paper,  several  assumptions 
were  made  in  the  sections  on  terminal  SAMs,  weapons  delivered,  and 
value  extracted.  These  assumptions  were  made  for  modelling  purposes  and 
it  is  recognized  that,  from  the  point  of  view  of  realism,  these  are 
probably  among  the  least  defensible  assumptions  in  this  model. 
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For  this  reason,  it  is  suggested  that  a  major  area  for  further 
analysis  on  this  model  is  the  terminal  SAM  and  weapons  delivered  section 
This  would  probably  have  to  be  combined  with  some  form  of  weapon  alloca¬ 
tion  scheme  for  the  model  which  would  in  itself  be  a  significant  improve 
ment  to  the  model . 

The  Schultis  treatment  of  TSAMs  and  value  extracted  might  provide 
a  good  starting  point  for  such  work;  however,  it  will  have  some  limita¬ 
tions  if  one  tries  to  apply  it  directly  to  this  model.  For  one  thing, 
the  Schultis  model  gives  expected  value  damage  calculations  based  on 
saturation  of  the  terminal  defenses.  This  may  not  be  an  optimum 
strategy  in  the  first  place  and  if  it  is,  it  still  gives  no  indication 
of  the  variability  in  the  results,  which  was  one  of  the  major  goals 
of  this  model . 

Another  problem  that  should  be  addressed  is  the  relationship  be¬ 
tween  the  terminal  SAM  defenses  and  the  depth  of  the  location  of  the 
target  inside  the  defended  area.  In  the  Schultis  model,  targets  are 
defended  strictly  according  to  their  value  and  their  relative  location 
is  ignored.  This  may  not  be  the  actual  case. 

One  last  problem  to  consider  is  the  uncertainty  that  arises 
because  weapons  are  allocated  to  targets  before  the  mission  instead  of 
just  before  the  weapons  reach  targets.  The  Schultis  model  gives  cal¬ 
culations  to  show  how  much  damage  could  be  extracted  by  a  particular 
number  of  weapons  that  survive  all  the  other  defenses,  if  the  weapons 
were  free  to  be  allocated  at  that  point.  Unfortunately,  weapons  must 
be  allocated  to  targets  before  the  mission  begins  and  weapons  bound  for 
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high  value  targets  have  as  much  chance  of  being  lost  to  early  defenses 
as  weapons  destined  for  lower  value  targets.  This  needs  to  be  considered 
when  weapons  are  being  allocated  and  should  be  incorporated  in  this 
model  when  calculations  are  being  done  at  the  mean  and  variance  of  the 
damage  extracted  by  an  attack 

OTHER  AREAS  FOR  STUDY 

There  are  several  other  areas  in  which  work  could  be  done  to 
improve  this  model. 

Base  Escape.  Two  possible  suggestions  to  improve  the  model  repre¬ 
sentation  of  the  base  escape  problem  are  offered.  One  would  be  to 
change  the  geometry  employed  in  the  problem.  This  could  include  sub¬ 
stituting  some  more  likely  shape  for  the  volume  within  which  the  airborne 
bombers  will  be  located.  A  more  important  change  might  be  to  overlap 
the  lethal  volumes  from  various  warheads  targeted  on  a  base. 

The  second  suggestion  is  related  to  the  first.  The  model  implicitly 
assumes  that  all  the  bombers  that  take  off  from  the  base  before  the 
warheads  arrive  have  equal  chances  of  surviving.  In  reality,  the 
bombers  that  take  off  first  will  almost  surely  survive  while  a  bomber 
which  has  just  left  the  runway  when  the  warheads  begin  to  detonate  will 
have  little  chance  of  surviving.  One  possible  way  to  account  for  this 
would  be  to  fix  the  volume  within  which  the  detonating  warheads  will 
probably  lie  and  then  find  the  number  of  bombers  that  took  off  but  are 
still  inside  that  potentially  lethal  volume  of  airspace.  These  would 
be  the  most  recent  aircraft  to  take  off  and  their  survival  probabilities 
could  be  found  as  in  this  model.  The  aircraft  that  had  taken  off  first 
would  be  outside  the  space  where  the  warheads  would  detonate  and  so 
would  be  assured  of  surviving. 
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Refueling.  The  refueling  part  of  the  mission  is  a  complex  problem 
in  logistics  complicated  by  the  uncertainties  that  arise  when  one  models 
a  pre-emptive  SL8M  strike  against  the  strategic  bomber  and  tanker  forces. 
Very  little  work  has  been  done  in  the  past  toward  analytically  modelling 
this  problem.  This  paper  dealt  with  the  refueling  problem  in  greater 
depth  than  most  analytic  models;  however,  there  remains  room  for  greater 
detail  in  the  model  . 

3DMs.  Bomber  defense  missiles  were  not  considered  in  this  model. 

At  this  time,  that  is  probably  not  a  serious  drawback  to  the  model; 
however,  the  flexibility  of  the  model  would  be  increased  if  they  were 
added . 

Decoys.  Bomber  decoys  were  not  explicitly  included  in  this  model 
although  the  user  can  consider  them  when  selecting  a  value  for  the 
intercept  density  parameter  (See  Ref  5,  Appendix  A).  Penetrating 
bombers  may  find  the  use  of  decoys  advantageous  in  the  forward  air 
defense  stage  of  this  mission.  At  least  one  approach  to  analytically 
modelling  decoys  can  be  found  in  the  PENEX  model  (Ref  2). 

Command  and  Control.  In  this  model,  all  defensive  fighter- 

inteceptors  were  assumed  to  be  close  controlled.  This  requires  an 

effective  defensive  command  and  control  system.  This  model  does  not 

2 

consider  ways  on  which  the  C  system's  effectiveness  could  be  degraded 
or  the  course  of  a  mission.  There  are  at  least  three  possibilities. 

First,  false  tracks  may  be  generated  on  the  defensive  radars. 
Interceptors  would  be  assigned  to  these  false  tracks,  which  would  act 
and  could  be  modelled  similar  to  bomber  decoys. 
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Another  possibility  is  the  failure  of  a  C“  system  to  detect  a 
bomber.  In  this  model,  this  possibility  is  accounted  for  in  the 
battle  over  the  heartland  by  the  distribution  of  the  time  to  detection; 
however,  it  is  not  considered  in  the  forward  air  defense  battle.  The 
PENEX  model  demonstrates  one  approach  to  this  problem. 

The  third  possibility  is  the  saturation  or  even  destruction  of 
some  elements  of  the  C“  system.  This  is  simulated  in  the  APM  model 
but  a  good  analytic  approach  would  be  valuable. 

Approximations .  A  few  approximations  were  made  for  convenience 
on  this  model.  It  would  perhaps  be  worth  investigating  them  to  deter¬ 
mine  under  what  conditions  they  are  reasonable.  In  addition,  improved 
or  new  approximations  for  some  of  the  calculations  may  be  desirable. 
This  could  become  important  if  changes  to  the  model  result  in  large 
increases  on  the  amount  of  computer  time  required  to  run  it. 

Validation.  At  some  time  in  the  future,  perhaps  after  many  of  the 
changes  to  the  model  described  above  have  been  made,  it  might  be  worth¬ 
while  to  compare  some  of  the  results  of  this  model  to  a  large  simula¬ 
tion  such  as  the  APM.  This  would  be  an  excellent  way  to  validate  the 
model  and  determine  when  its  results  provide  a  reasonable  description 
of  the  bomber  penetration  mission. 

CONCLUSIONS 

This  paper  has  described  an  analytic  model  of  the  entire  bomber 
penetration  mission.  The  model  has  been  developed  to  allow  calcula¬ 
tions  of  the  variance  of  the  model  results. 
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Although  the  original  suggestion  for  this  paper  came  for  the 
Aeronautical  Systems  Division,  Deputy  for  Strategic  Systems,  it  is 
hoped  that  some  of  the  advantages  of  this  model  will  make  it  useful 
to  others  as  well.  The  model's  two  biggest  advantages  are  its  scope, 
in  terms  of  the  entire  bomber  penetration  mission,  and  its  capability 
to  calculate  variances  for  its  results. 

The  advantage  to  modelling  all  the  stages  of  the  bomber  penetra¬ 
tion  mission  is  that  it  allows  the  results  to  reflect  the  impact  of 
all  of  the  events  that  could  effect  the  outcome  of  the  mission.  For 
example,  in  a  model  that  only  describes  bomber  penetration  of  air 
defenses,  the  number  of  bombers  that  survive  the  base  escape  and  suc¬ 
cessfully  refuel  must  be  input.  This  will  not  reflect  the  uncertainty 
that  is  inherent  in  the  estimates  of  the  number  of  bombers  that  escape 
and  refuel.  With  this  model,  if  one  desires  to  use  the  model  for  just 
the  penetration  of  air  defenses,  that  can  be  done  by  excluding  the 
base  escape  and  refueling  portions. 

The  advantage  to  knowing  the  variance  of  the  results  of  a  model  is 
that  it  provides  an  idea  of  the  confidence  of  the  results.  Another  way 
of  thinking  of  it  is  that  the  variance  gives  an  indication  of  the 
spread  of  values  that  are  being  represented  by  the  model’s  expected 
value  result.  This  in  turn  will  provide  a  better  understanding  of  the 
significance  of  the  results. 

Potential  users  of  this  model  will  each  have  their  own  ideas  about 
how  the  model  most  effectively  could  be  used.  However,  to  conclude 
this  paper,  two  possible  uses  will  be  suggested. 
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First,  this  model  might  be  used  to  perform  excursions  for  studies 
where  the  APM  is  the  primary  model.  Because  of  the  APM's  long  run 
time,  it  is  not  practical  to  use  it  to  perform  parametric  variations 
or  sensitivity  analysis.  As  an  alternative,  some  of  the  numerous  out¬ 
puts  from  the  APM  could  be  used  as  inputs  to  this  model  .  This  model 
could  then  be  used  to  do  the  parametric  variations  and  sensitivity 
analysis  for  the  APM  results. 

The  second  way  in  which  this  model  might  be  used  is  in  conjunc¬ 
tion  with  a  force  on  force  attrition  model  such  as  the  Arsenal  Exchange 
Model.  This  model  could  provide  the  AEM  with  a  more  detailed  descrip¬ 
tion  of  the  air  breathing  element  of  the  strategic  triad  than  is 
currently  included. 
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APPENDIX  A 
COMPUTER  PROGRAM 

This  appendix  give  a  listing  of  the  program  written  for  this 
model.  Since  the  program  itself  contains  no  comments,  a  brief  outline 
of  the  program  is  presented.  In  particular,  the  purpose  of  each  of 
the  functions  and  subroutines  will  be  explained.  In  the  main  program, 
the  purposes  of  the  various  00  loops  will  be  outlined.  The  necessary 
inputs  to  the  program  are  also  defined.  The  program  is  written  in 
FORTRAN  5. 

INPUTS 

Values  of  the  following  variables  must  be  read  in  the  program 
(see  the  READ  statements  for  correct  order;  format  is  freefield). 

I  NBA  -  initial  number  of  bombers  airborne; 

INBG  -  initial  number  of  bombers  on  the  groud; 

INTA  -  initial  number  of  tanker's  airborne; 

INTG  -  initial  number  of  tankers  on  the  ground; 

VCB  -  bomber  rate  of  climb  during  base  escape; 

VHB  -  horizontal  velocity  of  bomber  during  base  escape; 

VCT  -  tanker  rate  of  climb  during  base  escape; 

VHT  -  horizontal  velocity  of  tanker  during  base  escape; 

NWP  -  number  of  weapons  on  each  bomber; 

CM  -  number  of  cruise  missiles  on  each  bomber; 

RB  -  bomber  take  off  rate; 

RT  -  tanker  take  off  rate; 

M  -  number  of  SLBM  warheads  targeted  at  base; 
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TMAX  -  maximum  time  in  which  bombers  can  take  off; 

TMIN  -  minimum  time  in  which  bombers  can  take  off; 

RL  -  lethal  radius  of  SL3M  warheads; 

PREF  -  bomber's  probability  of  refueling  given  an  available  tanker; 

F  -  number  of  FAD  fignter-intercepters ; 

WL  -  PLD  of  SAM  site; 

ilBSAM  -  number  of  SSAM  sites; 

I'iRSAM  -  number  of  RASAM  sites; 

W  -  width  of  defended  area; 

D-  depth  of  defended  area; 

LAMBDA  -  intercept  intensity  parameter;  and 

VCM  -  velocity  of  a  cruise  missile. 

In  the  following  inputs,  the  last  letter  of  each  variable  will  be 
H,  L,  or  C,  corresponding  to  bombers  at  high  and  low  altitudes  and 
cruise  missiles  respectively. 

pdcwx  -  probability  of  detection  and  conversion  for  FAD  intercepter; 

pdclx  -  probability  of  detection  and  conversion  for  area  fighter- 
intercepter 

pkwx  -  probabiity  of  kill  given  detection  and  conversion  for  FAD 
intercepter; 

pklx  -  probability  of  kill  given  detection  and  conversion  for  area 
fi ghter-intercepter ; 

pksamx  -  SAM  probability  of  kill  given  penetrater  enters  SAM 
coverage; 

ptsamx  -  probability  of  surviving  TSAM  defense  at  target; 

alphax  -  detection  parameter;  and 
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pw^  -  probability  that  a  weapon  that  is  released  successfully 
destroys  its  target. 

The  following  variables  correspond  to  target  locations  (see 
Chapter  VI ) . 

RH(I)  -  these  are  the  x .  for  bomoers; 

NCMR  -  number  of  intervals  defended  area  has  been  divided  into 
to  do  numerical  integration  over  location  of  cruise  missiles  target; 

LMR(I)  -  the  possible  locations  (depths)  of  the  cruise  missile 
targets;  and 

PCMR(I)  -  probability  that  a  cruise  missile  target  is  at  depth 
CMR(I) . 

OUTPUTS 

There  are  six  outputs  for  the  model.  (This  does  not  include  some 
of  the  basic  inputs  which  are  listed  out  at  the  end  of  each  run  for 
convenience). 

MEANW  -  expected  number  of  weapons  delivered; 

VARW  -  variance  of  number  of  weapons  delivered; 

MEANBS  -  expected  number  of  bombers  that  survive  the  mission; 

VARBS  -  variance  of  number  of  bombers  that  survive  the  mission; 

MEANCM  -  expected  number  of  cruise  missiles  that  reach  their 
targets;  and 

VARCM  -  variance  of  the  number  of  CMs  that  reach  their  targets. 
INTERMEDIATE  VARIABLES 

Many  intermediate  variables  are  calculated  in  this  program.  Some 
of  the  more  meaningful  ones  are  listed  here. 
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P3SAM^  -  probability  a  penetrater  survives  the  BSAM  defense; 

PBD^  -  probability  a  weapon  is  released  over  a  target; 

PREC^  -  probability  a  bomber  survives  to  leave  the  defended  area 
and  recover;  and 

PCM  -  probability  a  CM  reaches  its  target. 

FUNCTIONS 

There  are  four  external  functions  employed  in  this  program. 

CGMB(N,K)  -  calculates  ( ^ ) ; 

XT0P(X,P)  -  finds  ,  with  0°  defined  as  1  for  purposes  of  the 
binomial  distribution; 

FP( X)  -  probability  density  function  for  the  Beta  distribution 
with  parameters  a=5  and  3=1;  and 

FE(X)  -  used  to  find  the  expected  value  of  the  Beta  distribution 
over  a  particular  interval. 

SUBROUTINES 

There  are  five  subroutines  called  in  the  program. 

TIME.  SUBROUTINE  TIME  takes  TMAX  and  TMIN  and  fits  a  Beta  distri¬ 
bution  between  them.  It  then  finds  the  distribution  for  the  number  of 
aircraft  (bombers  or  tankers)  that  can  become  airborne  before  the  SLBM 
warheads  at  their  base.  It  also  finds  the  expected  time  the  bombers 
had  in  which  to  take  off  given  the  number  of  bombers  that  did  take  off. 
(See  Chapter  IV,  the  section  on  base  escape.) 

ESCAPE.  SUBROUTINE  ESCAPE  takes  the  distribution  of  the  number  of 
aircraft  (bombers  or  tankers)  that  can  get  off  the  ground  during  an 
SLBM  attack  and  finds  the  probability 'that  a  specified  number  of  air¬ 
craft  survive  the  base  escape  problem. 
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REFUEL.  SUBROUTINE  REFUEL  takes  the  distribution  for  the  number 
of  tankers  that  have  survived  the  base  escape  problem  and  finds  the 
probability  that  a  specified  number  of  bombers  successfully  refuel  for 
a  given  number  of  bombers  that  have  survived  the  base  escape. 

FAD.  SUBROUTINE  FAD  finds  penetraters'  probabilities  of  surviving 
the  forward  air  defense.  In  the  case  of  bombers,  FAD  will  find  the 
probability  that  a  specified  number  of  bombers  at  a  given  altitude 
(high  or  low)  survive  given  the  total  number  of  penetraters  entering 
the  forware  air  defense.  For  cruise  missiles,  FAD  finds  the  probability 
that  any  cruise  missile  survives  the  forward  air  defense,  given  the 
total  number  of  penetraters. 

BOMBS.  SUBROUTINE  BOMBS  finds  the  probability  that  a  weapon  is 
released  over  a  target. 

MAIN  PROGRAM 

In  this  section,  the  purpose  of  each  CO  loop  will  be  briefly 
explained.  DO  loops  will  be  referred  to  by  their  CONTINUE  statement 
numbers.  The  index  variable  for  the  loop  may  be  included  in  parentheses. 

5,15  -  variable  initialization  loops; 

50  -  finds  the  probability  a  LRCM  penetrates  the  area  defenses  to 
reach  its  target; 

60,70  -  find  the  distribution  of  the  number  of  bombers,  including 
those  initially  on  airborne  alert,  that  survive  the  base  escape;  and 

80,90  -  find  the  distribution  for  the  number  of  tankers  that  sur¬ 
vive  the  base  escape. 

100(N)  -  the  number  of  weapons  released,  bombers  surviving,  and/or 
cruise  missiles  that  reach  targets.  N  will  be  incremented  from  1 
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to  the  maximum  of  the  total  number  of  weapons  or  cruise  missiles  in 
the  problem.  Fore  each  N,  the  probability  that  N  weapons  are  released, 
that  N  bombers  survive,  and  that  N  cruise  missiles  reach  their  targets 
will  be  calculated,  resulting  in  distributions  for  all  three  outputs. 

2Q0(JE)  -  conditions  the  distributions  described  above  on  the 
number  of  bombers  ( JE )  that  survive  the  base  escape;  and 

300(JR)  -  conditions  the  distribution  on  the  number  of  bombers 
that  successfully  refuel. 

400,500  -  for  given  JE  and  JR,  a  distribution  for  the  number  of 
low  altitude  bombers  that  survive  both  the  FAD  and  BSAM  defenses  is 
found  and  stored  in  array  PSPEN . 

600,700  -  for  given  JE  and  JR,  a  distribution  for  the  number  of 
high  altitude  bombers  that  survive  the  FAD  and  BSAMs  is  found  and 
stored  in  the  array  PSPEN. 

800,900,1000  -  find  the  probability  that  a  total  of  N  weapons  are 
delivered  by  the  high  and  low  altitude  bombers.  This  probability  is 
conditioned  on  the  number  of  weaons  released  by  the  low  altitude  bombers. 

100,1200,1300  -  find  the  probability  that  a  total  of  N  bombers  from 
both  the  high  and  low  altitude  groups  survive  the  mission.  The  proba¬ 
bility  is  conditioned  in  the  number  of  low  altitude  bombers  that  sur¬ 
vive. 

Two  last  points  will  be  made  about  the  main  program.  The  first 
concerns  lines  1620,  1670  and  1740.  These  lines  provide  a  time 
saving  approximation  to  the  program.  Although  most  of  the  calculations 
in  the  program  are  DOUBLEPRECISION,  it  will  frequently  be  the  case  that 
accuracy  to  one  or  two  decimal  places  will  be  sufficient.  In  this 
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program,  if  the  probability  that  a  particular  number  of  bombers  sur¬ 
vive  the  base  escape  is  very  small,  subsequent  calculations  in  the 
200  loop  will  be  skipped  since  they  will  not  contribute  any  significant 
change  to  the  final  result.  The  300  loop  will  be  skipped  if  the 
probability  that  a  particular  number  of  bombers  refuel  is  too  small 
to  contribute  to  the  final  results.  The  way  the  program  is  written, 
one  decimal  place  accuracy  can  be  had  with  up  to  100  bombers  in  the 
problem.  In  the  results  presented  in  Chapter  VIII,  where  only  15 
bombers  were  modelled,  the  accuracy  was  better  than  two  decimal  places. 
If  time  is  not  a  major  consideration  in  the  program,  these  lines  can 
be  deleted.  As  the  program  is  improved,  better  approximations  may  be 
found . 

The  second  point  to  be  made  concerns  the  cruise  missile  code.  It 
should  be  noted  that  the  probability  that  a  CM  fails  to  reach  its  tar¬ 
get  due  to  nondefense  related  reasons  (qc,  see  Chapter  VI)  has  not 
been  explicitly  included,  implying  that  it  is  always  1.  This  can  be 
remedied  by  reading  a  QC  and  altering  line  2300  to  read: 

2300=  PCM=PCMA*PRI*PB$AMC*PTSAMC*( 1-QC) 

This  concludes  the  description  of  the  computer  program  used  in 
this  model.  This  appendix  is  not  meant  to  provide  detailed  documenta¬ 
tion  of  the  program.  However,  it  is  hoped  that  with  the  help  of  this 
description,  the  program  can  be  understood  by  someone  familiar  with  the 
model  and,  in  particular,  the  integration  of  the  model  described  in 
Chatper  VII . 

A  listing  of  the  program  follows. 
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140=  PROGRAM  MODEL 

150=  INTEGER  BT»F,H»I» INBA , INBG » INTA » INTG ,J»JE,JR»K»L»M»N 

160=  INTEGER  NBSAM  * NRSAM » NWP 

170=  INTEGER  C » CM , NCMR 

180-  REAL  ALPHAH , ALPHAL , D , LAMBDA , PDCLH , PDCLL , PDCUH » PDCWL 

190=  REAL  PNLH  , PKLL » PNWH , PNWL  , PKS AMH » PKSAML  , PREF 

200*  REAL  RB»RL ,  RT » TMAX  ,  TMIN  f  VCB  ,  VCT  ,  UHB, UHT  , W, WL , VBH  ,  UBL 

210=  REAL  PDCWC , PDCLC  ,  PNWC  »  PNLC » PKSAMC 

220*  REAL  PUH»PUL 

230=  REAL  CMR< 1  ♦  10) »PCMR< 1*10) 

240=  DOUBLEPRECISION  PTSAMC » PTSAMH , PTSAML 

250=  DOUBLEPRECISION  AfFNC 

260=  DOUBLEPRECIS  I  ON  MEANBS  ,  MEANW  ,  VARBS .  VARW  .•  PBSAMH  r  PBSAML  ,PBDH  , 

PBDL 

270=  DOUBLEPRECISION  PH , PL , P JR  * P JE ,  PNB , PNW ,  PRECH » PRECL ,  P Rif  SUM 

280=  DOUBLEPRECISION  SUMA  , SUMB ,  SUMC , SUMD » SUMBEB ,  SUMBEW » SUMRB , SUM 

RWnSUMT 

290=  DOUBLEPRECISION  CA , CB » CC , PCM , PCMA 

300=  REAL  RN<0*25) 

310=  DOUBLEPRECISION  PSPEN<0J201> 

320=  DOUBLEPRECISION  BESCCOt  100)  ,TESC<0*.200) 

330=  DOUBLEPRECISION  BET<0 i 100 ) , BAIR ( 0 J 100 ) t TET ( 0 J 200 ) t TAIR ( 0 J 20 

0) 

340=  DOUBLEPRECISION  FP » FE 

350=  DOUBLEPRECISION  SUME,SUMF»SUMBEC,MEANCM,VARCM»PBSAMC 

360=  DOUBLEPRECISION  X » Y » Z , COMB , XTOP 

370=  READ*,  I NBA , I NBG , VCB , UHB , RB » UBH , UBL » NWP 

380=  READ*,  INTA, INTG, VCT,VHT,RT» PREF 

390*  READ*,  PUH , PUL 

400=  READ*,  M,RL, TMAX, TMIN 

410=  READ*,  F, PDCLH, PDCLL fPDCWH, PDCWL, PKLH, PKLL, PKWH,PKWL 

420=  READ*,  U,D 

430=  READ*,  NBSAM, NRSAM, PKSAMH , PKSAML, PTSAMH, PTSAML, WL 

440*  READ*,  ALPHAC, ALPHAH, ALPHAL, LAMBDA 

450=  -READ*,  <RN< I ) ,1=1 ,NWP ) 

460=  READ*,  PDCWC, PDCLC, PKWC,PKLC, PKSAMC, PTSAMC 

470*  READ*,  CM»NCMR» VCM 

480*  READ*,  <CMR<X> ,1=1, NCMR) 

490=  READ*,  <PCMR< I ) , 1=1 , NCMR ) 

500*  PRINT' <A> ' , ' 1 ' 

510*  PRINT '</A,I6,A,F10»4,AfF10»4,A,I10)'f 

520=  1  '  INBA*' , INBA, '  VCB=',VCB,'  VBL=',VBL»'  NWP=',NWP 

530*  PRINT' (/AfI6»3<A»F10.4)>'f 

540-  l  '  INBG*' , INBG, '  VHB-'fVHB,'  VBH=',VBH,'  RB=',RB 

550*  PRINT' </A,I6,A,F10.4,A,F9.4) ' , 

560*  J  '  INTA-', INTA,'  VCT-'fVCTf'  PREF*'fPREF 

570=  PRINT' (/AfI6,A,F10.4fAfF11.4>  ' , 

580=  J  '  INTG-' f INTG f '  VHT*',VHT,'  RT»'fRT 

590=  PRINT' </2<A,F10.6))'f'  PWH»'fPWHf'  PWL-'fPWL 

600=  PRINT' </A, 19) ' f '  F-'fF 

610*  PRINT' </4<A,F8.6) ) ' » 

620=  J  '  PDCWL*' , PDCWL, '  PDCLL* ', PDCLL , '  PKWL=' fPKWLf '  PKLL= 

' fPKWL 


V 
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6 30= 
640= 
' ,PKLH 
650= 
660= 
670= 
680= 
690= 
700= 

- 710=- 

720= 

730= 

740= 

750= 

760= 

770= 

780= 

790= 

800= 

810= 

820= 

830= 

840= 

850= 

860= 

870= 

880= 

890= 

900= 

910= 

920= 

930= 

940= 

950= 

960= 

970= 

980= 

990- 

lOCO- 

1010= 

1020= 

HAHf 

1030= 

1040= 

HALr 

1050= 

1060= 

1070- 


PRINT'  </4<A»F3. 6)  )  '  t 

:  '  PDCWH= ' t PDCWH  r '  PDCLH= ' » PDCLH » '  PKWH=' f PKWHf '  PKLH= 

PRINT' < /A » FI 1.6) 'f '  WL= ' » WL 
PRINT' ( /Af 16  f3 ( AfFS  *  6) )  '  f 
J  '  NBSAM=' fNBSAM* '  PKSAML= ' > PKSAML t 

\  '  PTSAML=' fPTSAMLf '  PNSAMC= ' f PKSAMC 

PRINT' </A> 16  f 3< Af F8.6 ) ) ' » 

J  '  NRSAM=' f NRSAM f '  PKSAMH= ' t PKSAHH , 

:  '  PTSAMH=' fPTSAMHf '  PTSAMC=' fPTSAMC 

PRINT' </AfI9fAfF10.5) ' f 
:  '  M='fMf'  TMAX* ' ? TKAX 

PRINT' </AfF3.2fAfF10»5> '  » 

1  '  RL= ' » RL  f '  TMIN= ' f TMIN 

PRINT' </2<A»F10.5) ) ' f 

1  '  ALPHAL=' fALPHALf '  ALPHAC= ' f ALPHAC 

PRINT' </2CAfF10.5) > ' f 

t  '  ALPHAH=' fALPHAHf '  LAMBDA= ' f LAMBDA 

PRINT' </AfF11.4fAfF12.4> ' r 
J  '  W='*Wr'  D='fD 

PRINT' </Af I10fAfF10.4) ' f 
1  '  CM=' »CM* '  VCM= ' » UCM 

PRINT' </2<A*F10.4) ) ' t '  PDCWC=' fPDCWCf '  PKWC=' fPKWC 
PRINT' <///A) ' f '  ' 

H=0 

L=1 

C=2 

MEANM=0 . 

MEANBS=0. 

VARBS=0. 

VARW=0« 

MEANCM=0 . 

VARCM=0 . 

BT-INBA+INBG 
X=1.-PKSAML*WL/W 
PBSAML=XTOP  <  X  f  NBSAM ) 

Y=1.-PKSAMH*WL/W 
PBSAMH=XTOP <  Y . NBSAM ) 

Z=1 « -PKSAMC*WL/W 
PBSAMC=XTOP ( Z  ? NBSAM ) 

CALL  BOMBS  <  NWP f  PKSAMH  f  WL  f  W  f  D  f  RN  f NRSAM f PDCLH* PKLH f LAMBDA f ALP 
1  VBH  t PTSAMH f  PBDH  > 

CALL  BOMBS  <  NWP  *  PKSAML f  WL  f  W  f  D  f  RN  f NRSAM f PDCLL f PKLL *  LAMBDA  *  ALP 

:  vblfPtsamlfPbdd 

PBDH»PBDH*PWH 

PBDL=PBDL*PWL  -  •  ~  ■  r.--/ 


1080=  SUMA=  <  PDCLL*PKLL*LAMBDA/  <  PDCLL*PNLL*LAMBDA-ALPHAL )  ) 

1090=  J  *EXP<-ALPHAL*D/UBL) 

1100=  SUMB=< ALPHAL/ <  PDCLL*PKLL*LAMBDA-ALPHAL ) ) 

1110=  l  *EXP<-PDCLL*PNLL*LAMBDA*D/VBL) 

1120=  SUMC=XT0P (Xf NRSAH ) 

1130=  PRECL=SUMC*  <  SUMA-SUMB >#XT0P(PTSAMLr  NUP ) 

1140=  SUMA=  <  PDCLH*PKLH#LAMBDA/ < PDCLH*PKLH*LAMBDA-ALPHAH) ) 

1150=  t  *EXP<-ALPHAH*D/VBH> 

1160=  SUMB= ( ALPHAH/ ( PDCLW*PKLH*LAMBDA-ALPHAH ) > 

1170=  i  *£XP<-PDCLH*PNLH*LAMBDA*D/VBH> 

1180=  SUMC=XTDP  <Y»NRSAM) 

1190=  PRECH=SUMC* ( SUMA-SUMB ) *XT0P ( PTSAMH » NUP ) 

1200=  DO  5  1=0 » 100 

1210=  BAIR ( 1 ) =0 . 

1220=  BET  < I )  =  1 . 

1230=5  CONTINUE 

1240=  DO  15  I =0  f 200 

1250=  TAIR(I)=0. 

1260=  TET ( I ) =1 . 

1270=15  CONTINUE 

1280=  CALL  TIME<TMAXrTMIN»RB*BAIR*BET  r INBG) 

1290=  CALL  TIME < TMAX »TMIN»RT»TAIRf TET » INTO) 

1300=  SUM=0 . 

1310=  DO  50  1=1 >NCMR 

1320=  CA=<PDCLC*PKLC#LAMBDA/<PDCLC*PKLC*LAMBDA-ALPHAL> ) 

1330=  J  *EXP<-ALPHAL*CMR<I)/VCM> 

1340=  CB= <  ALPHAL/ ( PBCLC*PKLC*LAMBDA-ALPHAL ) ) 

1350=  :  *EXP<-PDCLC*PKLC*LAHBDA*CMR<I>/UCH> 

1360=  2=1 . -PKSAMC*UL*CMR< I ) / ( W#D > 

1370=  CC=XTOP  <  Z  r  NRSAM ) 

1380=  SUM=SUM+CC*<CA-CB>*PCMR(I> 

1390=50  CONTINUE 

1400=  PCMA»SUM 

1410=  DO  60  1=0 » INBA-1 

1420=  BESC<I)=0. 

1430=60  CONTINUE 

1440=  DO  70  !=INBA» INBA+INBG 

1450-  K=I-INBA 

1460=  CALL  ESCAPE <K» INBG»M»RL» BAIR » BET f VHB»UCB»BESC< I) ) 

1470-70  CONTINUE 

1480=  DO  80  I=0»INTA-1 

1490=  TESC< I )=0. 

1500-80  CONTINUE 

•1510=  DO  90  I=INTAf INTA+INTG 

1520=  K-I-INTA 

1530=  CALL  ESCAPE <K» INTGrMrRLr TAIR> TET» VHTr VCT»TESC< I ) ) 

1540=90  CONTINUE 

1550=  SUMA-0 • 

1560=  SUMB-0 • 

1570=  SUMC-O. 

1580=  SUMD=0. 

1590=  SUME=0 . 

1600=  SIJMF=0. 
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.1  -- 


1610= 

1620= 

1630= 

1640= 

1650= 

1660= 

1670= 

1680= 

1690= 

1700= 

1710= 

1720= 

1730= 

1740= 

1750= 

1760= 

1770= 

1780= 

1790= 

1800=500 

1810= 

1820=400 

1830= 

1840= 

1850= 

1860= 

1870= 

1880=700 

1890= 

1900=600 

1910= 

1920= 

1930- 

1940= 

1950= 

PEN< J) 
1960=900 
1970= 
1980= 
1990= 

-  2000=  — 

+  I>* 

2010- 


DO  100  N=1»MAX<BT*NWP»BT*CM) 

A=1./<1000000.*N*N) 

SUMBEB-0. 

SUHBEU=0. 

SUMBEC=0» 

DO  200  JE  =  INBA r BT 
PJE=BESC< JE) 

IF(PJE.LE.A)  GO  TO  199 
SUMRB-O. 

SUMRW-O. 

SUMRC-0 « 

DO  300  JR=0 1 JE 

CALL  REFUEL < JR » JE  , INT A  r INTG  t TESC  r PREF  > P JR ) 

IF (PJE#PJR«LE  .A)  GO  TO  299 
DO  400  J  =  0 » JR 
SUM=0 . 

DO  500  I=J» JR 

CALL  FAD ( I f  JE » JR » JE*CM » L  >  F  > PDCUL  r PKWL » PRI  > 

SUM=SUM+COMBCI» J)*XTOP<PBSAML> J)*XTOP< 1 .-PBSAMLr I-J>*PRI 

CONTINUE 

PSPEN< J>=SUM 

CONTINUE 

DO  600  J=0» JE-JR 
SUM=0. 

DO  700  I=J» JE-JR 

CALL  F AD < I » JE f  JR » JE*CM »H»F» PDCUH » PKUH , PR I ) 

SUM-SUM+COMB  < I » J) *XTQP  <  PBSAMH f  J ) *XTOP ( 1 . -PBSAMH ,  I- J  > *PRI 
CONTINUE 

PSPEN< JR+J+1)=SUM 

CONTINUE 

SUMT=0» 

DO  800  I=0»JR*NWP 
SUM=0. 

DO  900  J=0» JR 

SUM-SUM+COMB  <  J»NWP > I ) *XTOP  <  PBDL » I ) *XTOP  < 1 . -PBDL f  J*NWP- 1 >  *PS 

CONTINUE 
PL=SUM 
SUM-0 . 

DO  1000  J=0» JE-JR 

SUM=SUM+COMB ( J*NWP  t  N- 1 ) *XTOP ( PBDH t N- 1 ) *XTOP  < 1 . -PBDH » J*NWP-N 
1  PSPEN<  JR+ J+l ) 


2020=1000 

2030- 

2040= 

2050-800 

2060- 

2070- 

2080= 


CONTINUE 

PH-SUM 

SUMT»SUMT+PL*PH 

CONTINUE 

PNW-SUMT 

IF<N«8T« JE)  00  TO  298 
SUMT-O. 
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2090= 

2100= 

2110= 

2120= 

2130=1200 

2140= 

2150= 

2160= 

2170= 

2180= 

2190=1300 

2200= 

2210= 

2220=1100 

2230= 

2240= 

2250=298 

2260= 

2270=299 

2280=300 

2290= 

2300= 

2310= 

2320= 

2330= 

2340= 

2350=199 

2360=200 

2370= 

2380= 

2390= 

2400- 

2410= 

2420= 

2430- 

2440- 

2450= 

2460- 

2470= 

2480= 

2490= 

2500* 

2510= 

2520-100 

2530- 

2540- 

2550= 

2560- 

2570= 

2580= 


DO  1100  1=0 » JR 
SUM=0 » 

DO  1200  J=I f  JR 

SUM-SUM+COMB  <  J » I )  *XTOP  <  PRECL » I >  *XTOP  <  1 .  -PRECL » J-I  >  *PSPEN  ( J ) 

CONTINUE 

PL=SUM 

SUM=0. 

DO  1300  J=0 1 JE-JR 

SUM-SUM+COMB < J » N- 1 ) *XTOP  < PRECH  t N- 1 ) *XTOP ( 1 . -PRECH  f J-N+I ) 
t  *PSPEN< JR+J+1 ) 

CONTINUE 

PH=SUM 

SUMT=SUMT+PL*PH 

CONTINUE 

PNB=SUMT 

SUMRB=SUMRB+PNB*PJR 

CONTINUE 

SUMRW=SUMRW+PNW*PJR 

CONTINUE 

CONTINUE 

CALL  FAD<0»JE»0»  JEfrCM*  C » F  » PDCUC  t PKWC » PRI ) 
PCM=PCMA*PRI*PBSAMC*PTSAMC 

PNC=COMB  <  JE*CM » N ) *XTOP ( PCM , N ) *XTOP ( 1 . -PCM » JE*CM  - N  ) 
SUMBEB=SUMBEB+SUMRB*PJE 
SUMBEW=SUMBEU+SUMRU*PJE 
SUMBEC=SUMBEC+PNC*PJE 
CONTINUE 
CONTINUE 
IF<N.LE.BT*NUP) 
l  THEN 

SUMA=SUMA+N*SUMBEU 

SUMB=SUMB+N*N*SUMBEU 

ENDIF 

IF(N.LE.BT) 
l  THEN 

SUMC=SUMC+N*SUMBEB 

SUMD=SUMD+N#N*SUMBEB 

ENDIF 

IF<N.LE*BT*CM) 
t  THEN 

SUME«SUME+N*SUMBEC 

SUMF=SUMF+N*N*SUMBEC 

ENDIF 

CONTINUE 

HEANU-SUMA 

VARW=SUMB-SUMA*SUMA 

MEANBS-SUMC 

VARBS»SUHD-SUMC*SUMC 

MEANCM-SUME 

VARCM»SUMF-SUME*SUME 
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— 2590= 
2600= 
2610= 
2620= 
2630= 
2640= 
2650= 
2660= 
2670= 
2680= 
2690= 
2700= 
2710= 
2720= 
2730= 
2740= 
2750= 
2760=10 
2770= 
2780= 
2790= 
2800= 
2810=20 
2820= 
2830= 
2840= 
2850=30 
2860= 
2870=40 
2880= 
2890= 
2900= 
2910= 
2920= 
2930= 
2940= 
2950= 
2960= 
2970= 
2980= 
2990= 
3000- 
3010= 
3020= 
3030- 
3040= 
3050= 
3060= 


PRINT' <///AfF25. 16) ' f '  MEANW= ' fMEANW 
PRINT' </AfF26.16> ' f '  VARW=' fVARW 
PRINT' </AfF24. 16) ' f '  MEANBS= ' fMEANBS 
PRINT' (/AfF25. 16) ' f '  UARBS= ' t VARBS 
PRINT' </AfF24. 16) ' f '  MEANCM=' fMEANCM 
PRINT ' ( /A f F25 ► 16 ) ' f '  VARCM= ' f VARCM 
END 

FUNCTION  COMB(NfK) 

INTEGER  IfKfN 
DOUBLEPRECISION  CfCOMBfR 
IF<K.LT.O.OR.K.GT.N)  GO  TO  30 
C=l. 

IF(K.LE.N-K) 

THEN 

DO  10  1  =  1 f  K 
R=( I+N-K)/I 
C=C*R 
CONTINUE 
ELSE 

DO  20  1=1 fN-K 
R=( I+K)/I 
C=C*R 
CONTINUE 

ENDIF 
COMB=C 
GO  TO  40 
CONTINUE 
C0MB=0. 

CONTINUE 

END 

FUNCTION  XTOP(XfP) 

INTEGER  P 

DOUBLEPRECISION  XfXTOP 
IF<P.LE.O) 

THEN 

XTOP-1. 

ELSE 

XTOP=X**P 

ENDIF 

END 

FUNCTION  FPU) 

DOUBLEPRECISION  XrFP 
FP»42.*<  <X**6. )/6.-<X*#7. )/7. ) 

END 

FUNCTION  FE  <  X ) 

DOUBLEPRECISION  XfFE 
FE=42.*<  <X**7. >/7»-<X**8. >/8, > 

END 
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3070- 

3080= 

3090= 

3100= 

3110= 

3120= 

3130= 

3140= 

3150* 

3160= 

3170= 

3180= 

3190= 

3200= 

3210= 

3220= 

3230= 

3240= 

3250= 

3260= 

3270= 

3280= 

3290= 

3300= 

3310= 

3320= 

3330= 

3340= 

3350= 

3360= 

3370= 

3380= 

3390= 

3400=105 

3410= 

3420= 

3430=205 

3440= 

3450= 

3460= 

3470* 

3480= 

3490= 

3500= 

3510= 

3520= 

3530=211 

3540= 

3550= 

3560= 

3570= 

3580= 

3590= 

3600= 

.3610= 


SUBROUTINE  TIME ( TMAX . TMIN > R > AIR , PET > INPG ) 

INTEGER  N » INPG 
REAL  TMAX  >  TMIN » R 

DOUBLEPRECISION  AIR < 0 t * ) t PET ( 0 J * >  f  ANfFEfFPfPTfET 
DOUBLEPRECISION  XfY 
IF<  TMAX .EQ.TMIN) 

THEN 

N=TMIN*R 
N=MIN<Nf INPG) 

AIR(N)=1 . 

PET  <  N )=TMIN 

ELSEIFC  TMAX-TMIN » LE . 1 . /R . AND . TMAX . NE .TMIN) 
THEN 

N=TMIN*R 

AN=N 

IF (N. GE . INPG )  GO  TO  105 
IF(AN+1.0.GE.TMAX*R> 

THEN 

AIR  ( N ) =1 . 

PET(N)=TMIN+0 . 75* ( TMAX-TMIN ) 

ELSE 

X=<  <AN+1. >/R-TMIN)/(TMAX-TMIN) 
PT=FP<X) 

ET=FE<X)/PT 

. — - . —  AIR<N)=PT 

PET  <N)=TMIN+ET*(TMAX-TMIN) 

X=<  <AN+1. >/R-TMIN)/(TMAX-TMIN) 

PT=1 »-FP<  X ) 

ET=<0.75-FE(X) )/PT 
AIR(N+1 ) =PT 

PET (N+l >=TMIN+ET*(TMAX-TMIN) 

ENDIF 
GO  TO  205 
CONTINUE 
AIR< INPG)=1 ♦ 

PET  < INPG  > -TMIN+0 . 75* ( TMAX-TMIN) 

CONTINUE 

ELSE 

N=TMIN*R 

AN=N 

IF (N . GE . INPG)  GO  TO  305 

X=< <AN+1. >/R-TMIN)/<TMAX-TMIN> 

PT=FP<  X ) 

ET=FE<X)/PT 

AIR(N)=PT 

PET <N ) =TMIN+ET* < TMAX-TMIN ) 

CONTINUE 

N-N+l 

AN=N 

IF<N*EQ * INPG)  GO  TO  405 

IF<  < AN+1 . ) /R .GE . TMAX )  GO  TO  111 

X»< (AN+1 . >/R-TMIN)/(TMAX-TMIN> 

Y=<AN/Rr-TMIN)/<TMAX-TMIN> 

PT =FP ( X  > -FP ( Y ) 

ET*<FE<X)-FE< Y>  >/PT 
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3620= 

3630= 

3640= 

3650=111 

3660= 

3670= 

3680= 

3690= 

3700= 

3710= 

3720=305 

3730= 

3740= 

3750= 

3760=405 

3770= 

3780= 

3790= 

3800= 

3810= 

3820=505 

3830= 

3840= 

3850= 

3860= 

3870= 

3880= 

3890= 

3900= 

3910= 

3920= 

3930= 

3940=101 

3950= 

3960= 

3970= 

3980= 

3990= 

4000= 

4010= 

4020= 

4030= 

4040= 

4050= 

4060= 

4070= 

4080=102 

4090= 

4100= 

4110= 

4120= 

4130=202 

4140= 

4150= 

4160= 

4170= 


AIR  <  N ) =PT 

PET ( N ) =  TMIN+£T*  <  TMAX-TMIN ) 

GO  TO  211 
CONTINUE 

Y=  <  AN/R-TMIN ) / <  TMAX-TMIN) 

PT=1 »-FP(Y) 

£T=<0. 75-FE ( Y ) )/PT 
AIR  ( N )  =F'T 

PET <N)=TMIN+ET*< TMAX-TMIN) 

GO  TO  505 
CONTINUE 
AIR ( INPG ) =1 . 

PET (INPG)=TMIN+0. 75* (TMAX-TMIN) 

GO  TO  505 
CONTINUE 

Y=< AN/R-TMIN)/ (TMAX-TMIN) 

PT=1 . -FP  <  Y ) 

ET=<0.75-FE(Y) )/PT 
AIR ( N ) =PT 

PET <N)=TMIN+ET*< TMAX-TMIN) 

CONTINUE 

ENDIF 

END 

SUBROUTINE  ESCAPE ( JE t INPG fM , RL f AIR f PET fUH f VC fPRQB) 

INTEGER  INPGf Jf JEf JRfH 
REAL  RLfUHfUC 

DOUBLEPRECISION  PET ( 0  J  *  > f AIR ( 0 1*  >  f PROB . COMB f XTOP  f P J 
PR0B=0. 

DO  101  J=JE » INPG 

PJ=<4.*M*<RL**3. ) >/<3.*VH*UH*UC*<PET<J>**3. > ) 

IF (PJ.GT .1.0)  PJ=i.O 

PROB=PROB+COMB ( J  t JE) *XTOP ( 1 . -P J f  JE ) *XTOP  <  P J  f  J- JE ) *AIR ( J ) 

CONTINUE 

END 

SUBROUTINE  REFUEL <  JR f JE t INTA f INTGf TESC f PREF f PROB > 

INTEGER  Jf JEf JRfINTAfINTG 
REAL  PREF 

DOUBLEPRECISION  TESCCOJ*) fXPREFfCOMBfXTOPfPROB 

XPREF=PREF 

IF< JE.LE. INTA) 

J  THEN 

PROB=COMB <  JE  f  JR ) *XTOP  <  XPREF f  JR ) *XTOP ( 1 . -XPREF f  JE- JR  > 

ELSE 

SUMF=0 . 

DO  102  J=JE-INTAfINTG 
SUMF =SUMF +TESC  <  J ) 

CONTINUE 
SUMS=0 » 

DO  202  J=MAX<JR-INTAfO)fJE-INTA-1 
SUMS=SUMS+COMB  <  J+ 1 NTA f  JR ) *XTOP  <  XPREF  f  JR ) 

X  *XTOP  < 1 » -XPREF  f  J+ INTA- JR ) *TESC ( J ) 

CONTINUE 

PROB-COMB ( JE  f  JR ) *XTOP  <  XPREF  f  JR ) *XTOP ( 1 . -XPREF  t  JE- JR ) 
J  •  *SUMF+SUMS 

ENDIF 
END 


4180=  SUBROUTINE  FAD  < I » JE  >  JR  >  C  > DUMMY  f  F  t  PDC  t  F’K  i  F'ROB ) 

4190=  INTEGER  I t J » JE > JR > C » F t DUMMY , BA » BB > N 

4200=  REAL  PDC»PK 

4210=  DOUBLEPRECISION  PSA>PSB»PROB 

4220=  DOUBLEPRECISION  COMB»XTOP 

4230=  DOUBLEPRECISION  P»X*Y 

4240=  IF ( JE+C . EG . 0 )  GO  TO  1303 

4250=  N=0 

4260=  X=1.-PDC*PK 

4270=203  CONTINUE 

4280=  IF ( F ► LT . JE+C )  GO  TO  103 

4290=  F =F- JE-C 

4300=  N=N+ 1 


4310= 

4320=103 

4330= 

4340= 

4350= 

4360= 

4370= 

4380= 

4390= 

4400= 

4410= 

R-I ) 

4420= 

4430=303 


GO  TO  203 

CONTINUE 

IF (DUMMY. EG. 0) 

THEN 

IF<I.GT. JE-JR) 

THEN 

PR0B=0. 

ELSE 

IF(F.LT. JE-JR)  GO  TO  303 
PSA=XTOP(X  »N+1 ) 

PROB=CQMB <  JE-JR  f I ) *XTOP ( PSA  1 1 ) *XTOP ( 1 . -PSA  >  JE- 

GO  TO  403 
CONTINUE 


4440= 

4450= 

4460= 

4470= 

4480= 

4490= 

4500=403 

4510= 

4520= 

4530= 

4540= 

4550= 

4560= 

4570= 

4580= 

4590= 

4600= 

4610= 

4620=603 

4630= 

4640= 


BA=JE— JR-F 
BB=F 

PSA=XTOP ( X  *  N ) 

PSB=XTOP( X  t N+l ) 

P=BA*PSA/ ( JE-JR ) +BB*PSB/ <  JE-JR ) 

PROB=COMB ( JE-JR  r I >  *XTOP ( P  r I ) *XTOP ( 1 . -P » JE- JR-I 
CONTINUE 

ENDIF 

ELSEIF<  DUMMY. EG . 1 ) 

THEN 

IF<I.GT.JR> 

THEN 

PR0B=0. 

ELSE 

IF<F.GT. JE-JR)  GO  TO  603 
PSA=XTOP<XrN) 

PROB=COMB <  JR » I ) *XTOP ( PSA » I ) *XTOP  < 1 . -PSA  * JR-I ) 
GO  TO  703 
-  CONTINUE 

IF(F.GT.JE)  GO  TO  903 
BA=F-JE+JR 


4650= 

4660= 

4670= 

4680= 

4690= 

4700= 


BB=JE-F 

PSA=XTOP  <  X  f N+l ) 

PSB=XTOP  <  X  r  N ) 

P=BA*PSA/JR+BB*P3B/JR 

PROB=COMB ( JR » I ) *XTOP <P» I ) *XTOP ( 1 . -Pf  JR-I ) 
GO  TO  703 
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4710=903 

4720= 

4730= 

4740=703 

4750= 

4760= 

4770= 

4780= 

4790= 

4800= 

4810= 

4820= 

4830= 

4840= 

4850= 

4860= 

4870= 

4880= 

4890= 

4900=1303 

4910= 

4920= 

4930=1403 

4940= 

4950= 

4960= 

4970= 

4980= 

4990= 

5000= 

5010= 

5020= 

5030= 

5040= 

5050= 

5060= 

5070= 

5080= 

5090= 

5100= 

5110= 

5120= 

5130= 

5140=104 

5150= 

5160= 

5170= 

5180= 

5190= 

5200= 

5210= 

5220=*E0R 


CONTINUE 
PSA=XT0P(X,N+1 ) 

PROB=COMB  <  JR » I ) *XTOPCPSA» I )*XTOP< 1 .-PSA> JR- 1 ) 
CONTINUE 

END  IF 
ELSE 

IF < F .LE  *  JE ) 
t  THEN 

PROB=XTOP ( X » N ) 

ELSE 
Y=F- JE 
PSA=Y/C 
Y= JE-F+C 
PSB=Y/C 

PROB=PSA*XTQP  <  X , N  +  l ) +PSB*XTQP <  X » N ) 

END  IF 

ENDIF 

F=F+N*< JE+C) 

GO  TO  1403 

CONTINUE 

PR0B=0. 

IF(I.EQ.O)  PR0B=1 . 0 

CONTINUE 

END 

SUBROUTINE  BOMBS <  N»PK » WL>  Wf  D»  Rr  NSAM t PDC r PF »L » A*  VB»PTSAM>  PRO 
INTEGER  J,N,NSAM 

REAL  R(0:25) »PTSAM»PK»WL»U»D»PF»A»L»UB»PDC 
DOUBLEPRECISION  PA ,  PB .  F'C * PD » PE » PROS  -  SUM 
DOUBLEPRECISION  X» Y»Z» XPTSAM * XTOP 
XPTSAM=PTSAM 
SUM=0 . 

DO  104  J=1 »N-1 

X=1 . -PK*UL*R ( J ) / <  U*D ) 

PA=XTOP<X,NSAM) 

PB= <  FDC*PF*L/ <  PDC*PF*L-A ) ) *EXP  < -A*R  <  J  > /VB ) 

PC=  <  A/ <  PDC*PF*L-A ) ) *EXP  < -PDC*PF#L*R  <  J  > /VB ) 

PD=PA*  <  PB— PC ) 

Y=1 .-PK*WL*R< J+l >/<W*D> 

PA=XTOP( Y  fNSAM) 

PB= ( PDC*PF*L/ ( PDC*PF*L-A ) >  *EXP  < -A*R ( J+l  > /VB ) 

PC=  ( A/  <  PDC*F’F*L-A  > )  *EXP  < -PDC*PF*L*R ( J+l )  /UB ) 

PE=PA*<PB-PC> 

SUM=SUM+ J*XTOP  ( XPTSAM  >  J )  *  <  F'D-PE*XPTSAM ) 

CONTINUE 

Z=1.-PK*UL*R<N>/<W*0> 

PA=XTOP<Z,NSAM) 

PB=  <  PDC*PF*L/ <  PDC*PF*L-A ) >  *EXP  < -A*R ( N  > /VB  > 
PC=<A/<PDC*PF*L-A> >#EXP<-PDC*PF*L*R<N)/VB) 

SUM=SUM+N*XTOP ( XPTSAM  r N  >  *PA* ( PB-PC ) 

PROB=SUM/N 

END 
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APPENDIX  B 


ESTIMATION  OF  THE  TIME-DEPENDENT 
INTERCEPT  DENSITY  FUNCTION 


This  appendix  describes  the  technique  developed  by  R.  H.  Monahan 
of  the  Stanford  Research  Instititue  for  estimating  x^,  the  time 
dependent  intercept  density  function,  for  the  COPEM-1  bomber  penetra¬ 
tion  model.  Time  constraints  and  model  scope  resulted  in  this  paper 
assuming  a  constant  intercept  parameter;  therefore,  the  procedure 
described  below  is  not  meant  to  be  directly  applicable  to  this  model. 

The  procedure  is  presented  because  it  gives  some  insight  into  the 
factors  that  should  be  considered  when  selecting  a  value  of  X.  It  may 
suggest  ideas  for  an  approach  to  the  development  of  a  similar  procedure 
for  estimating  a  time-independent  intercept  parameter.  It  may  also 
suggest  a  means  for  increasing  the  realism  of  this  model  by  extending 
it  to  the  case  of  a  time-dependent  intercept  parameter. 

In  reviewing  Monahan's  derivation,  one  minor  conceptual  error  was 
discovered.  In  this  appendix,  the  method  for  estimating  xt  is  presented 
as  derived  by  Monahan;  the  error  is  pointed  out  and  both  Monahan's  and 
the  corrected  results  are  given. 

ESTIMATION  OF  Xt 

Recall  that  for  a  Poisson  distribution  with  parameter  X,  the 
waiting  time  to  the  first  event  is  distributed  exponentially  with  the 
same  parameter  X.  Recall  also  that  the  expected  value  of  the  exponen¬ 
tial  is  l/x.  Therefore,  if  an  estimate  can  be  found  of  the  expected 
waiting  time  to  the  first  intercept  of  a  bomber  detected  at  time  t,  the 


159 


reciprocal  of  this  will  in  turn  be  an  estimate  of  the  parameter  >t. 

To  find  an  estimate  of  the  waiting  time  to  the  first  intercept, 
the  penetration  corridor  is  divided  into  r  rows  and  c  columns  for  a 
total  number  I=rc  cells  which  are  labelled  sequentially  from  left  to 
right  and  top  to  bottom  from  1  to  I  as  in  Figure  B.l. 


□1 

a 

a 

H 

BI 

a 

□1 

H 

BI 

m 

BI 

BI 

n 

B 

B! 

BI 

B 

■ 

BI 

BI 

■: 

BI 

■ 

■ 

■ 

B 

fli 

■ 

B 

<9 

BI 

B 

B 

m 

■ 

1 

B 

m 

_ 

B 

— 

■ 

B 

r~ 

_ 

C 

■ 

■ 

B 

B 

B 

B 

g 

■ 

■ 

B 

B 

B 

B 

■ 

■ 

■ 

■ 

B 

_ 

_ 

_ 

5 

■ 

■ 

_ 

_ 

■ 

■ 

r 

■ 

| 

■ 

IB 

B 

B 

IB 

■ 

IB 

_ — 

■ 

_ 

IB 

IB 

B 

IB 

□ 

IQ 

>81 

i -  V/  - » 


s 


FIGURE  B.l 

Penetration  corridor  partitioned  for  calculation  of 


Define  the  following  variables: 

n^  -  the  total  number  of  fighter-interceptors  at  airbase  k; 
n^(t)  -  expected  number  of  interceptors  from  base  k  that  are  air¬ 
borne  and  unassigned  at  time  t; 

p.k  -  probability  that  an  arbitrary  airborne  interceptor  from  base 
k  is  in  cell  i ; 
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M 


Pik(t)  -  probability  that  there  is  an  interceptor  from  base  k 
in  cell  i  at  time  t; 

q.(t)  -  probability  that  there  are  no  interceptors  available  to 

J 

be  assigned  to  a  bomber  in  cell  j  at  time  t; 

q.j(t)  -  probability  that  an  interceptor  assigned  to  a  bomber  in 
cell  j  at  time  t  is  from  cell  i; 

Q.  .(t)  -  probability  that  there  is  an  interceptor  in  cell  i  that 
’  j 

can  be  assigned  to  a  bomber  in  cell  j  at  time  t; 

R^  -  ranking  of  cell  i  relative  to  cell  j  in  terms  of  the  time  it 
would  take  a  fighter  from  cell  i  to  intercept  a  bomber  first  detected 
in  cell  j ; 

Tc  -  bomber  flight  time  between  cells; 

T.  .  -  expected  time  to  the  first  intercept  given  that  an  inter- 

*  J 

ceptor  from  cell  i  is  assigned  to  a  bomber  in  cell  j; 


indicator  variable  equalling  1  if  an  interceptor  from  air 


base  k  in  cell  i  is  capable  of  intercepting  a  bomber  in  cell  j  and  0 


if  not. 


Both  bombers  and  interceptors  are  assumed  to  always  be  in  the 
center  of  whatever  cell  they  are  occupying.  Bombers  fly  straight  down 
the  corridor  at  a  constant  velocity  while  interceptors  can  fly  in  any 
direction  at  a  constant  velocity  greater  than  the  bombers  velocity. 

With  these  assumptions,  the  calculation  of  T. .,  the  time  it  takes  an 

*  J 

interceptor  from  cell  i  to  intercept  a  bomber  detected  in  cell  j,  can 
be  calculated  with  basic  algebra  and  trigonometry.  This  will  not  be 
done  here.  These  T^  can  be  used  to  define  the  R^.:  for  a  specific  j, 
will  equal  1  for  the  i  which  has  the  minimum  T^.,  2  for  the  i  having 
the  next  least  T^,  and  so  on. 


161 


How  consider  the  P.^it).  For  exactly  ( t )  interceptors  from 
base  k  airborne  and  unassigned  at  time  t,  one  would  have 


P* 

i  k 


(t) 


1-(1-Pik)Nk(t) 


However,  it  is  assumed  that  only  E['l^(t)]  =  n^(t)  is  known  so  one 
attempts  to  find 

Pik(t)  =  E[PJk(t)]  =  E[1-{1'pik)N|<(t)]  3-1 

If  N. (t)  is  assumed  to  be  binomially  distributed  with  parameters 
P=n^(t)/nk  and  N=nk,  then  it  can  be  shown  that 

Pik(t)  =  l-[l-Pfknk(t)/nkink  8.2 


With  this  result,  one  can  now  find  the  probability  that  there  is 
an  interceptor  in  cell  i  that  can  be  assigned  to  a  bomber  detected  in 
cell  j  at  time  t; 

Qij(t>  ■  i-k:1[1-*ukpik(t>i*  B-3 

where  K  is  the  total  number  of  interceptor  bases. 

To  continue,  let  i  denote  the  value  of  i  for  which  R.  .=m  for  a 
m  i  j 

specific  j.  One  can  now  find  the  probability  that  the  interceptor 
assigned  to  a  bomber  in  cell  j  is  from  cell  i  (Note:  it  is  possible 
that  more  than  one  interceptor  will  be  assigned  to  a  bomber  in  cell  j 
at  time  t  but  one  is  only  concerned  here  with  the  interceptor  which 
will  make  the  earliest  interception  of  the  bomber.) 
m-1 

q1  At)  =  [1-  E  q.  .(t)}Q.  .(t),  with  q.  ,(t)-Q.  (t)  B.4 

mJ  x=l  xJ  ’mJ  ’l0 
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From  this,  the  probability  that  there  are  no  interceptors  available 
to  be  assigned  to  a  bomber  in  cell  j  at  time  t  is 

I 

q •  ( t )  =1-1  q,  ,(t)  3.5 

J  x«l  V 

Now  consider  the  expected  time  to  the  first  intercept  of  a  bomber 
detected  in  cell  j  at  time  t..  This  value  may  be  broken  into  two 
parts  and  expressed 

E[Tj(t)]  =  E: [Tj(t)]  +  E2£TJ(t) ]  9.6 

where  one  has  conditioned  on  the  availability  or  unavailability  of  an 
interceptor  to  be  assigned  to  the  bomber.  The  first  term  represents 
the  expected  value  of  the  time  to  be  first  intercept  given  that  there 
is  an  interceptor  available  to  be  assigned  to  the  bomber  times  the 
probability  that  there  is  an  interceptor  available.  The  second  term 
represents  the  expected  time  to  the  first  intercept  given  that  there 
are  no  interceptors  immediately  available  to  be  assigned  to  the 
bomber  times  the  probability  that  there  are  no  interceptors  availabH-. 
If  A  is  the  event  that  there  is  an  interceptor  available  to  be 
assigned  to  a  bomber  detected  in  cell  j  at  time  t,  then  the  above 
can  be  expressed  mathematically  as 

E1[Tj(t)]  =  E[Tj(t) | A] P ( A )  and  B.7a 

E2[Tj(t)l  =  E[Tj(t)|AC]P(AC)  B  .7b 

(Note:  the  E1  and  E2  notation  may  seem  a  litte  confusing  but  is  being 
used  here  because  that  is  how  the  derivation  is  presented  in  Ref  5.) 
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The  expected  value  of  the  time  to  the  first  intercept  if  there 
is  an  interceptor  immediately  available  can  be  expressed  in  terms  of 
quantities  already  expressed, 

El[Tj(t)I  ■  B'3 

It  remains  to  find  the  expected  value  of  the  time  to  the  first  inter¬ 
cept  if  there  are  no  interceptors  available  to  be  assigned  to  the 
bomber  detected  in  cell  j  at  time  t. 

To  approximate  this  quantity,  Monachan  suggests  the  following 
appraoch.  Allow  the  bomber  to  continue  to  the  next  cell,  consuming  a 
time  T  .  Attempt  to  assign  an  interceptor.  If  there  are  still  none 
available,  allow  the  bomber  to  continue  to  the  next  cell,  repeating 
until  a  time  is  reached  when  an  interceptor  is  available  to  be 
assigned  to  the  bomber.  To  find  the  expected  time  to  intercept  if 
there  are  no  interceptors  available  at  t,  one  can  condition  on  the 
number  of  cells  the  bomber  flies  before  an  interceptor  becomes  avail¬ 
able  to  be  assigned  to  it.  The  following  expression  is  proposed: 

E2tTj{t)]  =  qj(t)[l-qj(t)HTc+E1(Tj(t))] 

+  qjCt)2[l-qj(t)][2Tc+E1(Tj(t))] 

+  . . . 

+  qj(t)m[l-qj(t)][mTc+E1(Tj(t))] 

+  ...  '  .  8.8 

Several  things  should  be  noted.  First,  both  the  T^  and  the  q^ 
may  be  functions  of  time.  Furthermore,  each  time  the  bomber  moves 
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to  a  new  cell,  the  value  of  j  is  incremented  by  c,  the  number  of 
columns.  With  these  two  observations,  the  second  term  in  the  above 
sum  should  look  1  ike 

and  the  other  terms  would  be  modified  similarly.  Monahan's  only  com¬ 
ment  on  the  approximation  is  that  he  "assumes  a  semblance  of  station- 
arity  (Ref  5:46) 

The  second  thing  to  be  noted  is  that  there  are  only  a  finite 
number  of  cells  and  if  the  bomber  penetrated  beyond  the  end  of  the 
corridor  he  would  presumably  be  safe  from  interception.  The  above 
sum,  therefore  should  be  finite,  rather  than  infinite.  Mo  comment 
is  made  on  this. 

There  is  nothing  really  wrong  with  these  two  approximations  as 
long  as  they  are  recognized  as  such.  However,  even  if  one  accepts 
the  approximations,  there  is  an  error  in  the  sum  given  in  Equation 
B.8. 

Recall  that,  using  A  as  previously  defined, 

E1£Tj(t)]  «  E [Tj ( t ) j  A]P(A)  B  .10 

Rearranging  the  above  expression  and  noting  that  P(A)  can  be  expressed 
as  l-q.(t) ,  one  has 

J 

ElTj(t)  |  A]  =  E:1ITj(t)]/(l-qj(t))  B.ll 

Now  consider  what  Is  happening  when  is  conditioned  over  the 
number  of  cells  the  bomber  flies  before  an  Interceptor  becomes  avail- 


able  to  be  assigned  to  it.  The  m'"  term  of  the  sum  gives  the 
expected  waiting  time  for  the  first  intercept  given  that  an  inter¬ 
ceptor  has  become  available  after  the  bomber  has  flown  exactly  m  cells 
times  the  probability  that  the  first  interceptor  becomes  available  at 
that  time.  The  probability  is  geometric  and  is  given  by  the  expression 

P(Sm)  =  qj(t)fn[l-qj(t)]  3.12 

where  the  approximations  discussed  above  are  being  used. 

The  expected  waiting  time  given  that  an  interceptor  becomes  avail¬ 
able  after  the  bomber  flown  m  cells  is  just  the  time  it  took  the 
bomber  to  fly  the  m  ce.’s  .0  us  the  expected  time  to  the  first  inter¬ 
ceptor  given  that  there  is  an  interceptor  availabl e  to  be  assigned. 


E(Wm)  *  mTc  +  E[TJ(t){A]  B  .13 

(In  the  above  expressions,  Bm  is  the  event  that  the  first  inter¬ 
ceptor  to  be  assigned  to  the  bomber  becomes  available  after  the  bomber 
has  flown  exactly  m  cells  and  Wm  is  the  waiting  time  to  the  first 
intercept  given  B^.) 

The  mth  term  of  the  sum  for  E2  should  therefore  read: 
mtfl  term  *  P(Bm)E(Wm) 

■  qj(t)m[l-q4(t>] [mTc  +  E(Tj(t) ( A] 

»  qJ(t)m[l-qj(t) 1 CmTc  +  E1(Tj(t))/(l-qj(t))]  3.14 

This  differs  from  the  expression  in  Equation  B.8  by  the  factor  of 
l-qj(t)  in  the  denominator  of  Ej. 
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The  sum  can  easily  be  reduced  by  using  the  identities 

oo  . 

Z  x1  =  x/ ( 1 -x )  and 

i=l 

I  ix1  =  x/(l-x)2 

i=l 

Monahan  finds 

E2[Tj(t))  =  qJ(t)E1[Tj.(t)  ]  +  TcqJ(t)/[l-qJ( t)  ]  B  .15 

which  can  be  added  to  to  give 

E[Tj(t)l  =  [l+qj(t) ]Ej [Tj ( t } ]  +  Tcqj(t)/[l-qj(t)]  B.16 

The  corrected  expressions  are 

E2[Tj(t)]  =  qj(t)E1[Tj(t)]/[l-qj(t)]  +  Tcq J ( t ) / [ 1 - q ^ ( t ) ]  3.17 

and 

E[Tj(t)l  =  E1tT.(t)]/tl-qJ(t)]  +  Tcqj(t)/[l-qj(t)]  3.18 

Because  of  the  assumption  of  "semblance  of  stationarity,"  an 
easy  method  exists  to  check  the  above  expression.  If  one  conditions 
the  expected  time  to  first  intercept  on  the  availability  or  unavail¬ 
ability  of  an  interceptor  to  be  assigned  to  the  bomber,  one  has 

E[Tj(t) ]  =  Ej  £T j ( t ) ]  +  qj[Tc  +  E(Tj(t))]  8.19 

This  is  a  simple  linear  equation  E[T.(t)]  which  can  be  solved  to  give 

J 

E[Tj(t)]  =  E1lT1U)3/[l-qJ(t)]  +  Tcqj(t)/[l-qj(t)l  B.20 
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as  above. 


To  conclude  the  original  problem,  which  was  to  estimate  the 
parameter  one  has 


t  - 

H.  ( T)  =  z  P.(t)E[T.(t)] 
r  j-1  3  3 


B .  21 


where  I  is  the  time  to  the  first  intercept  and  Pj(t)  is  the  probability 
that  a  bomber  detected  at  time  t  is  in  cell  j.  The  reciprocal  of 
E .  ( t )  is  the  desired  estimate  of  \t. 

ESTIMATION  OF  nfe(t) 

The  next  problem  is  to  find  a  suitable  estimate  for  nk(t),  the 
expected  number  of  fighters  from  airbase  k  available  for  assignment 
as  a  function  of  time.  Define  the  following  variables: 


At  -  interval  between  successive  estimations  of  X^.; 

Tj  -  ith  time  interval  for  the  calculation  of  Xt;  T.  is  the 
time  interval  [(i-l)At,  iAt]; 

nfc  -  total  number  of  fighters  stationed  at  base  k; 

)  -  the  number  of  bombers  detected  in  the  time  interval  T. ; 
N ^ ( T i )  -  the  number  of  fighters  from  base  k  assigned  to  bombers 
detected  in  the  time  interval  ; 

P k ( T i )  -  the  probability  that  a  bomber  detected  in  the  time 
interval  will  have  a  fighter  from  base  k  assigned  to  it; 

-  total  time  a  fighter  can  remain  airborne  before  recycling; 
Tr  -  the  minimum  time  for  a  fighter  to  recycle; 
r$  -  the  scramble  rate  at  which  fighters  take  off,  measured  in 
fighters/unit  time. 
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The  approach  is  to  first  construct  a  fighter  availability  curve 
for  base  k  assuming  that  no  bombers  are  detected  or  intercepted.  The 


curve  is  then  iteratively  modified  at  discrete  time  intervals  to 

account  for  the  change  in  fighter  availability  arising  from  the 

detection  and  interception  of  bombers  in  each  interval. 

Assume  that  before  the  battle,  the  number  of  airborne  fighters 

is  maintained  at  some  constant  level,  n^(0).  The  maximum  possible 

★ 

value  for  that  constant,  n^O),  is 

"k(0)  ’  rmr  "k  B-22 

t  r 

At  time  t=0,  all  the  fighters  are  flushed  and  the  number  of  air¬ 
borne  fighters  increases  linearly  according  to  the  scramble  rate,  r  . 
All  the  fighters  will  be  airborne  at  time  t^ , 


At  some  time  airborne  aircraft  will  begin  to  recycle  and  the 
number  of  airborne  fighters  will  now  decrease  according  to  the 
scramble  rate  for  a  length  of  time  T  ,  the  time  taken  by  one  fighter 
to  recycle.  At  this  time,  the  first  fighters  to  return  for  recycling 
will  be  taking  off  again  and  the  airborne  force  level  will  remain 
constant  until  the  last  fighter  to  be  recycled  has  landed.  The  air¬ 
borne  force  level  will  then  begin  to  increase  again  according  to  the 
scramble  rate  and  the  process  will  be  repeated.  The  process  is 
illustrated  in  Figure  B.2  for  n(((0)=20,  n^O,  rs=l/minute,  T^-50 
minutes,  and  Tr*20  minutes. 
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Figure  8.2  gives  the  number  of  airborne  fighters  from  base  k  as 

a  function  of  time.  To  find  the  number  of  airborne  fighters  available 

for  assignment,  it  is  necessary  to  subtract  off  the  number  of  fighters 

assigned  as  a  function  of  time.  The  method  for  doing  this  will  be 

demonstrated  by  continuing  the  example  presented  in  Figure  B.2.  For 

this  example,  esimates  (a,  )  will  be  computed  afor  the  intercept 

i 

density  function  A^,  at  10-minute  intervals.  An  estimate  (user 
specified  of  Afl  is  assumed  to  be  such  that  Aq"*,  the  expected  time 
to  first  intercept,  is  14  minutes. 
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Using  a  convolution  of  the  detection  time  and  arrival  time  dis¬ 
tributions  (Pj(t^)  and  p  (t  )),  an  expected  number  of  bombers 
detected,  denoted  N^(T^),  can  be  found  for  the  first  19-minute  inter¬ 
val.  Say  this  is  7  bombers.  It  is  assumed  that  when  making  an 
original  estimate  of  xn,  an  unconditional  probability,  ?,(!,),  that  a 
bomber  detected  in  this  time  interval  will  have  a  fighter  from  air¬ 
base  k  assigned  to  it  can  be  calculated.  For  this  example,  this 
probaoility  is  assumed  to  be  5/7.  From  this,  one  finds  the  expected 
number  of  fighters  from  airbase  k  assigned  to  attack  bombers  in  this 
interval  is 

MV  ■  MVMV  B-24 

which  in  this  example  is  5. 

Thus  at  time  t=I0  minutes,  the  number  of  fighters  from  base  k 
available  for  intercept  has  decreased  by  5.  This  point  is  labelled 
A  in  Figure  B.3.  It  is  assumed  for  this  example  that  the  probability 
of  kill,  given  an  intercept,  for  a  fighter  against  a  bomber  is  0.7. 

The  expected  time  to  a  successful  intercept  (when  the  bomber  is  killed) 
is,  therefore,  equal  to  14/0.7,  or  20  minutes.  Thus  at  time  t=20 
minutes,  all  the  fighters  assigned  to  bombers  in  the  first  10-minute 
interval  will  rejoin  the  unassigned  fighter  force.  This  is  indicated 
by  the  line  segment  connecting  points  B  and  C  in  Figure  B.3. 

This  completes  the  modification  to  the  availability  curve  based 
on  the  activity  during  the  first  10-minute  interval;  the  remainder 
of  the  curve  remains  unchanged  from  Figure  B.2.  At  this  point,  xT 

i 

can  be  calculated. 
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FIGURE  B . 3 

Fighter  availability  curve  modified  for  activity  in  T^ . 

This  procedure  is  repeated  for  T2»  T3,  ....  Tn*  At  each  stage, 

the  previous  estimate  of  the  intercept  density  function  (,\T  )  is 

1-1 

used  to  find  the  expected  time  to  a  successful  intercept.  The  avail¬ 
ability  curve  can  then  be  modified  for  the  activity  in  that  time 
interval  and  a  value  for  Ay  can  be  calculated. 

In  performing  the  above  calculations,  it  has  been  assumed  that 

XT  was  constant  over  each  interval.  Specifically,  A.  =  at  for 
'i  th  Z  1-1 

the  entire  ic  interval.  However,  once  the  calculations  are  complete, 

and  estimates  have  been  made  for  all  the  AT  ,  values  for  at 

i  1 

interior  points  on  the  intervals  may  be  found  by  interpolation  if 
desired.  This  completes  the  derivation. 
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